AMO  84t>B5 


'•  *.Vr"  H  ,4  '  f  J  '  J*:" 

■  *.  Y.  •  •  \  ■  '/"  jf  )% 

■S  'vf  V  '  fe  •'■-7 1 

\=,  6 ■• /'-L.  •  ;  0  ,£0 

TECHNICAL  REPORT  GL-liO  2  '  -'fS^ 


VMD  GRAPHICAL  METHODS  FOR 
OF  SLOPES  IN  ROCK  MASSES 


f'L  Or.,  c.  vi.  Co  rainy,  *nd  /\.  K.  Aiyer 
Uepnrtmenfc  of  Civil  Engine*:, lag 
Versify  of  Illinois,  Urban?.,  HI.  61801 


March  1950 

inf.  of  NCG  Technical  Report  No.  36 


i'iov.-i  Foi  Public  !?cl-:.T>e;  Uktribution  (Min 


'MS*r  ’V'V.Cv’  . 


,  X  '  -i  •  f  A’  >  t 

i  V  .  L5-6v  't 6 ; .  ■ 

0  ■'\>v,vV  ...  '  ,  »'•  ' 

rXV  ,  -  ...  ,  w- 


.  n  .  -  •*  -J  1  y 

i,  V  -  '  •  ,*  •  *  V”  "*  '-<«>’  *  AT  J  -  -  - 

S3  .  -  4  1  ?!  4  f  1  S*  8J 


*11  J  .  S,  4  >if,  Jk  M  h  ,1 

-4.  -  S*c  :>  44  V  .  9  ‘  • 

*ki«.  4  ,/  ~  "4  ' 

:  ...<■;>.  ■pioL.iL^iSid.pS;  :LV  .  •  1  . 


X  ''X  " 
-fV 


Army  Engineer  Nuclear  Coloring  Group 
Livermore,  Calif.  94550 

L  nkracfc  No.  DACW39-67  C-0097 


nit,  ,  id  by  Geotechnical  Labor*lr>iy 
Engineer  Waterways  Experiment  Station 
.  Ho;<  631,  Vicksburg,  Miss.  89180 


O  f 


£>.  rv 

'*■  .*  J? 


i\”  V1*,  ;*A 

%  ..  • 


ii  ropt.rt  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


in  this  report  are  not  to  be  construed  as  an  official 
C*nt  of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 


SECURITY  CLASSIFICATION  of  THIS  FAQE  (When  Dot*  Entered) 


REPORT  DOCUMENTATION  PAGE 


^ANALYTICAL  AND  GRAPHICAL  JjETHODS  FOR  THE 
"^ANALYSIS  OF  JSLO'PES  IN  JOCK_MASSES  .  - - 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


s.  type  op  report  a  period  covered 

Reprint  of  NCG  Technical 
Report  No.  36 


*■  PERFORMING  ORO.  REPORT  NUMBER 


~r&s?pr  1  ” '  ■"  I  "  —  I 

c  *■■**/ Hendron,  Jr.  E.  J./ Cording j  and  /* 
A.  K.  /Aiyer  f" '  1  f  .  — *  (j 

S.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Department  of  Civil  Engineering  ^ 

University  of  Illinois,  Urbana,  Ill.  6l801 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ~ 

U.  S.  Army  Engineer  Nuclear  Cratering  Group  f 
Livermore,  Calif.  9^550  \ 


t.  CONTRACT  OR  GRANT  NUMBERS 


;5>  J  DACW  39-67-0-^97/ 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


«.  MONITORING  AGENCY  NAME  A  AODRESSfff  dttlmrmnt  from  CtmlnWnl  Otllem)  IS.  SECURITY  CLASS,  (ol  Mm  import) 

U.  S.  Army  Engineer  Waterways  Experiment  Station  . 

Geotechnical  Laboratory  Unclassified 

P.  0.  Box  631,  Vicksburg,  Miss.  39180  ts«.  declassification/’ 00 wngrai 


s«.  declassification/oownsraoing 
SCHEDULE 


I  IS.  DISTRIBUTION  STATEMENT  (ol  Mm  Rmporl) 


Approved  for  public  release;  distribution  unlimited. 


1 X7 


M7.  DISTRIBUTION  STATEMENT  (ot  th«  obotroct  onfredlk 


>-  a-.  rPf<  - 


WE 


Dynamic  slope  stability 
Graphical  methods 
Rock  masses 
Rock  mechanics 


Slope  stability 
Stereonet 
Vector  analysis 


ACT  fCifftwiM  mm 


1  aM  H  nimmy  Ml  identity  by  block  number) 


In  this  report  the  methods  of  analyzing  the  static  stability  of  rock 
slopes  cut  by  a  three-dimensional  network  of  discontinuities  are  given.  The 
general  use  of  vector  analysis  to  solve  these  problems  analytically  is  de¬ 
scribed  and  a  method  utilizing  stereonets  to  solve  these  problems  graphically 
is  also  given.  For  both  the  graphical  and  analytical  methods  the  general 
analysis  of  slopes  cut  by  one,  two,  or  three  sets  of  discontinuities  is  pre¬ 
sented  which  can  take  into  account  the  porepressures  acting  on  the  ,  . 


DO  /A 


EDITION  OF  >  MOV  At  IS  OBSOLETE 


_ Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  ( 


1  Dmto  Bntorod) 


20.  ABSTRACT  (Continued). 

discontinuities  and  external  forces  acting  on  the  slope.  Detailed  examples 
are  given  to  illustrate  both  the  graphical  and  vector  methods  of  analysis. 

The  dynamic  stability  of  slopes  is  also  treated  in  this  report.  It  is 
shown  that  the  dynamic  resistance  of  a  three-dimensional  rock  slope  can  be 
calculated  by  either  the  graphic-stereonet  method  or  the  analytic  vector 
analysis  method.  The  dynamic  resistance  can  then  be  used  to  estimate  the 
movement  of  the  slope  under  dynamic  loading  using  a  procedure  given  by 
Ifewmark  (1965).  A  criterion  is  then  given  for  determining  if  the  calculated 
movement  of  the  rock  slope  is  acceptable  or  harmful. 


Unclassified 

SCCUNITV  CLASSIFICATION  OF  THIS  FAOCOWlF" 


FOREWORD 


This  study  was  performed  by  the  Department  of  Civil  Engineering, 
University  of  Illinois,  un^T  contract  to  the  U.  S.  Army  Engineer  Water* 
ways  Experiment  Station,  Vicksburg,  Mississippi  (WES)  for  the  U.  S.  Army 
Engineer  Nuclear  Cratering  Group,  Lawrence  Radiation  Laboratory,  Livermore, 
California  (NCG).  The  study  was  performed  under  Contract  No.  QACW  39-67- 
C-0097,  "Evaluation  of  Analytical  Methods  of  Determining  the  Stability  of 
Rock  Slopes,"  and  was  funded  by  NCG  IA0  2-63,  "Engineering  Properties  of 
Nuclear  Craters."  The  contract  was  negotiated  on  16  June  1967. 

This  report  was  prepared  by  Messrs.  Hendron,  Cording  and  Aiyer;  and, 
was  reviewed  by  Don  C.  Banks,  Chief,  Rock  Mechanics  Section,  Soils  Division, 
WES  and  by  Major  Richard  H.  Gates,  C.E.,  Chief,  Engineering  Geology  Division, 


The  contract  was  monitored  by  Don  C.  Banks,  WES.  The  Contracting  Officer 
at  the  time  of  publication  was  COL  Ernest  D.  Peixotto,  CE,  Director  of  WES. 
Technical  Director  of  WES  was  Mr.  Fred  R.  Brown.  The  Director  of  NCG 
was  LTC  Robert  L.  LaFrenz,  CE,  and  the  Technical  Director  was  Mr.  Walter  C. 


Accession  For 

i  itTiij  I 

I  P.O  I.  B 
i  v  ... -.crd 


■ r. 


ABSTRACT 
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is  also  given.  For  both  the  graphical  and  analytical  methods  the  general 
analysis  of  slopes  cut  by  one,  two,  or  three  sets  of  discontinuities  is 
presented  which  can  take  into  account  the  porepressures  acting  on  the  dis¬ 
continuities  and  external  forces  acting  on  the  slope.  Detailed  examples  are 
given  to  illustrate  both  the  graphical  and  vector  methods  of  analysis. 
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is  shown  that  the  dynamic  resistance  of  a  three-dimensional  rock  slope 
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CHAPTER  ONE 


INTRODUCTION 


1.1  General 

The  design  and  analysis  of  rock  slopes  is  somewhat  different  than 
the  design  and  analysis  of  slopes  in  soil  because  of  the  patterns  of 
discontinuities  in  the  rock  mass.  The  spatial  orientation  of  these 
discontinuities  and  the  shearing  resistance  along  them  govern  the 
stability  of  rock  slopes.  Thus  the  method  of  analysis  used  must  take 
into  account  the  three  dimensional  intersection  of  the  joint  sets  with 
each  other  and  intersection  of  these  discontinuities  with  the  face  or 
surface  of  the  rock  slope.  Limit  equilibrium  methods  of  analysis  have 
recently  been  developed  to  analyze  these  problems  in  three  dimensions 
which  will  be  explained  and  illustrated  in  this  report. 

In  all  methods  of  limit  equilibrium  analysis  the  shape  of  the 
potential  failure  is  assumed  at  the  outset.  In  the  limit  equilibrium 
methods  used  for  soil  slopes,  sections  of  log  spirals  or  circles  are 
normally  chosen  to  represent  the  failure  surface.  Although  displacements 
are  ignored  in  limit  equilibrium  methods,  it  must  be  kinematically 
possible  for  the  displacements  to  take  place  in  the  direction  assumed 
along  the  failure  surface  chosen.  Surfaces  composed  of  sections  of 
circles  or  log  spirals  pose  no  kinematic  difficulties.  In  rock  slopes 
the  potential  system  of  failure  surfaces  already  exist  in  the  mass  but 
the  kinematics  of  sliding  must  be  checked  to  delineate  the  possible 
directions  and  surfaces  on  which  it  is  physically  possible  for  sliding 
to  take  place. 
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After  the  potential  failure  surface  is  assumed  in  either  the  rock 
or  soil  slope  stability  analysis,  the  next  step  in  the  limit  equilibrium 


method  is  to  calculate  the  shearing  resistance  required  along  the 
potential  failure  surface  to  keep  the  potential  sliding  mass  in  equilibrium. 
This  portion  of  the  analysis  is  basically  an  exercise  in  statics. 

After  the  shearing  resistance  required  for  equilibrium  has  been 
found,  it  is  compared  with  the  available  shearing  resistance.  This 
comparison  is  usually  expressed  in  terms  of  a  factor  of  safety,  which 
must  be  defined  very  carefully.  Finally  the  slip  surface  giving  the 
lowest  factor  of  safety  is  found.  In  soils  this  is  usually  an  iterative 
process  with  failure  surfaces  of  the  same  shape  but  with  different 
sizes  and  orientations.  But  in  rock  slopes  there  may  only  be  several 
potential  failure  wedges  to  consider,  each  having  a  different  shape 
goverened  by  various  intersections  of  the  sets  of  discontinui ties. 

1.2  Scope 

In  this  report  the  methods  of  analyzing  the  static  stability  of 
rock  slopes  in  three  dimensions  are  given  and  a  method  is  suggested  for 
assessing  the  dynamic  stability  of  rock  slopes.  The  methods  of  static 
analysis  for  three  dimensional  wedges  are  based  primarily  on  the  work 
of  Wittke  (1964,  1965a,  1965b,  1966),  and  Londe  (1965).  Since  vector 
analysis  is  used  in  the  analyses  of  Wittke  and  Londe,  a  review  of  vector 
operations  commonly  used  in  slope  stability  calculations  is  given  in 
Chapter  2.  The  notation  used  for  expressing  strikes,  dips,  etc.  in 
terms  of  vectors  is  also  given  in  Chapter  2. 

In  Chapter  3  various  combinations  of  the  vector  analyses  of  Wittke 
and  Londe  are  presented  for  determining  the  static  factor  of  safety 
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of  rock  slopes.  The  cases  treated  include  slopes  in  rock  masses  containing 
one,  two,  or  three  sets  of  joints.  Example  calculations  are  given  for 
determining  the  factor  of  safety  of  several  typical  problems  by  these 
methods.  The  concept  of  the  dynamic  resistance  of  rock  slopes  is  also 
introduced  in  Chapter  3.  The  method  given  in  Chapter  3  for  computing 
the  dynamic  resistance  of  a  rock  slope  in  three  dimensions  is  original 
with  this  report.  The  dynamic  resistance  can  be  used  for  predicting 
dynamic  displacements  due  to  earthquake  motions  in  the  method  of  dynamic 
analysis  given  by  Newmark  (1965). 

In  Chapter  4  procedures  are  given  for  performing  graphical  solutions 
of  three  dimensional  rock  slope  stability  problems  by  the  use  of 
stereonets.  The  principles  of  the  equal  angle  and  equal  area  projec¬ 
tions  are  reviewed  in  this  chapter  and  the  equal  angle  projection  is 
used  in  this  report  for  the  three  dimensional  analysis  of  rock  slopes. 

The  methods  for  analyzing  rock  wedges  bounded  by  one,  two,  and  three 
joint  planes  are  similar  to  those  given  by  John  (1968)  and  example 
problems  are  illustrated.  In  cases  where  the  static  factor  of  safety 
is  greater  than  unity  a  method  is  also  shown  for  computing  the  magnitude 
and  direction  of  the  limiting  dynamic  resistance  of  a  rock  slope  in 
three  dimensions  by  the  use  of  stereonets.  In  cases  where  the  factor 
of  safety  is  either  less  than  unity  or  less  than  the  desired  value 
a  method  is  also  shown  for  determining  the  optimum  direction  and  magnitude 
of  rock  anchor  or  rock  bolting  forces  required  to  achieve  the  desired 
factor  of  safety. 

In  Chapter  5  procedures  are  given  for  estimating  the  dynamic 
displacement  of  rock  slopes  by  utilizing  the  method  proposed  by 
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Newmark  (1965).  The  minimum  dynamic  resistance  for  rock  slopes  as 
developed  In  Chapter  3  is  used  in  these  calculations.  Guidelines  are 
also  given  for  determining  If  the  dynamic  displacement  calculated  is 
harmful  to  the  stability  of  the  slope. 

In  Chapter  6  a  summary  and  conclusions  are  given. 
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CHAPTER  TWO 


FUNDAMENTALS  OF  VECTOR  ANALYSIS 


2.1  General 

In  this  chapter  the  elements  of  vector  analysis  used  In  three  dimen¬ 
sional  slope  stability  analyses  are  reviewed  to  serve  as  a  ready  reference 
for  the  reader.  Then  the  system  used  in  this  report  for  describing  the 
three-dimensional  orientation  of  joint  planes,  the  line  of  intersection 
of  different  joint  sets,  and  the  resolution  of  forces,  is  introduced  in 
terms  of  vector  notation. 

2.2  Fundamental  Vector  Operations 

A  vector  is  a  quantity  which  possesses  both  a  magnitude  and  a 
direction.  Velocity,  force,  and  momentum  are  examples  of  vector  quan¬ 
tities.  Vectors  of  unit  length  may  also  be  used  to  describe  certain 
reference  directions  such  as  a  normal  to  a  plane  or  the  direction  of 
any  line  with  respect  to  a  set  of  orthogonal  axes.  A  vector  A  may 
be  described  by  the  set  of  its  directional  components  (A  ,  A  >  A  ) 
parallel  to  the  rectangular  Cartesian  axes  (x,y,z).  Thus, 

A  =  (Ax,  Ay)  Az)  (2-J> 

A  vector  may  also  be  expressed  in  terms  of  its  components.  For  example, 

A  -  Ta  +  7a  +  kA  (2-2) 

x  J  y  z 

where  T,  j,  and  k  are  unit  vectors  directed  along  positive  (x,y,z) 
axes  respectively. 

The  magnitude  of  a  vector  A  is  given  by  its  absolute  value  denoted  by 


A 


(A  +  A  + 

'  V  u 


A?) 


1/2 


(2.3) 


'ifei 


Vectors  may  be  added  simply  by  sunning  the  components  in  the  x,  y, 

and  z  directions.  Thus  if  C(C  »C  ,C  )  represents  the  sum  of  two  vectors 

a  y  z 

^(Ax»Ay,Az)  and  §(Bx,By,8z),  then  It  follows  that 

C  -  (lAx  +  jAy  +  kAz)  +  (7Bx  +  JBy  +  kBz) 

-  7(Ax  +  Bx)  +  ](Ay  +  By)  +  t(Az  +  Bz) 

-  TCX  +  ]Cy  +  kCz  (2.4) 

Equating  the  components  in  the  x,  y,  and  z  directions. 

C  ■  A  +  8  ,  C  *  A  +  B  ,  C  -  A  +  B  (2.5) 

x  x  xy  y  y  z  z  z 

The  scalar  product  or  the  dot  product  of  two  vectors  S  and  B  is  denoted 

in  the  form  A*B  and  has  a  magnitude  given  by 

A-B  ■  A  B  +  A  B  +  A  B  (2.6) 

xxyyzz  '  ' 

■  AB  cos  0  (2.7) 

where  9  denotes  the  angle  formed  by  the  vectors  A  and  B  (Fig.  2.1).  The 

scalar  project  is  frequently  used  to  obtain  the  component  of  a  vector  in 

a  given  direction.  For  example,  if  i  is  a  unit  vector  in  the  x  direction, 

A. 7  yields  Ax  ■  A  cos  a1  where  a*  is  the  direction  angle  between  the 

vector  A  and  the  positive  x  axis.  Similarly,  Ay  =  A  cos  B* ,  Az  *  cos  Y* 

where  g1 ,  and  y'  denote  direction  angles  between  the  vector  A  and  the 

positive  y  and  z  axes  respectively.  Substitution  of  these  expressions 

into  Eq.  (2.3)  yields 

cos^a'  +  cos^fj  1  +  cos2y'  “1  (2.8) 

Thus  the  cosines  of  the  direction  angles  (direction  cosines)  of  vector 
A  are  not  independent;  they  must  satisfy  Eq.  (2.8) 

A  vector  product  or  cross  product  of  two  vectors  A  and  §  Is  defined 
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FIG.  2.1  SCALAR  PRODUCT  OF  TWO  VECTORS 
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(2.9) 


l  *  -  ■J- 


to  be  a  third  vector  C  whose  magnitude  is  given  by  the  relation 
C  =  AB  sin  0 

where  9  denotes  the  angle  between  vectors  A  and  §.  The  direction  of  C 
is  perpendicular  to  the  plane  formed  by  vecto-s  A  and  B  as  shown  in  Fig. 
2.2.  The  vector  product  of  A  and  B  is  denoted  in  the  form 

C  =  A  X  B  (2.10) 

where  X  denotes  vector  product  or  cross  product.  The  sense  of  C  is 
such  that  it  is  in  the  direction  a  right  hand  threaded  screw  perpendicular 
to  the  plane  formed  by  A  and  §  would  move  if  A  were  rotated  into  B. 

In  determinant  notation  the  vector  product  is  given  as 


C  =  A  X  B  = 


B 

x 


B 

y 


B 


z 


It  should  be  noted  that 


(2.11) 


AXB  =  -8XA 


(2.12) 


2 . 3  Vector  Operations  Used  in  Three  Dimensional  Analysis  of  Slopes 
2.3.1  Unit  Vectors  Defining  the  Orientation  of  Joint  Planes  and  the 
Line  of  Intersection  of  Joint  Sets 
The  orientation  of  joints  and  planes  of  weakness  are  normally  re¬ 
ported  by  the  field  geologist  in  terms  of  strike  and  dip.  In  this  re¬ 
port,  the  system  given  by  Wittke  (19S4)  will  be  used  to  describe  the 
orientation  of  the  discontinuities  in  relation  to  the  slope  face  being 
investigated.  According  to  this  system,  as  shown  in  Fig.  2.3,  the  x 
axis  is  parallel  to  the  strike  of  the  slope  surface,  the  positive  z 
axis  is  upward  and  the  positive  y  axis  is  directed  toward  the  slope. 
The  strike  of  a  plane  of  weakness  is  given  by  the  angle  8  measured  in 
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a  horizontal  plane  In  a  counterclockwise  direction  from  the  positive 
x-axis  as  shown  In  Fig.  2.3.  The  value  of  8  can  range  between  0  and 
180  degrees.  The  dip  of  a  plane  with  the  horizontal  is  denoted  by  the 
angle  y  In  a  direction  at  90  degrees  to  the  strike.  The  dtp,  y»  can 
range  from  0  to  180  degrees  and  Is  measured  downward  from  a  horizontal 
line  directed  at  an  angle,  \  equal  to  6-90°  to  the  positive  x  axis.  An 
example  of  the  use  of  this  notation  to  describe  the  orientation  of  two 
planes  is  shown  in  Fig.  2.4.  The  strike  and  dip  are  described  by  the 
unit  vectors  u  and  v  respectively,  and  are  written  in  terms  of  the 
angles  6  and  Y  as  shown  in  Fig.  2.3,  i.e., 
u  *>  cos  6  T  +  s 1 n  8  j 


v  >*  cos  Y  s  I  n  8  T  -  cos  Y  cos  8  J  -  s  1  n  Y  k 


or 

u  =  (cos  8 ,  sin  8 ,  0)  (2.13) 

v  -  (cos  Y  sin  8  ,  -  cos  Y  cos  8  ,  -  sin  y)  (2.14) 

Since  the  strike  and  dip  are  at  90  degrees,  the  scalar  product  u’v 
should  be  zero.  u*v  ■  +  cos  8  cos  y  sin  8  -  sin  8  cos  y  cos  8  ■  0. 

Thus  Eqs .  (2.13)  and  (2.14)  satisfy  the  orthogonal  relationship  required 
for  the  unit  vectors  describing  the  strike  and  dip. 

The  cross  product  of  u  and  v  gives  a  unit  vector  w  which  Is  perpen¬ 
dicular  to  both  u  and  v  and  thus  directed  normal  to  the  plane  described 
by  u  and  v.  The  vector  w  Is  obtained  by  expanding  the  determinant 
given  in  Eq.  (2.15) . 

T  j  k 
w  ■  u  X  v  ■  ux  uy  uz 

vx  vy  vz  (2-15) 
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FIG.  2.4  CO-ORDINATE  SYSTEM  FOR  DESCRIBING  STRIKES 
AND  DIPS 
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The  direction  of  w  is  normal  to  the  plane  of  u  and  v  in  the  direction  of 
advance  of  a  right  hand  screw  in  turning  from  u  to  v  through  the  smallest 
angle  between  these  vectors  (<  180°) .  The  magnitude  of  w  is  equal  to  the 
quantity  uv  sin  0  which  assumes  the  vtilue  of  unity  because  6  =  90°  and 
u  and  v  are  unit  vectors.  The  sense  of  w  for  the  two  planes  shown  in 
Fig.  2.4  is  shown  in  Fig.  2.5.  Note  that  for  plane  1  the  direction  of 
Wj  is  normal  to  plane  1  and  directed  downward  into  the  slope  and  w2 
is  normal  to  plane  2  directed  upward  out  of  the  slope.  The  specification 
of  the  unit  vector  w  normal  to  a  plane  is  sufficient  to  completely 
describe  the  orientation  of  that  plane. 

The  direction  of  the  line  of  intersection  of  two  joint  planes 
(planes  1  and  2)  is  given  by  a  vector  Xj2  having  the  direction  of  the 
cross  product  of  the  normal  unit  vectors  to  the  two  planes.  Thus  for 
planes  1  and  2  shown  in  Fig.  2.5  a  vector  along  the  line  of  inter¬ 
section  is  given  by 

x12  “  w2  *  W1  (2.16) 
where  x^2  is  directed  downward  along  the  line  of  intersection  as  shown 
in  Fig.  2.5. 

2.3.2  Resolution  of  Forces 

The  component  of  a  force  R  in  the  direction  given  by  a  unit  vector  n 


is  given  by 


R*n  ■  R  cos  0 


(2.17) 


where  0  is  the  angle  between  R  and  n.  Thus,  for  example,  the  component 
of  a  force  R  normal  to  a  plane  is  given  by  R>w  and  is  given  by 

Rjj  =  R-w  -  Rxwx  +  Rywy  +  Rzwz  (2.18) 
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Similarly  the  component  of  R  In  the  direction  of  the  line  of  Intersection 
-  *12 

of  two  planes  Is  R-~ —  * 

*12 

The  obliquity  of  a  force  R  on  a  given  plane  is  the  angle  which 

the  force  R  makes  with  the  normal  to  the  plane  w  as  shown  in  Fig.  2.6, 
tan  «  Rj/Rfl 

where  R^  and  R^  are  the  ccmponents  of  R  normal  and  tangential  to  the  plane, 
respectively.  Note  that  R^  -  R*w  -  Rx*x  +  R  w  +  R^  *  R  cos  4*  and 


-  |r  X  w|  -  R  sin 

[  (Ryw2"R2wy)2  +  (Rzwx-Rxwz)2  +  (Rxwy-Rywx)2  ]I/2 


Z  X  X  z 

Therefore  the  obliquity  of  a  force  on  a  plane  is  given  by 
rt  [iKwnw^  '  “  '2 

T  I  '  V  7  7  N/‘ 

tan 


(2.19) 


,  S  [(RvVRzwv)  +  (Rxwy-Rywx)  ] 

(Rwx  +  «Twy^R2wz) 


2  1/2 


(2.20) 


^  'xx  y  y 

The  vector  RT  may  also  be  given  by  R-V  which  is  given  by 

RT  -  (rx-Vx)T  +  (Ry"VV)5  *  (Rz-Vz>‘k 

since  Ruw  -  Ruw  7  +  Rww  j  +  Ruw  k.  Thus  the  obliquity  may  also  be  given  as 
N  NX  N  y  N  z 


tan  4 


,  *  <yy/  *  w!  i 1/2 


R  w  +  R  w  +■  R  w„ 
xx  y  y  z  z 


(2.21) 


2.3.3  Line  of  Application  of  a  Force  and  Point  of  Intersection  of  Two 
Forces 

In  order  to  analyze  rotational  stability,  the  point  of  application 
of  a  force  and  its  direction  must  be  known.  If  the  vector  OS  from  the 
selected  origin  of  coordinates  0  to  a  point  5»  on  the  line  of  action  of 
the  force  W  is  known,  the  line  of  action  of  W  may  be  expressed  as  the 
line  joining  the  tips  of  the  set  of  radius  vectors  given  by 
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(2.22) 


r  -  OS  +  XW 
w 

and  shown  in  Fig.  2.7. 

In  a  three-dimensional  problem  the  set  of  applied  forces  will  not 
in  general  intersect,  and  the  moment  of  each  force  about  a  particular 
axis  of  rotation  may  be  considered  separately,  or  the  forces  may  be 
moved  parallel  to  the  axis  of  rotation  about  which  moments  are  being 
surnned  until  the  forces  intersect.  For  example,  in  analyzing  the  rota¬ 
tion  of  a  wedge  as  shown  in  Fig.  2.8  around  the  axis  defined  by  the 
unit  vector  d,  for  the  external  force  P  applied  at  point  N  and  the 
weight  W  applied  through  the  center  of  gravity  S,  either  of  these  forces 
may  be  shifted  any  distance  K  parallel  to  d  without  changing  the  manent 
about  d.  Thus  the  force s  may  be  moved  in  this  manner  until  their  lines 
of  action  intersect.  If  the  line  of  action  of  W  is  defined  by 

r  -  OS  +  XW  (X-  constant)  (2.23) 

w 

and  the  line  of  action  of  P  is  defined  by 

r  -  ON  +  5 P  (5  *  constant)  (2.24) 

P 

the  resultant  R  of  P  and  W  can  be  considered  to  act  at  a  point  of  inter¬ 
section  I  by  setting 

r  -  ?_  +  K-d  (2.25) 

w  p 

Substitution  of  Eqs.  (2.23)  and  (2.24)  in  Eq.  (2.25)  yields 

OS+"5ii-ON+5P+Kd  (2.26) 

If  three  equations  are  written  from  Eq.  (2.26)  in  terms  of  the  x, 
y,  and  z  components  of  P  and  W  they  can  be  solved  simultaneously  for 
X  ■  Xj  >  5  ■  fij  and  K  •  Kj 

which  locate  the  point  of  intersection  I.  The  vector  from  the  origin  0, 
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FIG.  2.8  MOMENT  CAUSED  BY  TWO  FORCES  WITH 
DIFFERENT  POINTS  OF  APPLICATION 


FIG.  2.9  MOMENT  OF  A  FORCE  ABOUT  A  GIVEN  AXIS 


18 


(2.27) 


to  I»  the  point  of  application  of  R,  is  thus 
OT  =  OS  +  XjW 

2.3.4  Moment  about  an  Axis 

The  magnitude  of  the  moment  about  axis  3  through  point  A  as  shown  in 
Fig.  2.9  caused  by  the  force  R  acting  at  point  I  is 

Md  =  (M  X  R)  •  d  (2.28) 

where 


AI  -  AO  +  01 

2.3.5  Point  of  Intersection  of  a  Force  and  a  Joint  Plane 

The  point  of  intersection  of  a  force  and  a  joint  plane  is  found  by 
equating  the  line  of  action  of  the  force  and  the  equation  of  a  plane. 

The  equation  of  a  plane  is  given  by 

rp  •  w  »  constant  (2.29) 

where  rp  Is  a  radius  vector  from  the  origin  to  a  point  In  the  plane, 
and  w  the  unit  vector  normal  to  the  plane.  If  the  vector  OF  from  the 
origin  to  any  point  F  in  the  plane  is  known,  then  the  constant  is  deter¬ 
mined  and  the  equation  of  the  plane  is: 

Tp  .  w  =  (OF-w) 

The  point  where  the  force  P  intersects  the  plane  is  thus  given  by 
solving  simultaneously  the  equation  for  the  line  of  action  of  the  force 
and  the  equation  of  the  plane  giving 

(ON  +  $P)-w  =  (OF-w)  (2.30) 

The  solution  yields  6q>  the  value  of  5  defining  the  piercing  point  Q. 
of  the  force  P  on  the  plane  p  as  shown  in  Fig.  2.10. 

2.3.6  Geometry  of  a  Triangle 

The  area  of  the  triangle  OKL  shown  in  Fig.  2.11  is  given  by 


19 


GEOMETRY  OF  A  TETRAHEDRON 
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(2.31) 


A  =>  1/2  (OK  X  OL| 

and  the  vector  from  0  to  the  centroid,  S,  Is  given  by 

OS  =  1/3  (OK  +  OL)  (2.32) 

2.3.7  Geometry  of  a  Tetrahedron 

The  volume  of  a  tetrahedron  as  shown  In  Fig.  2.12  is  given  by 

V  =  1/6  (OB'  X  DC  |  (hj  +  h2)  (2.33) 

The  centroid  at  point  S  may  be  described  by  the  vector  from  the  origin, 
OS,  given  by 

OS  =  1/4  (OD  +  OC  +  OB)  (2.34) 

The  components  of  OS  are  thus  the  coordinates  of  the  centroid. 
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CHAPTER  THREE 


ANALYSIS  OF  ROCK  SLOPES  BY  VECOTR  METHODS 


3.1  General 

In  this  chapter  analytical  methods  are  presented  for  determining 
the  static  factor  of  safety  of  rock  slopes.  The  cases  covered  include 
rock  slopes  cut  by  one,  two,  or  three  joint  sets.  Example  problems  are 
given  where  various  combinations  of  the  vector  analyses  of  Wittke  and 
Londe  are  utilized.  A  typical  example  problem  is  also  worked  by  common 
engineering  graphics.  The  notion  of  dynamic  resistance  is  also  introduced 
in  this  chapter  for  rock  slopes  and  example  calculations  of  the  minimum 
dynamic  resistance  are  illustrated.  The  methods  given  in  this  chapter 
for  computing  the  dynamic  resistance  of  a  rock  slope  in  three  dimensions 
is  original  with  this  report  and  is  intended  to  be  used  for  predicting 
dynamic  motions  under  earthquake  loadings  in  conjunction  with  Newmark' s 
method  of  analysis  for  the  dynamic  stability  of  slopes. 

3.2  Stability  Calculations  by  Vector  Analysis  for  Sliding  >n  One  Plane 
3.2.1  Calculation  of  Factor  of  Safety  for  Static  Loads 

The  simplest  special  case  of  a  rock  slope  stability  problem  is  where 
the  strike  of  one  of  the  planes  of  weakness  is  parallel  to  the  strike  of 
the  slope  face  as  shown  in  Fig.  3.1.  For  the  coordinate  system  adopted 
in  Chapter  2,  this  condition  can  be  expressed  when  the  unit  vector  u  in  the 
direction  of  the  strike  has  components  of  zero  in  the  y  and  z  directions 
i.e., 

u  -  (ux,  uy,  uz)  =  (1,  0,  0) . 

Then  the  unit  vector  v  in  the  direction  of  the  dip  has  its  x  component 
equal  to  zero,  i.e. , 
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Joint  Plana 


Section  A-A 


FIG.  3.1  SLIDING  ON  ONE  PLANE  -  STRIKE  OF  PLANE 
PARALLEL  TO  STRIKE  OF  SLOPE  FACE 
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The  inclination  of  the  line  of  fall  of  the  plane  ,  e^,  will  determine 

the  kinematic  possibility  of  sliding.  The  angle  of  fall  e  is  given  by 

v  v 

tan  ev  =^T  =  tan(y)  (3.1) 

y 

where  y  is  the  angle  of  dip  of  the  plane.  In  order  for  sliding  to  be 
kinematically  possible,  must  be  smaller  than  a  if  0  <  a  <  n  as  shown 
in  Fig.  3.1(b).  If  a  =  «,  then  must  be  smaller  than  for  the  sliding 
to  be  possible. 

For  a  slope  acted  upon  only  by  gravity  and  the  plane  of  weakness 
striking  parallel  to  the  slope  face,  the  sliding  will  occur  parallel  to 
the  unit  vector  v  in  the  direction  of  the  dip.  The  magnitude  of  the 
component  T  of  the  weight  W  acting  parallel  to  v  may  be  obtained  from 

T  =  W  •  v  (3.2) 

where  W  =  (0,0,  -W) .  The  vector  T  is  given  by 

T  =  Tv  (3.3) 

The  magnitude  of  the  component  of  W  normal  to  the  direction  of  sliding  is 
N  =  W  •  w 

where  w  is  the  unit  vector  normal  to  the  plane  of  sliding  as  given  by 
u  x  v.  The  magnitude  of  the  available  resisting  force  is  given  by  N  tan  <j, 
where  4>  is  the  angle  of  shearing  resistance  between  the  joint  surfaces 
in  the  direction  of  sliding.  The  factor  of  safety  against  sliding  is 
the  quotient  of  the  resisting  and  the  driving  force  in  the  direction  of 
sliding  and  is  given  by 


For  the  case  shown  In  Fig.  3.1,  the  unit  vector  in  the  direction  of 

the  strike  is  given  by  u  =  i  u  =  i  and  the  dip  is  given  by  v  =  jv  + 

x  y 

kvz-  The  unit  vector  w  normal  to  the  plane  of  weakness  is  given  by 

i  j  k  =  -jv  +  kv 

J  J  Z  y 

1  0  0 

0  v  v 

y  z 

Thus  the  magnitude  of  the  component  of  the  weight  in  the  direction  of 
sliding  is  given  by 

T  =  W  •  v  =  -Wvz  (3.5) 


reduces  to 


F.S. 


tan  it 
tan  y 


+  K 


tan  <t> 


tan  <t> 
tan  y 


-  K 


tan  4> 


s  i  n  y 


(3.9) 


where  Y  Is  the  dip  of  the  potential  failure  plane  and  vz  =  -sin  y. 


The  case  may  also  be  considered  where  sliding  can  take  place  on  one 
joint  or  bedding  plane  as  shown  in  Fig.  3.2(a)  or  3.2(b).  In  the  general 
case  the  potential  sliding  wedge  can  be  acted  on  by  its  weight  W,  the 
porewater  force  U  actingon  the  plane  of  sliding,  and  an  external  force  ft 
which  may  be  applied  by  a  structure,  such  as  a  dam.  In  many  cases  where 
we  are  concerned  with  large  slopes,  however,  the  weight  W  will  be  large 
compared  with  ft.  In  the  analysis  of  sliding  on  one  plane  with  the  forces 
W,  U  and  ft  acting  on  the  wedge,  the  forces  are  added  vectorially  into  a 
resultant  R  which  is  given  by 

R  =  W  +  U  +  ft  (3.10) 

The  resisting  reaction  in  plane  a  b  c  as  shown  in  Fig.  3.3  is  R'  and  is 
equal  and  opposite  to  R.  Thus  the  direction  of  sliding  is  in  the  direction 
of  the  projection  of  R  on  plane  a  b  c  and  not  necessarily  in  the  direction 
of  the  dip.  The  angle  of  friction  mobilized,  ^ 1 ,  by  the  force  R  is  given 
by  Eq.  2.20  for  sliding  on  one  plane  as 


tan  $ 


[  (R  w  -R  w  )2  +  (R  w  -R  w  )2  +  (R  w  -R  w  )21  1/2 
,  =  1  y  z  z  y _ z  x  x  z _  xy  y  x  J 


Rw  +  Rw  +  Rw 
xx  y  y  z  z 


Thus  the  factor  of  safeyt  for  this  case  is  given  by 


F.S. 


tan  j> 


(3.11) 


(3.12) 


tan 

3.2.2  Calculation  of  Dynamic  Resistance 

It  should  also  be  noted  that  Wittke  (1965)  has  treated  an. earthquake 
loading  as  an  equivalent  static  load  applied  in  a  horizontal  plane  and 
parallel  to  the  projection  of  the  unit  vector  in  a  direction  of  the  dip. 
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(a)  Cutslope 


(b)  Natural  Slope  —  Bedding  Planes  Dipping  Toward  Valley  —  Sliding  Block 
Isolated  from  Mass  by  Gully  on  Each  Side 


FIG.  3.2  SLIDING  ON  ONE  PLANE -STRIKE  OF  PLANE  NOT 
PARALLEL  TO  STRIKE  OF  SLOPE  FACE 


v,  in  a  horizontal  plane.  Thus,  according  to  the  analysts  of  Wittke  (1965) 
for  earthquake  loading,  the  problem  is  simply  analyzed  as  for  the  general 
case  presented  above  where 

R=W+U+Q+H 

v  v 

where  H  =  i  -  ,  k.W  +  j  - * -  ,k.w 

(v  2  ♦  v  2)T  (v  2  ♦  v  2)i 

x  y  x  y 

The  seismic  coefficient  kj  is  taken  between  0  and  0.2  depending  on 
the  intensity  of  the  earthquake  motion  expected,  and  the  force  H  is  in  a 
horizontal  plane  and  parallel  to  the  projection  of  the  unit  vector  in  a 
direction  of  the  dip,  v,  in  a  horizontal  plane.  The  factor  of  safety  is 
as  given  by  Eq.  3.12.  This  approach,  however,  is  not  recommended  by  the 
authors  since  it  is  considered  as  being  an  unduly  conservative  approach  to 
earthquake  stability.  The  approach  proposed  in  this  report  for  assessing 
the  dynamic  stability  of  rock  slopes  will  essentially  follow  the  concepts 
proposed  by  Newmark  (1965),  which  are  presented  in  Chapter  5.  In  order  to 
ise  the  Newmark  method  of  analysis,  however,  it  is  necessary  to  establish 
the  resistance  available  to  resist  dyanmic  loads.  This  dynamic  resistance 
is  the  resistance  which  is  available  in  addition  to  the  resistance  required 
for  static  stability.  The  dynamic  resistance  is  denoted  by  NW  where  W  is 
is  the  weight  of  the  potential  sliding  block  and  N  is  a  coefficient  to 
be  determined  in  the  following  manner.  The  force  NW  is  that  force  applied 
to  the  potential  sliding  block  which  is  necessary  to  just  make  the  block 
slide  (i.e.  F.S.  =  1).  Depending  on  the  direction  in  which  NW  is  applied, 
its  magnitude  will  vary.  The  magnitude  of  NW  appropriate  for  design  or 
analysis  is  the  magnitude  of  NW  applied  in  such  a  direction  as  to  make  NW 
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a  minimum.  For  a  potential  failure  of  a  block  sliding  on  one  plane  as 
shown  in  Fig.  3.4(a)  NW  should  be  applied  in  a  direction  6  to  the  hori¬ 
zontal  which  will  give  the  minimum  value  of  NW  to  just  cause  the  block 
to  slide.  The  direction  and  magnitude  of  the  minimum  value  of  NW  can 
be  determined  as  shown  in  Fig.  3.4(b).  The  direction  and  magnitude  of 
the  weight  W  is  known  and  the  direction  of  a  resultant  R  is  known  and  is 
inclined  at  $  to  the  normal  of  the  plane  of  sliding  when  sliding  begins 
to  take  place.  Then  the  magnitude  of  the  vector  NW  is  minimum  when  it  joins 
the  tip  of  the  weight  vector  W  in  a  direction  which  makes  an  angle  of  90° 
with  the  resultant  R.  Thus  from  geometry,  the  minimum  magnitude  of  NW  is 
given  by 

NW  =  W  sin  ($  -  y) 


or 

N  =  sin  (<{>  -  y)  (3.13) 

where  <p  is  the  angle  of  shearing  resistance  and  y  is  the  dip.  Thus  the 
minimum  value  of  N  occurs  when  NW  is  in  the  same  direction  of  the  horizontal 
projection  of  the  dip  but  is  inclined  upward  from  the  horizontal  at  an 
angle  of  6  =(<)>-  y)  and  N  has  a  magnitude  of  sin  ($  -  y)  for  the  case  of 
sliding  on  one  plane.  Using  this  minimum  value  of  NW  as  the  dynamic 
resistance  is  a  conservative  estimate  because  it  is  assumed  that  the 
earthquake  motions  are  in  the  most  unfavorable  orientation  for  the  slope 
being  investigated. 

3.3  Example  Problems  of  Sliding  on  One  Plane  by  Vector  Analysis 

Slope  stability  calculations  by  vector  analysis  are  performed  in 
this  section  for  several  cases  Involving  sliding  on  one  plane.  In  some 
of  these  cases  the  same  answer1 could  be  arrived  at  quickly  by  means  of 
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****** 


(b) 


FIG.  3.4  FAILURE  OF  A  BLOCK  SLIDING  ON  ONE  PLANE 


conventional  analysis,  in  others  it  would  be  more  difficult  to  do  by 
conventional  means. 

For  example,  consider  a  wedge  of  rock  sliding  on  a  plane  which  strikes 
East-West  and  dips  30°  South  and  which  has  a  friction  angle  of  40°. 

Consider  the  positive  x  direction  to  be  East,  the  positive  y  direction  to 
be  North  and  the  positive  z  direction  to  be  upwards.  The  unit  vector  in 
the  direction  of  the  strike  is  given  by 
u  =  (1,  0,  0) 

The  unit  vector  in  the  dip  direction  is  given  by 
v  =  (0,  -0.866,  -0.500) 
and  the  unit  vector  normal  to  the  plane  is 

w  =  u  x  v  =  (0,  0.500,  -0.866) 
and  w  is  directed  downward  normal  to  the  plane. 

Case  I 

First  consider  the  factor  of  safety  of  the  block  acted  on  by  its  own 
weight  only.  In  this  case  the  resultant  force  R  acting  on  the  block  is 
equal  to  the  weight  W  given  by 

R  =  W  =  (0,  0,  -W) 

The  magnitude  of  the  component  of  R  normal  to  the  plane  is  given  by 
N  =  R  •  w  =  0.866  W 
then  N  =  0.866  Ww  =  (0,  0.433W,  -0.750W) 
the  tangential  component  of  R  on  the  plane  of  sliding  is 
T  =  R  -  N  =  (0,  -0.433W,  -0.25W) 
the  magnitude  of  T  is  given  by 

T  =  V  [(-0.433)2  +  (-0.25) 2  ] 1/2  =  0.50W 
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then  F.S. 


N  tan  +  0.866W  tan  40° 

T  *  0.50W 


1 .455 


Check:  By  Eg.  3.7 


F.S. 


tan  40° 
tan  30° 


1.455 


Case  II 

Consider  now  that  a  force  A  acts  on  the  wedge  in  addition  to  the  weight 
of  the  wedge  W.  The  force  A  acts  parallel  to  the  strike  (East)  and  has  a 
magnitude  of  0.20W. 

Thus  R  -  A  +  W  «  (0.20W,  0,  -W) 

and  N  «  R  •  w  -  0.866W 

and  N  «■  Nw  *  (0,  0.433W,  -0.75W) 

The  component  of  R  tangential  to  the  plane  is 

T  -  R  -  N  -  (0.200W,  -0.433W,  -0.250M)  and  T  -  0.540W. 

•  r  e  N  tan  4  0.866W  tan  40°  ,  ««. 

•  •  F,S'  ”  T  "  0.54W  ,3 

Note  that  in  this  case  sliding  does  not  take  place  down  the  dip  but  in  the 
direction  of  the  vector  T. 


Case  III 

Consider  now  that  the  wedge  Is  acted  on  by  its  own  weight  and  a  force 
A  having  a  magnitude  of  0.20W  and  acting  in  a  direction  parallel  to  the 
unit  vector  in  the  direction  of  the  dip.  In  this  case  the  normal  component 
is  still  given  by  N  *  W  •  w  ■  0.866W  as  given  in  Case  I.  The  magnitude 
of  the  driving  tangential  force  T  is  the  sum  of  the  magnitudes  of  A  and  the 
tangential  component  of  the  weight  on  the  plane.  The  factor  of  safety 
is  therefore  given  by 


x. 
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1.05 


f  c  0.866W  tin  40° 

“  0.50W  +  0.20W  “ 


Case  Itf 

Now  consider  that  the  plane  under  the  wedge  of  weight  W  as  in  Case 
I  is  acted  upon  by  a  porewater  force,  U,  which  increases  until  the  factor 
of  safety  decreases  from  1.455  to  1.0.  The  porewater  force  does  not  change 
the  driving  force  T.  Therefore,  as  in  Case  I, 

T  -  0. SOW 


The  magnitude  of  the  normal  force  N  as  given  in  Case  I  is  reduced  by  the 
magnitude  of  the  porewater  force,  LI.  That  is 


and  F.S. 


1.0 


N  -  0.866W  -  U 

0.866W  -  U 
0.50W 


tan  40 


and  solving  for  U 


U  *  0.271V 


Check:  By  Eg.  3.9 


F.S. 


1.0 


tan  4> 
tan  Y 


K  tan  4>  tan  40°  tan  40° 

sin  y  tan  30°  ~  K  sin  30u 


Solving  K  =  0.271 

,*.U  -  0.271W 


Case  V 

Consider  the  same  wedge  to  be  acted  on  by  its  own  weight  W,  a  pore¬ 
water  force  U  of  magnitude  0.44W  acting  normal  to  and  on  the  plane  of 
sliding,  and  a  force  A  having  a  magnitude  A  =  0.60W  and  acting  in  a 
direction  S  4S°U  at  a  dip  of  10°.  Then 
W  =  (0,  0,  -W) 

U  =  0 . 44W ( -w)  *>  (0,  -0.22W,  0.371W) 
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The  unit  vector,  a,  in  the  direction  of  force  A  is  given  by  Eq.  2.14  with 
Y  -  170°  and  0  =  135°,  thus, 

a  ■  [-cos  10°  sin  45°,  -(-cos  10°) (-cos  45°),  -  sin  10° ] 
a  =  (-0.696,  -0.696,  -0.174) 


Thus 


and 


A  -  Aa  -  (-0.418W,  -0.418W,  -0.105W) 

R-W  +  U+  A-  (-0.418W,  -0.638W,  -0.734W) 


The  magnitude  of  the  component  of  R  normal  to  the  plane  of  sliding  Is 
given  by 

N  -  R  *  w  -  0.316W 

.*.  N  -  Nw  *=  (0,  0.1S8W,  -0.274W) 

The  component  of  R  tangential  to  the  plane  of  sliding  is 

T  -  R  -  N  -  (-0.41ft#,  -0.796W,  -0.460W) 

T  -  1.009W 

N  tan  40°  0.316W  tan  40° 

“  T  1.009W 

-  0.262 


Case  VI 

Consider  the  slope  acted  on  by  only  its  own  weight  as  in  Case  I.  * 

It  is  desired  to  calculate  the  magnitude  of  the  minimum  dynamic  resistance  NW. 
This  is  simply  given  by  Eq.  3.13  as 

NW  -  W  sin  (4  -  y)  -  W  sin  10° 
or  N  ■  sin  10°  *  0.174 


i 


i 

i 
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Sets  of  Joint  Planes. 


3.4.1  Calculation  of  Factor  of  Safety  for  Static  Loads 

3. 4. 1.1  Description  of  Geometry  and  Loads 

The  general  case  of  two  systems  of  joint  planes  is  as  shown  in  Fig. 

3..S  where  planes  1  and  2  denote  the  joint  planes,  planes  3  and  4  denote 

the  planes  defining  the  faces  of  the  slope,  Yj  and  denote  the  dip  angles 
of  plaries  1  and  2,  Sj  and  6^  denote  the  strike  angles  of  planes  1  and 
2  measured  counterclockwise  from  the  positive  x  direction  and  a  and  6 
denote  the  inclination  of  planes  3  and  4  with  the  horizontal.  The  unit 
vectors  in  the  direction  of  the  strike  planes  1  and  2  are  given  by  Eq.  2.13 

Uj  ■  (cos  8j,  sin  6j,  0)  u2  =  (cos  82>  sin  62,  0) 

and  the  unit  vectors  in  the  direction  of  the  dip  for  planes  1  and  2  are 
given  by  Eq.  2.14; 

Vj  =  (cos  Yj  sin  8j,  -  cos  Yj  cos  Bj,  -  sin  Yj) 

v2  -  (cos  Y2  sin  B2,  -  cos  Y2  cos  82,  -  sin  y2) 

The  unit  vectors  normal  to  each  plane  are  given  by 

Wj  -  Uj  X  Vj 

and  w2  “  u2  x  v2 

Note  on  Fig.  3.5  that  Wj  Is  directed  downward  into  plane  1  and  w2  is 

directed  outward  from  plane  2  when  the  normals  are  defined  in  this  manner. 

Also  note  that  the  plane  designated  as  plane  1  is  the  one  with  the  lowest 
value  of  8,  In  the  case  where  the  strikes  of  two  planes  are  the  same  the 

plane  designated  as  plane  1  is  the  one  with  the  smallest  value  of  y-  This 

convention  is  necessary  to  maintain  the  proper  sign  convention  for  the 
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The  loading  of  the  slope  consists  of  (1)  dead  load  W  acting  at  the 
center  of  gravity  of  the  mass  considered,  (2)  live  load  Q  applied  at  any 
point,  (3)  porewater  forces  Uj  and  U2  acting  on  planes  1  and  2  respectively 
and  (4)  dynamic  loads  induced  by  ground  motions  from  earthquakes  or  nuclear 
detonations.  The  resultant  R  of  the  loads  in  any  given  case  can  be 
determined,  and  let  the  point  of  appl icat ion  of  the  resultant  be  at  point  I. 
3.4. 1.2  Determination  of  the  Mode  of  Sliding  Failure 

For  the  case  of  a  tetrahedron  bounded  by  two  base  planes  which  may 
be  'intersect ing  joint  sets,  failure  may  occur  by  sliding  along  the  line  of 
intersection  of  the  two  planes  or  by  sliding  on  either  one  of  the  two 
planes . 

The  first  step  in  determining  the  mode  of  failure  is  to  check  if  the 
distrubing  forces  tend  to  lift  the  tetrahedron  from  either  or  both  of 
the  supporting  planes.  Thus  considering  the  rock  wedge  OBCD  (Fig.  3.5), 
the  resultant  force  R  tends  to  break  the  contact  between  the  tetrahedron 
and  planes  1  and  2  respectively  if 
R  •  w,  <  0 

1  (3.14) 

and  R  •  W2  >  0 

If  Eqs.  3.14  show  that  the  resultant  force  R  tends  to  lift  the  tetrahedron 
off  of  both  supports,  then  equilibrium  is  not  possible  unless  the  joints 
can  take  tension  or  rock  bolts  are  added  to  take  the  computed  tension. 
Normally  this  will  not  happen  for  large  slopes  acted  on  by  their  own 
weight  and  porepressures,  but  could  occur  for  small  tetrahedrons  near 
the  surface  of  steep  or  overhanging  slopes.  If  Eqs.  3.14  show  that  lifting 
occurs  off  of  one  of  the  supporting  planes  then  we  can  definitely  say 
that  sliding  cannot  occur  on  that  plane. 
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pJSt 


If  Eqs.  3.14  show  that  lifting  off  of  the  wedge  from  the  supporting 
planes  does  not  occur,  i.e. 


R  •  Wj  >  0 

R  •  w2  <  0 


(3.14a) 


then  we  must  make  further  kinematic  tests  to  see  whether  sliding  takes 
place  on  plane  1  only  or  plane  2  only  or  along  the  line  of  intersection 
of  planes  1  and  2. 

In  order  to  evaluate  the  mode  of  sliding  it  is  necessary  to  define 
two  new  vectors  jSj2  *nd  2^12  w*'*c*'  are  given  by 


1 S 12  “  x12  *  W1 


2S 12  ”  x12  x  w2 


(3.15) 


and  are  as  shown  in  Fig.  3.5(b).  The  vector  jSj2  is  in  plane  1  perpendicular 
to  the  line  of  intersection  xJ2  and  the  vector  2§12  is  in  plane  2  per¬ 
pendicular  to  the  line  of  intersection  Xj2- 

If  sliding  is  to  occur  along  the  line  of  intersection  Xj2>  then 
Eqs.  3.16,  3.17  and  3.18  must  be  satisfied  simultaneously. 


R  •  jS12  > 

0 

(3.16) 

R  ■  2S12  > 

0 

(3.17) 

€  <  a  if 

X 

0  <  a<  n  and  <  S  if  a  ■  it 

(3.18) 

-1 

e  -  tan 

(~ ) 

(3.19) 

X 

x12y 

and  x12y’  x12z  "  y  an<*  z  comPonents  of  vector  Xj2 

The  vector  Xj2  along  the  line  of  intersection  is  defined  in  Chapter  2  and 

is  given  by  Eq.  2.16  as 

x12  “  w2  x  W1 
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If  sliding  is  to  occur  on  plane  1  only,  then  both  the  following  equations 


must  be  satisfied: 

R  '  Wj  >  0  (3.20) 

and  R  •  <  0  (3.21) 

Similarly  if  sliding  is  to  occur  on  plane  2  only,  then  Eqs.  3.22  and  3.23 
must  be  satisfied. 

R  •  w2  <  0  (3.22) 

R  •  2S12  <  0  (3.23) 

The  physical  interpretation  of  Eqs.  3.16  -  3.23  may  be  made  as  follows. 
Eq.  3.16  is  satisfied  only  if  the  resultant  force  R  has  a  component  which 
tends  to  push  the  wedge  on  plane  1  toward  the  line  of  intersection  Xj2- 
Similarly  Eq.  3.17  is  satisfied  only  if  the  resultant  force  R  has  a  com¬ 
ponent  pushing  the  wedge  on  plane  2  toward  the  line  of  intersection  Xj2. 

Thus  Eqs.  3.16  and  3.17  ensure  that  the  resultant  force  R  wedges  the 
tetrahedron  between  the  two  plane  so  that  sliding  can  only  take  place  on 
both  the  planes  along  the  line  of  intersection.  In  order  for  sliding 
along  the  line  of  intersection  to  be  kinematically  possible,  it  should  also 
be  ensured  that  the  line  of  intersection  does  not  plunge  into  the  rock 
slope  and  this  check  is  provided  by  Eq.  3.18.  Thus  when  all  the  three 
kinematic  conditions  specified  by  Eqs.  3.16  through  3.18  are  satisfied 
simultaneously  sliding  can  occur  along  the  line  of  intersection.  The 
tendency  to  slide  will  be  downhill  if  R  *  Xj2  >  0  and  uphill  if  R  •  Xj2  <  0 
(Fig.  3.5). 

Eq.  3.21  indicates  a  component  of  R  on  plane  1  tending  to  move  the 
block  away  from  plane  2  by  sliding  on  plane  1  and  Eq.  3.20  establishes 
the  condition  for  contact  on  plane  1.  Thus  Eqs.  3.20  and  3.21  are  suf- 
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flcient  and  necessary  conditions  for  sliding  to  occur  on  plane  1.  Similarly 
Eqs.  3.22  and  3.23  specify  the  conditions  for  sliding  on  plane  2. 

3. 4. 1.3  Calculation  of  the  Factor  of  Safety  for  Sliding 

If  the  kinematic  tests  discussed  above  show  that  sliding  takes 
place  on  only  plane  1  or  on  only  plane  2,  then  the  factor  of  safety  can 
be  computed  from  Eq.  3.4  for  sliding  on  one  plane.  Thus  for  sliding  on 
plane  1  the  factor  of  safety  may  be  computed  as 


F.S. 


Nj  tan 


(R  •  Wj)  tan 


(3.24) 


where 


f, .  s  -  (R  •  -  tIxt  ♦  r  j  *  Tli; 


Thus  Eq.  3.24  becomes 


F.S. 


which  may  be  written  as: 


(R  ‘  *j)  tan 


2  2  T\ 

Tlx  +  Tly  +  TIz 


1/2 


(3.25) 


F.S. 


tan  4,[R  w,  +  R  w,  +  R  w,  ] 
1  x  lx  y  ly  z  lzJ 


2  2  21/2 
((Ryw,2-Rzwiy)  +  «V*lx-Rxwlz)  +  (Rxwly-Rywlx)  ] 


(3.26) 


For  sliding  on  plane  2  only,  the  factor  of  safety  is  given  as 


F.S.  - 


Nz  tan  4 j  -(R  '  *2)  tan 


(3.27) 


'2  2 

The  minus  sign  appearing  in  Eq.  3.27  is  due  to  the  direction  of  the  unit 
normal  w2  as  shown  in  Fig.  3.5. 


Eq.  3.27  can  be  expanded  to  yield 
F.S. 


tan  »2[ (-Bxw2x-ay»2y-Rzw2z)] 


[(RyW2z‘RzW2y)  +  (RZw2x_Rxw2z)  +  ] 


2,  1/2 


(3.28) 
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If  the  kinematic  tests  of  Eqs.  3.16,  3.17  and  3.18  are  satisfied 
and  sliding  takes  place  on  both  planes  1  and  2  along  the  line  of  inter¬ 
section  xj2»  then  the  factor  of  safety  may  be  computed  in  the  following 
manner. 

The  first  step  is  to  compute  the  magnitude  of  the  driving  force,  Tj2, 
shown  in  Fig.  3.6,  in  the  direction  of  sliding.  This  is  simply  given  by 

R  •  x, 


'12 


'12 


'12 


(3.29) 


where  X|2  represents  the  magnitude  of  the  vector  x^. 
is  in  the  same  direction  as  x12  and  is  given  by 

-  T12x12 


The  vector  Tj2 

(3.30) 


It  is  convenient  to  define  the  vector  NJ2,  normal  to  the  line  of  inter¬ 
section  which  is  given  by 

Nj2  *  R  -  Tj2  (3.31) 

In  order  to  evaluate  the  frictional  resistances  on  planes  1  and  2,  it  is 
necessary  to  determine  the  components  Nj  and  N2  of  NJ2  acting  normal 
to  planes  1  and  2  respectively.  The  relationship  of  the  vectors  R,  Tj2> 
^12’  *^1  an{*  N2  are  s^10wn  ’n  Sections  AA  and  BB  of  Fig.  3.6.  From  Fig. 

3.6  it  is  obvious  that 

NjWi  +  N2(-w2)  =  fi12  (3.32) 

where  Nj  and  N2  represent  the  magnitudes  of  the  two  component  vectors  Nj 
and  N2  respectively. 

Thus 


Nlwlx  '  N2W2x  =  N12x 
Nlwly  -  N2W2y  *  N12y 

Nlwlz  '  N2u2z  ‘  N12z 


(3.33) 

(3.34) 

(3.35) 
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Any  two  of  Eqs.  3.33,  3.34,  and  3.35  can  be  used  to  determine  Nj  and  N2 
and  the  third  equation  can  be  used  to  check  the  numerical  values  of 
and  Nj.  After  Nj  and  N2  are  obtained  the  factor  of  safety  for  sliding 
on  both  planes  may  be  determined  from  the  equation 

N,  tan  +  N,  tan  4>, 

F.S.  -  — - i-r: - - - -  (3.36) 

12 

3.4. 1.4  Calculation  of  Static  Factor  of  Safety  for  Rotations 

In  addition  to  the  previously  investigated  sliding  movements  which 
endanger  stability,  the  rock  wedge  OBCD  may  rotate  about  the  support 
edges,  OC  or  00,  or  about  the  axes  at  point  0  perpendicular  to  planes  1 
and  2,  when  the  resultant  load  exerts  an  overturning  moment  about  these 
axes  (Fig.  3.7).  Even  though  all  the  above  modes  of  failure  by  rotation 
are  conceivable,  under  normal  conditions  the  most  probably  axes  of  rotation 
are  djQ  and  d2Q  (Fig.  3.7)  and  therefore  consideration  is  given  only  to 
rotations  about  these  two  axes  in  this  section.  The  treatment  of  rotation 
about  OC,  OD,  djQ  or  d2g  is  similiar  and  is  not  developed  in  this  report. 
The  axes  of  rotation  djg  and  d2Q  pass  through  0  and  are  perpendicular 
to  planes  1  and  2  respectively.  In  a  rotation,  say  about  the  djg  axis, 
all  points  of  the  wedge  in  the  region  of  the  area  ODB  move  tangential  to 
plane  1  while  the  surface  OCB  of  the  rock  wedge  separates  from  plane  2. 

The  equations  of  the  djg  and  d2g  axes  are  obtained  as  follows: 

d1Q  -  -Wj  (3.37) 

d20  *  -w2  (3.38) 

In  the  analysis  for  rotations,  it  is  necessary  to  know  the  points 
of  application  of  the  various  forces  acting  on  the  rock  wedge  OBCD  so 
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1*  1 i-.-t  .  .1 


FIG. 3.7  ROTATIONAL  STABILITY  OF  A  WEDGE  BOUNDED  BY 
TWO  JOINT  PLANES 
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that  the  point  of  application  I  (Fig.  2.8)  of  the  resultant  force  R 
can  be  determined.  The  weight  W  acts  vertically  downwards  at  the  center 
of  gravity  S  of  the  rock  wedge  as  shown  in  Fig.  2.8.  The  vector  Us  as 
shown  in  Fig.  2.8  can  be  obtained  from  geometrical  considerations  as 

OS  =  1/4(0D  +  OC  +  OB)  (3.39) 

_  _  _  (also  Eq.  2.34) 

where  the  vectors  ODi  OC  and  OB  are  given  by  the  following  equations: 


=  ^tan  a  tan  Bj  tan  Yj  sin  Bj  ’  tan  a  *  hl^ 

__  hj  hj 

°C  "  ^tan  a  tan  02  tan  sin  B  2*  tan  a’  ^  (3-41) 

(3.42) 


OB  -  (h.  +  hj 

x12z  1 


tan  a  -  tan  e 


x  tan  6 


2  tan  -  tans  tan  a  1 


(3.43) 


where  hj,  h2>  a,  6,  Yj»  Bp  Y2>  8nd  &2  8re  def!ned  in  Fig.  3.9. 


The  weight  of  the  rock  wedge  can  be  determined  from  its  volume  V  which  is 
given  by: 

V  »  7  | OB'  x  DC  |  (h,  +  h2)  (3.44) 

(also  Eq.  2.33) 

where 

DC  =  OC  -  OD  (3.45) 

OB'  =  OB'  -  OD  (3.46) 

OB'  -  hj  (3.47) 

x12z 

The  point  of  application,  I,  of  the  resultant  force  R  is  determined 
from  the  known  magnitudes  and  lines  of  action  of  the  component  forces 
by  using  Eq.  (2.27)  and  the  principles  of  vector  analysis  as  explained 
in  Chapter  2. 
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For  a  rotation  to  be  possible  about  the  djQ  axis,  the  resultant 
force  R  must  have  a  positive  scalar  component  of  moment  about  the  xJ2  and 
djQ  axes  as  evaluated  by  Eq.  (2.28);  i.e., 

■  moment  of  R  about  Xj2  =  Xj2'(OIxR)>  0  (3.48) 

and  ^dlO  *  m<>ment  of  R  about  djQ  =  djQ-(oixR)>  0  (3.49) 

Similarly  the  moments  of  R  about  the  Xj2  and  d^Q  axes  have  to  satisfy 

Eqs.  (3.50)  and  (3.51)  if  a  rotation  is  to  occur  about  the  d2Q  axis. 

Mx  =  *12  ‘  (°*x^  4  0  (3.50) 

and  Md20  *  ^20  *  >  0  (3.51) 

In  addition  a  few  kinematic  tests  must  also  be  satisfied  and  these 
tests  are  dependent  on  the  magnitude  of  the  angles  n ,  kjg  and  k2g  which 
are  defined  as  follows: 

n  -  wedge  angle  between  planes  1  and  2 

=  cos-1  (wrw2)  0<  n  <  n  (3.52) 

ki0  "  DOB  =  cos"1  (^||  0<  k1Q<  n  (3.53) 

k20“  COB  =  cos'1  0<  k20<  n  (3.54) 

The  range  of  angles  <1,  kjQ  and  k2Q  for  which  a  rotation  about  the  dJ0  and 
d2Q  axes  is  kinematically  impossible,  is  given  in  Table  3.1. 

The  analyses  for  determining  the  static  factor  of  safety  for  rotations 
about  dj0  and  d2^  axes  are  similar  in  principle  and  therefore  the  details 
of  the  analysis  will  be  given  only  for  the  case  of  rotation  about  the  d^Q 
axis. 

The  resultant  R  is  first  resolved  into  components  Nj  and  Tj  at 
its  point  of  intersection.  Q,  with  plane  1,  as  shown  in  Fig.  3.7(a).  Thus 
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Axis  of  Rotation:  d 


10 


n 

k10 

o 

r>4 

Supplementary 

Cond i t i on 

0<rj<jr 

>»r/2 

>jt/2 

- 

0<n<* 

<«/ 2 

>n/2 

- 

<«/2 

>jt/2 

<n/2 

tan  k2Q 

Axis  of  Rotation:  d2Q 

n 

k10  k20 

Supplementary 

Condition 

0<n<n 

>x/2  >«/2 

0<n<it 

>x/2  <r/2 

Mn  k10  ,  % 

- - 7 - ; - r->  seclx-n) 

tan(jt-k20) 

<«/2 

<«/2  >n/2 

Table  3.1  Range  of  Angles  for  which  a  rotation 
is  kinematically  Impossible. 
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Nj  =  (R-Wj)^  (3.55) 

and 

Tj  =  R  -  Nj  (3.56) 

The  component  Tj  tangential  to  plane  1  is  now  resolved  into  components  Tr 
and  Fig.  3.7 (b) | .  The  force  T^  has  the  direction  of  the  vector  Q.0  and 
the  force  T  has  the  direction  of  the  tangent  to  the  rotation  which  Q 
executes  in  the  case  of  a  rotation  about  djQ.  The  force  T  is  thus 
the  only  component  of  the  loading  which  exerts  an  overturning  moment 
about  the  djQ  axis.  The  resolution  of  force  Tj  into  its  components  T 
is  done  as  follows: 

Tj  -  Tf  +  Tt  -  Cj(-OQ)  +  c2(OC  x  Wj)  (3.57) 

In  Eq.  3.57,  -OQ  and  OQ  x  Wj  are  vectors  in  the  direction  of  Tr  and  Tt« 

By  equating  the  x,  y  and  z  components  of  Tj  as  given  by  Eqs.  3.56  and 
3.57,  the  values  of  the  two  coefficients  Cj  and  c2  may  be  determined. 

Eqs.  3.56  and  3.57  give  three  equations  for  the  two  unknowns  Cj  and  c2 
and  therefore  one  of  these  equations  can  be  used  to  check  the  calculations 
for  Cj  and  c2>  With  Cj  and  c2  known,  Tr  and  Tt  are  obtained  as  follows: 

Tr  w  -Cj  OQ  (3.58) 

Tt  «  c2(Od  x  wj)  (3.59) 

The  magnitude  of  the  overturning  moment  M^jq  can  be  obtained  by  the  relation 

MdlO  "  Tt  0,1  (3-60) 

The  magnitude  of  the  restoring  moment  of  the  frictional  force  on  plane  1 
due  to  the  normal  component  Nj  is  obtained  as 

MrdlO  "  N1  tan  *1  0(1  (3-61) 
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The  factor  of  safety  against  rotation  can  now  be  obtained  as  the  ratio  of 
the  restoring  moment  to  the  overturning  moment 

N ,  tan  <(>,  OQ. 

F.S. {against  rotation)  = — - — ^ -  —  Njtan<J>1/Tt  (3.62) 

The  factor  of  safety  against  rotation  about  the  d^Q  axis  can  also  be 
determined  in  a  similar  manner.  Moments  M^q  and  Md2Q  are  very  often 
negative  and  in  these  cases  only  the  sliding  stability  need  be  analyzed. 

3.4.2  Calculation  of  Dynamic  Resistance  Against  Sliding  on  Two  Planes 
The  direction  and  magnitude  of  the  minimum  dynamic  resistance  NW 
which  is  necessary  to  just  make  the  potential  block  slide  on  the  two 
base  planes  may  be  found  by  the  following  procedure. 

A  unit  vector  rj  is  first  defined  in  the  direction  of  the  resultant 
reaction  Rj  on  plane  1  (Fig.  3.8).  In  the  limiting  state  of  equilibrium 
Rj  is  inclined  at  an  angle  4>j  to  the  upward  normal  -Wj  to  plane  1  and 
tends  to  oppose  the  downward  movement  along  the  line  of  intersection  xJ2. 
Therefore , 

rl  =  "w!  cos  ^1  “  X12  s'n  ^l^x12  (3.63) 

Similarly  a  unit  vector  r2  is  defined  in  the  direction  of  the  resultant 
reaction  R2  on  plane  2.  From  Fig.  3.8  it  can  be  seen  that 

r2  =  w2  C0S  ^2  *  x12  s,rl  <*>2^x12  (3-64) 

where  $2  is  the  angle  of  friction  on  plane  2. 

The  magnitude  of  the  dynamic  resistance  vector  NW  will  be  a  minimum 
when  the  vector  NW  is  normal  to  the  plane  containing  Rj  and  Rj.  Therefore 
a  unit  vector  n  in  the  direction  of  NW  may  be  obtained  by  the  equation: 

n  =  (r{  x  r2) / 1  r-j  x  r2  |  (3.65) 
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r,  * -w,  •  COS  p!rj  sin*,  rz  ■  Wj-cos**-^  sin*2 


FIG.  3.8  STABILITY  OF  A  WEDGE  BOUNDED  BY  TWO 
JOINT  PLANES 
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where  jrj  x  r2 |  represents  the  magnitude  of  the  vector  (rj  x  ?2) . 

The  magnitude  of  the  minimum  dynamic  resistance,  NW,  may  now  be  deter 
mined  by  the  equation 

NW  =  R  •  n  (3.66) 

where  R  is  the  resultant  of  all  static  loads  acting  on  the  sliding  rock 
wedge.  From  Eq.  3.66  it  follows  that 

N  =  (3.67) 

w 

3.4.3  Example  Problems  for  Slopes  with  Two  Intersecting  Planes  of  Dis¬ 
continuity  Worked  by  Vector  Analysis. 

Problem  1 

Determine  the  factor  of  safety  of  the  rock  wedge  OBCD  shown  in  Fig. 
3.9.  Also  estimate  the  direction  and  magnitude  of  the  minimum  dynamic 
resistance  NW  which  is  necessary  to  just  make  the  potential  block  OBCD 
slide. 


Plane  1 

Plane  2 

0 1  =  36° 

02  =  94° 

V,  -  62° 

y  2  =  121' 

4>  J  =  20° 

*2  =  40° 

a  =  70° 

6  =  20° 

Solution 


Static  Factor  of  safety  against  sliding 

According  to  Eqs.  2.13,  2.14  and  2.15,  for  plane  1, 


FIG.  3.9  STABILITY  OF  A  ROCK  WEDGE  BOUNDED  BY  TWO 
JOINT  PLANES 
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(2.13) 


Gj  =  (0.809,  0.588,  0.000) 


"V1  “ 

:  (0.276,  -0.380,  - 

■0.883) 

(2. 

.14) 

"1  ' 

.  (-0.519,  0.714,  - 

0.469) 

(2. 

.15) 

and 

for  plane 

2, 

“2  = 

■  (-0.070,  0.998, 

0.000) 

(2. 

.13) 

v2  = 

■  (-0.514,  -0.036, 

-0.857) 

(2. 

14) 

w2  = 

(-0.855,  -0.060, 

0.515) 

(2. 

15) 

The 

only 

load 

that  enters  the  calculation 

in  this  problem  is 

the  weight 

W  of  the 

rock 

wedge  OBCD 

i  acting  vertically 

1  downwards  in  the  - 

■z  di recti 

ion. 

Therefore 

the 

resultant 

load  R  may  be  expressed  as 

R  = 

(0,  0,  -W) 

R  • 

Wj  =  0.469W  >  0 

and 

R  • 

w2  =  -0.515W  <  0 

(3. 

14a) 

Therefore 

lift 

:ing  off  of 

the  rock  wedge  from  the  support  planes  does 

not 

occur 

• 

*12 

=  W2  X  Wj 

=  (-0.340,  -0.669, 

-0.642) 

(2. 

16) 

and 

*12 

=  0.987 

1 S 1 2 

*  *12  *  W1 

(3. 

15) 

.  (0.722,  0.172,  - 

0.590) 

and 

2  S 1 2 

=  Xj2  X  W2 

(3. 

15) 

=  (-0.383,  0.725, 

-0.551) 

R-jS 

12  =  0.590W  >  0 

(3. 

16) 

r-2s 

12  -  0.551W  >  0 

(3. 

17) 
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-1  /  12zv 


-  I  / 

€  “  tan  ^ J 
x  X12y 


“1  /-0.642,  oC 

tan  (^OS9}  '  43 -8 


(3.19) 


<$  <  <  <  a 

x 


(3.18) 

Thus  according  to  Eqs.  3.16,  3.17,  and  3.18,  sliding  is  kinematically 
possible  only  along  the  line  of  intersection  Xj2*  Since  R'Xj2  “  0.642W  > 

sliding  tends  to  occur  down  the  line  of  intersection. 


T12  =  R  •  ^12/*i2  “  0.642W/0.987 


(3.29) 


0.650W 


and 


T. _  -  0.650W  —  -  (-0.223W,  -0.440W,  -0.420W)  (3.30) 

U  x12 


’12 


R  -  T12  -  (0.223W,  0.440W,  -0.580W) 


Nlwl  +  N2^"*2^ 


(3.31) 

(3.32) 


*  ^(-0.519  ,  0.714,  -0.469)  + 


N2 (0.855,  0.060,  -0.515) 


Sol v i nq 


Nj  -  0.565W,  N2  *  0.605W 


F.S. 


Nj  tan  ♦j  +  N2  tan  $2 


(3.36) 


1  12 


0.565W  tan  20°  +  0.605W  tan  40° 
0.650W 


1.10 


Stability  against  rotation. 

According  to  Eqs.  3.39  through  3.43 


0D  =  (-0.404hj ,  0.346hj,  hj) 


(3.40) 


[jjlj 


(3.41) 


OC  «=  (0.602hj,  0.364hj,  hj) 

OC  =  1.220^ 

OB  -  (0. 741hj ,  1.460h1(  1.398hj)  (3.42) 

06  =  2. 155hj 

OS  *  (0. 235hj ,  0.547hj ,  0.850hj)  (3.39) 


In  order  to  apply  the  kinematic  tests  for  rotation,  it  is  necessary  to 
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20 


the  angles  kjg, 

k2Q  and  n 

• 

-1 

=  cos 

0D-0B 

(ocRo® 

=  48.1°  < 

n/2 

(3.53) 

-1 

=  cos 

OC-OB 

(oc5To5 

*  25.3°  < 

n/2 

(3.54) 

-1 

=  cos 

(wj  •  w2) 

=  80.9°  < 

n/2 

(3.52) 

=  x12 

(OIxR) 

(3.48) 

=  XI2 

(OS  x  R) 

=  Q.03Whj 

>  0 

;10’  k20 

and  M  , 

X 

a  rotation 

about  axis  dj0 

i  s 

kinematically  possible.  However  the  rotation  can  occur  only  if  M^jq  >  0. 

djQ  =  -Wj  =  (0.519,  -0.714,  0.469)  (3.37) 

01  =  OS  =  (0. 235hj ,  0.5471^,  0.850hj) 


Mdl0  =  B10  '  (01  x 


(3.49) 


*  -0.452Whj  <  0 

Therefore  rotation  about  djg  axis  does  not  occur. 


Minimum  Dynamic  Resistance 

The  unit  vector  rj  in  the  direction  of  the  resultant  reaction  on 
plane  1  is  given  by  Eq.  3.63  as 

^  =  (0.607,  -0.439,  0.663)  (3.63) 
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The  unit  vector  r2  In  the  direction  of  the  resultant  reaction  on  plane 
2  can  be  obtained  in  a  similar  manner  from  Eq.  3.64 

r2  -  (-0.434,  0.388,  0.810)  (3.64) 

The  unit  vector  in  the  direction  of  the  minimum  dynamic  resistance  vector 
NW  is  then  given  by 

?i  x  r2 

n  ,  -i. - -  »,  (-0.616,  -0.785  ,  0.046)  (3.65) 


Solut ion 
Case  (a)  P  ■  0 

Static  Factor  of  Safety  against  sliding 
For  plane  1 , 

Gj  -  (0.955,  0.292,  0.000)  (2.13) 
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Vj  =  (0.146,  -0.478,  -0.866) 

(2.14) 

Wj  -  (-0.253,  0.827,  -0.499) 

(2.15) 

For  plane  2, 

u2  -  (0.454,  0.890,  0.000) 

(2.13) 

;2  *  (0.155,  -0.079,  -0.985) 

(2.14) 

w2  -  (-0.877,  0.447,  -0.174) 

(2.15) 

When  P  =  0 

R  *  (0,  0,  -W) 

where  W  =  weight  of 

the  rock  wedge  OBCD 

*12  *=  «2  x  "l 

(2.16) 

=  (-0.079,  -0.394,  -0.594) 

x12  «  0.717 

1^12  "  x12  *  W1 

(3.15) 

=  (0.688,  0.110,  -0.165) 

2§12  =  *12  *  *2 

(3.15) 

=  (0.334,  0.507,  -0.381) 

R-jS12  =  0.165W  >  0  and 

(3.16) 

R.2312  =  0.381W  >  0 

(3.17) 

-1  ,*12z,  .  -1  /-0.594\ 

€x  -  tan  (— )  =  tan  (.0<394) 

(3.19) 

=  56,4° 

.  .  6  <  <  <  a 

(3.18) 
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Thus  according  to  Eqs.  3.16,  3.17  and  3.18,  sliding  is  kinematically 


possible  only  along  the  line  of  intersection 

Since  R’X^  “  0-594W  >  o,  sliding  tends  to  occur  down  the  line  of 
sect  ion. 

T12  **  R  •  *l2/xl2  15  0.S49W/0.717 

=  0.828W 

x12 

T„  =  0.828W  ~ 

1 1  x12 

«=  (-0.091W,  -0.45SW,  -0.696W) 

NJ2  =  R  -  Tj2 

=  (0.091W,  0.455W,  -0.314W) 

=  NjWj  +  N2(-w2) 

=  Nj (-0.253 ,  0.827,  -0.499)  + 

N2(0.877,  -0.477,  0.174) 

Solving  Nj  -  0.733W,  N2  =  0.314W 

0.733  tan  30°  +  0.314W  tan  30° 

F-5-  “  0.S28W 

=  0.73  <  1 

Stability  against  Rotation. 

00  =  (-23.70,  0,  12.00) 

00  =  26.60 

OC  *  (-2.40,  0,  12.00) 

OC  =  12.25 

OB  =  _^594  (-0.079,  -0.394,  -0.594) 

-  (1.60,  8.00,  12.00) 
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inter- 

(3.29) 

(3.30) 

(3.31) 


(3.40) 

(3.41) 

(3.42) 


(3.39) 


06  =  14.50 

OS  =  |  (OB  +  OC  +  OD) 

-  (-6.10,  2.00,  9.00) 

kjO  -  cos  (qD-OB^  =  74  <  «/2 

k20  =  C0S  (oToB)  *  38  n/2 

n  =  cos  l(wj  •  w2)  =  47.4°  <  «/2 
=  Xj2  *  (OS  x  R) 

=  2.561W  >  0  indicating  thereby  that  the 
R  intersects  plane  1.  For  these  values  of  n,  k^,  k2Q,  and  a 
about  axis  djQ  is  kinemat ical ly  possible. 

djo  =  -Wj  =  (0.253,  -0.827,  0.499) 

Mdl0  =  510  •  G*  *  R) 

=  4.539W  >  0 

Therefore  a  rotation  can  occur  about  the  djQ  axis  and  the  factor 
against  rotation  can  be  determined  as  follows: 

Nj  =  R  •  Wj  =  0.499W 

Nj  =  (R  •  Wj)wj  =  (-0.126W,  0.413W,  -0.250W) 

Tj  =  R  -  Nj  =  (0.126W,  -0.413W,  -0.750W) 

«*  Cj(-OQ)  +  c2(0Q  x  Wj) 

OQ  *=  01  +  fR  -  OS  +  YR 
*  [-6.1,  2.0,  (9.0-hW)] 

6) 


(3.53) 

(3.54) 
(3.52) 

(3.48) 

resul  tant 
rotation 

(3.37) 

(3.49) 

of  safety 

(3.55) 

(3.56) 

(3.57) 


Since 


OQ  and  Wj  are  mutually  perpendicular,  0(J  •  w  ■  0 

.  .  (6.1  x  0.253)  +  (2.0  x  0.827)  -0,499  (9.0  -  1  w)  «  0 

.  .  (9.0  -  W)  =  6.40 

JQ  =  (-6.10,  2.00,  6.40) 

OQ  x  Wj  =  (-6.29,  -4.67,  -4.54) 


Sovl ing 


Tj  =  (0.126W,  -0.413W,  -0.750W) 

=  Cj (6.10,  -2.00,  -6.40)  +  c2(-6.29,  .4.67,  -4.54) 
c j  =  0.078W;  c2  =  0.055W 

Tt  =  c2(0Q  x  Wj)  (3.59) 

=  (-0.346W,  -0.2S7W,  -0.250W) 


Tt  -  0.498W 


F.S. 


Nj  tan 


0.499W  x  tan  30° 
0.498W 


0.58 


(3.62) 


Note:  It  may  be  noted  that  all  the  lengths  in  the  above  case  are  expressed 
in  feet-units. 

Case  (b)  P  =  (0,  10,  0) 

In  Case  (a)  the  only  force  in  the  system  is  the  weight  W  of  the  rock 
wedge  and  it  is  not  necessary  to  know  the  magnitude  of  W  for  estimating 
the  factor  of  safety  of  the  rock  wedge.  But  in  Case  (b)  there  is  an 
additional  external  force  P  of  magnitude  10  tons  acting  in  the  positive 
y-direction  through  the  center  of  gravity,  S,  of  the  rock  wedge  and  there¬ 
fore  it  becomes  necessary,  in  the  present  case,  to  compute  the  magnitude 
of  W. 
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DC  *  OC  -  00 

=  (21.30,  0,  0) 

_  *12 

OB-  =  -f-  h 

x12r 

=  (1.60,  7.96,  12.00) 

DB-  =0B'  -  00 

«  (25.30,  7.96,  0) 

V  -  ||DB-  x  DC| (h|  +  h2) 

»  339.1  ft3 

W  a - 2000 -  27,13  tOOS 

W  *  (0,  0,  -27.13) 

P  =  (0,  10,  o) 

R  =  W  +  P  =  (0,  10,  -27.13) 

R  «  jSu  =  5.58  >0 

R  •  =  15.41  >  0 


6  <  e  <  a 

x 


(3.43) 

(3.45) 


(3.47) 


(3.46) 

(3.46) 


(3.44) 


(3.16) 

(3.17) 


7Tie  above  values  show  that  sliding  is  kinematically  possibly  only  along 
the  line  of  intersection  of  planes  1  and  2.  Since  R  *  *12  *  12>48  >  °» 
sliding  tends  to  occur  down  the  line  of  intersection. 


(3.29) 


T12  =  R  •  *12/*i2  =  12.48/0.717  =  17.41 
Tj2  =  17.41  ;j2/x12  =  (-1.918,  -9.567,  -14.423)  (3.30) 

N i 2  =  R  -  Tu  =  (1.918,  19.567,  -12.707)  (3.31) 

=  N^-0.253,  0.827,  -0.499)  + 

N2(0.877,  -0.447,  0.174) 

Solving  Nj  =  29.4  tons  N2  =  10.7  tons 

_  29.4  tan  30°  +  10.7  tan  30° 

'  17.41 

=  1.33 

Stability  against  Rotation. 

Mx  =  X12  ‘  (°s  x 
=  112.84  >  0 

.  .  The  resultant  R  intersects  plane  1  as  in  Case  (a), 
n  =  47.4°  <  n/2 

k10  =  74°  <  n/2 

k2Q  =  38°  <  n/2 

For  these  values  of  n,  kjQ,  k^,  and  M^,  a  rotation  about  djg  is  kine¬ 
matically  possible. 

d10  =  "*1  =  <°-253>  -°-827>  0.499)  (3.37) 

Mdl0  =  d10  ’  ^0S  *  ^ 

=  70.0  >  0 

Therefore  a  rotation  tends  to  occur  about  the  5jq  axis.  The  factor  of 
safety  aginst  rotation  can  be  determined  as  follows. 

Nj  =  R  .  Wj  =  21,83 
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e  *£vm-  »■  • 


Nj  =  Nj  Wj  =  (-5.523,  18.053,  -10.893) 

Tj  =  R  -  Nj  =  (5.523,  -8.053,  -16.237) 

*  Cj(-OQ)  +  c2(00  x  Wj) 

OQ  -  01  +HR  =.  OS  +  W 

=  [-6.1,  (2.0  +  10 'l) ,  (9.0  -  27. m)] 


(3.55) 

(3.56) 

(3.57) 


00  and  Wj  are  mutually  perpendicular 

.  .  OQ  .  Wj  *  0 

.  .  (6.1  x  0.253)  +  (2.0  +  10  ^  0.827  -  0.499(9.0  -  27.13'f)  -  0 

Solving  H' =  0.0593 

.  .  00  =  (-6. 10,  2.59,  7.39) 

00  x  Wj  =  (-7.41,  -4.91,  -4.39) 

Tj  =  (5.523,  -8.053,-16.237) 

*  Cj (6. 10,  -2.59,  -7.39)  + 

c2(-7.41»  -4.91,  -4.39) 

Solving  Cj  =  1.77  c2  «  0.71 

Tt  *  c2(W  X  ;,)  (3.59) 

=  0,71 (-7.41 ,  -4.91,  -4.39) 

=  (-5.261,  -3.496,  -3.117) 

T  =  7.04  tons 


F.S.  = 


N1  tan  *1  21.83  tan  30° 


Thus  the  provision  of  the  lateral  force  P  increases  the  stability  of  the 
wedge  OBCD  against  both  sliding  and  rotation. 

Note:  In  case  (b)  all  the  forces  are  in  ton-units  and  all  the  lengths  are 
in  feet-units. 
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3 . 5  Analysis  for  Sliding  on  Two  Planes  by  Engineering  Graphics 

The  factor  of  safety  of  a  rock  wedge  sliding  on  the  two  base  planes 
can  also  be  determined  graphically  by  using  the  principles  of  engineering 
descriptive  geometry.  To  illustrate  the  procedure.  Problem  1  of  section 
3.4.3  will  be  solved  using  this  method.  The  details  of  this  graphical 
•elut'en  *r»  siiown  in  Fig.  3.11. 

Th*  lines  DB'  and  CB1  represent  the  strikes  of  planes  1  and  2  inclined 
at  angles  Bj  and  with  the  front  of  the  slope,  I-I.  An  edge  view  of 
each  plane  is  drawn  as  an  auxiliary  elevation  to  locate  the  position  of 
a  point  0  common  to  both  the  planes  situated  at  any  depth,  d,  below  the 
horizontal  plane  OCB1 .  Since  B1  is  also  a  point  common  to  both  the  planes, 
B'O  represents  the  line  of  intersection  of  planes  1  and  2.  A  side  elevation 
parallel  to  8'd  gives  the  true  dip  of  the  line  of  intersection  B'O.  The 
weight  vector  W  is  then  resolved  into  components  Nj2  and  Tj2,  respectively 
normal  and  parallel  to  the  line  of  intersection  B'O  as  shown  in  Fig.  3.11. 

An  auxiliary  elevation  of  the  two  planes  looking  in  the  direction  OB'  is 
obtained  and  the  components  Nj  and  N2  of  NJ2  normal  to  planes  1  and  2 
respectively  are  then  determined.  Once  the  magnitudes  of  Nj  and  N2  are 
known,  the  factor  of  safety  is  computed  using  the  relationship 

N,  tan  $ .  +  N,  tan  4, 

F.S  .=— - 7 - - - -  (3.36) 

'12 

3 .6  Method  of  Stability  Analysis  for  Rock  Slopes  with  Three  Intersecting 

Joint  Sets, 

In  this  section,  the  stability  against  sliding  of  a  tetrahedral  volume 
of  rock,  ABCO  bounded  by  three  planes  1,  2  and  3  and  an  exterior  surface 
ABO,  is  investigated  by  using  Londe's  method  of  analysis  (F  g,  3.12). 
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W  =  66 

V  ^2 

N,  =  38  fan  20°  =  .364 

N8=  40  tan  40°  =.840 


F.  S. 


47,4 

42.0 


1.12 


38  x, 364  =  13.8 

40  x  .840  =  33.6 
47.4 


FIG.  3.11  GRAPHICAL  SOLUTION  OF  SLIDING  STABILITY  OF 
A  ROCK  WEDGE  BOUNDED  BY  TWO  JOINT  PLANES 
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FIG.  3.12  FORCES  ON  A  ROCK  WEOGE  BOUNDED  BY  THREE 
INTERSECTING  JOINT  PLANES 


Sliding  failure  of  the  tetrahedral  rock  mass  ABCD,  can  occur  by  separation 
from  one  or  two  of  the  three  bounding  planes.  There  are  thus  six  possible 
modes  of  sliding  failure  as  shown  in  Fig.  3.13.  The  mode  of  failure  in  a 
given  case  will  depend  on  the  geometry  of  the  problem  and  the  magnitude 
and  direction  of  the  resultant  of  the  applied  forces,  R,  as  defined  by 


the  equation: 


w  +  i  +  Oj  +  u2  +  u3 


(3.68) 


Where  W  ■  (W  ,  W  ,  W  )  =  total  weight  vector  of  the  tetrahedral  volume  of 
x  y  7  3 


Q  ■  (Qx»  Qy,  Qz)  ■  any  externally  applied  force  on  the  rock  wedge 


ul’  U2’  U3 


hydrostatic  uplift  or  porewater  forces  that  act  on 
planes  1,  2  and  3  resepect i vely 


The  first  step  in  the  stability  analysis  of  the  rock  wedge  ABCD,  is  to  deter¬ 
mine  the  mode  of  sliding  failure  for  a  given  set  of  input  conditions.  This 
can  be  done  as  explained  in  the  following  section. 

3.6.1  Determination  of  the  Mode  of  Sliding  Failure  (Fig.  3.14) 

Let  Wj ,  w2  and  represent  unit  vectors  normal  to  planes  1,  2  and  3 
respectively,  directed  towards  the  inside  of  the  rock  volume.  The  resultant 
force  R  lifts  the  tetrahedron  from  all  three  contact  faces  if  all  the  three 


following  equations  are  satisfied  simultaneously 


R  •  wj  >  0 


R  •  w2  >  0 


R  •  w3  >  ° 


(3.69) 

(3.70) 


(3.71) 


In  such  a  case  equilibrium  is  not  possible  unless  the  joints  can  take  tension 
or  rock  bolts  are  provided  to  resist  the  tensile  forces  across  the  faces. 


>**!■.  V 


.T  • 


NATURE  OF 
SLIDING 


in  Plana  I 
Direction  Between 
CA  and  CB 


CONTACT 

FACES 


OPEN 

FACES 


DIAGRAM 


In  Plone  2 
Direction  Between 
CD  and  CB 


FIG.  3.13  MODES  OF  SLIDING  FAILURE  OF  A  ROCK  WEDGE  BOUNDED 
BY  THREE  INTERSECTING  JOINT  SETS 
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If  Eqs.  3.69  through  3.71  show  that  lifting  off  of  the  wedge  from  the  support 
planes  does  not  occur,  then  further  kinematic  tests  must  be  made  to  determine 
the  mode  of  sliding  failure. 

The  vectors  Xj2>  *23  an(*  *31  a*on9  the  lines  of  intersection  CO,  CB 
and  CA  are  given  by  the  following  equations: 

(3.72) 

(3.73) 

(3.74) 

Let  us  now  define  two  new  vectors,  jSjj  and  2^12  orthogonal  to  Xj2  and  lying 
in  planes  1  and  2  respectively  as  follows: 

(3.75) 


x j 2  -  w2  x  Wj 


*23  “  w3  x  w2 


x3 1  “  W1  x  w3 


1512  ”  x12  X  "l 


2S12  “  w2  x  x12 


(3.761 


Similarly  the  vectors  2^3  and  3S23  normal  to  x23  and  lying  in  planes  2 
and  3  respectively  are  given  by 


2S23  *  x23  x  W2 


3S23  “  w3  *  x23 


(3.77) 

(3.78) 


The  vectors  and  3§3j  normal  to  x^j  and  lying  in  planes  1  and  3  re¬ 

spectively  are  similarly  given  by 

(3.79) 

(3.80) 

The  orientations  of  all  the  vectors  defined  by  Eqs.  3.72  through  3.80  are 
shown  in  Fig.  3 . 14. 


3S3 1  “  x3 1  x  w3 


1 S3 1  *  W1  x  x3 1 
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the  resultant 


If  sliding  is  to  occur  along  the  line  of  intersection  Xj2» 

R  of  the  applied  forces  must  have  a  component  along  Xj2  which  tends  to  open 
up  plane  3.  This  condition  can  be  represented  vectorially  by  the  equation 


R  •  Xj2  _>  0  (3.81) 

In  addition  the  components  of  R  on  plane  1  and  2  along  vectors  and 

2§j2  must  be  directed  towards  the  line  of  intersection  Xj2-  In  other  words, 


R  •  tS12>  0 

(3.82) 

R  *  2S12  1  0 

(3.83) 

Eqs.  3.81  through  3.83  must  be  satisfied  simultaneously  if  the  rock  wedge 

ABCD  is  to  slide  along  Xj2  with  face  3  open. 

The  conditions  to  be  satisfied 

for  sliding  to  occur  along  x23  and  x^  can  be 

They  are: 

obtained  in  a  similar  manner. 

For  sliding  along  x23.  R  *  x23 

(3.84) 

*  •  2*23^  0 

(3.85) 

S  •  3'S23 

(3.86) 

For  sliding  along  x3j,  R  *  x^j  0 

(3.87) 

*  •  iSii  0 

(3.88) 

5  •  3^31—  0 

(3.89) 

If  sliding  is  to  occur  on  any  one  plane  only, 

say  on  plane  1,  then  R  must 

have  a  component  normal  to  plane  1  directed  towards  the  outside  of  the  rock 

wedge  ABCD.  In  other  words 


In  addition  the  components  of  R  on  plane  1  along  and  must  be 
directed  away  from  x ^  and  x^j.  In  other  words 

(3.91) 

(3.92) 


1S12  i° 


1^31  ~  ° 


The  corresponding  equations  for  cases  of  sliding  on  planes  2  and  3  are  as 
follows. 

For  sliding  on  plane  2  only: 


R  '  w2  1  0 

(3.93) 

R  •  2^12  —  0 

(3.94) 

*  '  2§23  i  0 

(3. 95) 

plane  3  only: 

R  •  w3  i  0 

(3.96) 

*  •  3*23  i  0 

(3.97) 

5  •  3531  -  0 

(3.98) 

3.6.2  Calculation  of  the  Factor  of  Safety  for  Sliding 

After  deciding  on  the  mode  of  sliding  failure  based  on  the  kinematic 
tests  mentioned  above,  the  next  step  is  to  estimate  the  factor  of  safety 
against  sliding  under  the  given  conditions.  The  procedure  for  estimating 
the  factor  of  safety  is  basically  the  same  as  that  explained  in  Section  3. 4. 1.3 
of  this  chapter  for  the  case  of  a  rock  wedge  bounded  by  two  joint  planes. 

Three  example  problems  have  been  added  to  illustrate  the  method  of  analysis. 

In  the  preceding  analysis,  of  Section  3.6.1,  however,  it  has  been  tactily 
assumed  that  the  critical  •'ock  wedge  is  bounded  by  all  the  three  joint  planes 
and  the  exterior  slope  face  as  shown  in  Fig.  3.12.  When  the  field  conditions 
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are  such  that  this  assunption  is  valid,  the  method  of  stability  analysis 
presented  above  is  directly  applicable.  But  in  a  majority  of  cases,  it 
is  likely  that  the  critical  rock  wedge  is  bounded  by  two  (rather  than  by  all 
the  three)  Joint  planes.  Under  these  conditions  the  stability  analysis  has 
to  be  performed  as  explained  in  Section  3. 4. 1.3. 

3.6.3  Example  Problems  for  Slopes  with  Three  Intersecting  Planes  of 
Discontinuity  Worked  by  Vector  Analysis 
Problem  1 

Determine  the  factor  of  safety  against  sliding  of  the  rock  wedge  ABCD 
shown  in  Fig.  3.12.  Also  estimate  the  direction  and  magnitude  of  the  minimum 
dynamic  resistance  MW  which  Is  necessary  to  just  make  the  potential  block 
ABCD  slide. 

wj  -  (0.00,  0.72,  0.69) 

w2  -  (0.63,  -0.12,  0.77) 

w3  -  (0.00,  0.00  1.00) 

W  -  (0,  0,  -36.5  tons)  0  -  (0,  0,  0) 

Uj  23.6  tons  U2  ■  8.0  tons  Uj  »  5.7  tons 

-  40°  4 2  -  40°  *3  -  40° 


Solut ion 


r-W  +  Q  +  Uj'+Uj  +  Uj  (3.68) 

-  W  +  Q  +  UjWj  +  U2w2  +  U3w3 
■  (5,05,  16.04,  -8.34)  all  in  ton-units 


R  -  18.8  tons 

;12  -  (-0.638,  -0.435,  0.454) 


'12 


0.895 


(3.72) 
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x23  -  (0.120,  0.630,  0) 

(3.73) 

*23  ”  0-6^0 

;31  -  (0.720,  0,  0) 

(3.74) 

*31  “  °-720 

jS12  -  (-0.626,  0.440,  -0.459) 

(3.75) 

2S12  -  (0.280,  -0.777,  -0.351) 

(3.76) 

2§2g  -  (0.485,  -0.093,  -0.410) 

(3.77) 

3§2g  =  (-0.630,  0.120,  0) 

(3.78) 

gS31  -  (0,  -0.720,  0) 

(3.79) 

S  =  (0,  0.497,  -0.518) 

(3.80) 

1  31 

R  '  Wg  *=  -8.34  <  0 

(3.96) 

Thus  plane  3  Is  closed  and  a  failure  by  lifting  from  the  base  planes  does 
not  occur. 


R  •  gS23  =  -1.25  <  0 

(3.97) 

R  '  3S31  =  -11.55  <  0 

(3.98) 

Eqs.  3.96  through  3.98  thus  indicate  that  sliding  can  occur  only  on  plane  3. 
Ng  =  R  •  (-Wg)  »  8.34  tons 

Ng  -  Ng(-Wg)  -  (o,  0,  -8.34) 

T3  -  R  -  Ng  -  (5.05,  16.04  ,  0) 

Tg  -  16.8  tons 

Ng  tan  $3  g .34  tan  40° 

F.S. - j- -  -  - TO 

-  0.42 
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Problem  2 


Work  out  Problem  1  with  the  following  changes: 

Uj  -  12.0  tons  U2  ■  2.0  tons  llg 


Solution 

R  -  W  +  Q.  +  Uj  +  U2  +  Ug 

“  W  +  Q  +  UjWj  +  U2W2  +  UgWg 

-  (1.26,  8.40,  -24.68) 

R  •  w3  «=  -24.68  <  0 


2.0  tons 


(3.68) 


(3.71) 


and  therefore  lifting  off  of  the  rock  wedge  from  all  the  base  planes  is  not 
possible. 


R  •  x23  "  5.45  >  0 

(3.84) 

R  •  2S23  *  9.94  >  0 

(3.85) 

R  •  2S23  -  0.22  >  0 

(3.86) 

The  above  equations  show  that  sliding  can  occur  only  along  the  line  of 
intersection 

T23  »  R  •  ^23^x23  “  ®,5°  tons 


I  ,  ■  T„  •  ^  =  (1.60  ,  8. 35  ,  0) 

“  x23 

N23  -  R  -  T23  -  (-0.34,  -.05,  -24.68) 
■  H2(-w2)  +  N3(-w3) 

-  N2(-0.63,  0.12,  -0.77)  + 

n3  (0,  0,  -1.00) 
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24.26  tons 


Solving  Nj  “  0.54  tons  ■ 


Problem  3 

A  rock  cut  slope  runs  East-West  and  the  three  major  joint  sets  inter¬ 
secting  the  slope  have  the  following  orientations: 

Joint  Plane  Strike  D i P 

1  N47°E  44°SE 

2  N20°W  83°SW 

3  16°SW 


The  angle  of  shearing  resistance  on  all  the  three  joint  planes  is  estimated 
to  be  20°.  Determine  the  factor  of  safety  of  the  slope  against  a  sliding 
f  a  i  1  u  re . 


Solution 

Consider  the  positive  x  direction  to  be  East,  the  positive  y  direction 
to  be  North  and  the  positive  z  direction  to  be  upwards.  Then  the  three  joint 
planes  have  the  following  strike  and  dip  angles. 

Plane  1  6  j  =  47°  Yj  “  44° 

Plane  2  =  Yj  *  ®7° 

Plane  3  ®3  “  *59°  Y3  * 

The  unit  normals  to  planes  1,  2  and  3  can  be  defined  by  Eqs.  2.13,  2.14 
and  2.15.  When  these  normals  are  oriented  such  that  they  are  directed 
toward  the  interior  of  the  rock  wedge  they  are  defined  by  the  following 
equations: 
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Wj  =  (0.474,  -0,508,  0.719) 
w2  =  (-0.933,  -0.339,  0.122) 
w3  =  (-0.099,  -0.257,  0.961) 


x12  =  (-0. 182,  0.729,  0.635) 

(3.72) 

=  0.983 

x23  =  (0.295,  -0.885,  -0.206) 

(3.73) 

x23  =  0.955 

i3l  =  (-0.303,  -0.526,  -0.172) 

(3.74) 

x3l  =  0631 

jS12  =  (0.847,  0.432,  -0.253) 

(3.75) 

2S12  =  (-0.304,  0.570,  -0.741) 

(3.76) 

2S23  =  (-0.178,  -0.157,  -0.925) 

(3.77) 

3S23  =  (0.903,  0.263,  0.163) 

(3.78) 

3§31  =  (-0.550,  0.309,  0.026) 

(3.79) 

jS31  =  (0.466,  -0.137,  -0.403) 

(3.80) 

R  =  (0,  0,  -W) 

R  •  Wj  =  -0.719W  <  0 


R  •  w2  =  -0.122W  <  0 


R  •  w3  =  -0.961W  <  0 

Therefore  failure  by  lifting  off  of  all  the  base  planes  is  not  possible  as 
shown  by  comparison  of  the  above  three  equations  with  Eqs.  3.69,  3./0, 
3.71.  It  can  easily  be  verified  that  all  kinematic  tests  are  satisfied 
only  for  sliding  on  plane  3.  In  other  words 


(3.96) 


R  •  w3  =  -0.961W  <  0 

R  *  3S23  =  -0.I63W  <  0  (3.97) 

R  ‘  3S3 1  =  -0.026W  <  0  (3.98) 

N3  =  R  *  (-w3)  =  0.961  W 

N  =  N0(-w,)  =  (0.095W,  0.247W,  -0.924W) 

3  J  J 

T3  =  R  -  N3  =  (-0.095W,  -0.247W,  -0.076W) 

T3  =  0.275W 

0.961W  tan  20°  _ 
l"-b-  “  0.275W  * 

The  preceding  calculations  have  been  carried  out  under  the  assumption 
that  the  critical  rock  wedgeis  bounded  by  all  the  three  joint  planes.  As 
has  been  pointed  out  earlier,  in  a  majority  of  cases,  there  exists  a  rock 
wedge,  bounded  by  only  two  joint  planes,  which  is  more  critical  than  the  one 
considered  in  the  preceding  analysis.  As  a  matter  of  fact,  in  the  present 
problem,  the  rock  wedge  bounded  by  planes  1  and  2  has  a  lower  factor  of 
safety  with  respect  to  sliding.  The  determination  of  the  mode  of  failure 
and  the  factor  of  safety  against  sliding  can  be  done  as  explained  in  section 
3.4. 1.3.  The  details  of  this  analysis  will  not  be  given  here  except  the 
fact  that  the  sliding  tends  to  occur  down  the  line  of  intersection  of  planes 
1  and  2  and  that  the  factor  of  safety  is  0.58  as  compared  to  the  previous 
value  of  1.27. 

3 . 7  Computer  Techniques 

The  stability  of  rock  slopes  bounded  by  two  or  three  joint  sets  can 
also  be  analyzed  using  digital  computer  techniques.  This  method  avoids 
lengthy  hand-calculations  and  is  particularly  useful  when  there  is  a  need 
for  solving  a  whole  series  of  stability  problems. 
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The  basis  of  the  procedure  is  the  same  as  explained  in  the  previous 
sub-sections.  The  essential  steps  In  this  procedure  are  as  follows  (see 
flow  chart,  Fig.  3.15): 

1.  Using  the  input  data  calculate  all  the  required  directional  vector 
quant i ties. 

2.  Check  to  see  if  failure  by  lifting  off  the  base  planes  of  the  rock  wedge 
is  possible. 

3.  If  not,  determine  the  probably  mode  of  sliding  failure. 

4.  Calculate  the  factor  of  safety  for  this  mode  of  sliding  failure. 

5.  Check  for  stability  against  the  possible  mode  of  rotation. 

6.  Print  the  results. 

A  documentation  and  listing  of  the  computer  programs  using  Fortran  IV 
language  is  given  in  Appendix  A. 


Hi 


Calc.  Req'd 

Directional 

Victor* 


Determine  Mode 
of  Sliding 
Failure 


Cole.  FS 


Print  Resulte 


k 
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Chapter  Four 

GRAPHICAL  SLOPE  STABILITY  ANALYSIS  BY  USE  OF  STEREONETS 

4. 1  Properties  of  Spherical  Projections 

4.1.1  General 

The  orientation  (strike  and  dip)  of  planes  or  lines  in  space 
can  be  represented  by  the  intersection  of  the  plane  or  line  with  the 
surface  of  a  reference  sphere  through  whose  center  the  plane  or  line 
passes.  As  can  be  seen  in  Figure  4.1,  the  intersection  of  a  plane  with 
the  sphere  is  a  great  circle,  while  a  line  which  parallels  the  plane 
will  plot  as  two  points,  180  degrees  apart,  on  the  great  circle.  A 
plane  can  also  be  represented  by  the  intersection  of  its  normal  with 
the  sphere  (the  pole  of  the  plane),  which  will  plot  as  a  point  located 
90  degrees  from  the  great  circle,  in  both  the  upper  and  lower  hemispheres 
of  the  sphere. 

To  communicate  this  information,  a  two-dimensional  representation 
of  the  spherical  projection  is  necessary.  Several  types  of  projection 
can  be  used  to  transfer  great  circles  and  points  from  the  spherical 
surface  to  the  equatorial  plane  of  the  sphere. 

The  equal  angle  projection  (termed  a  Wulff  net  or  stereonet)  is 
the  method  used  in  this  report  because  of  the  simplicity  in  plotting 
the  projections.  Each  great  circle  on  the  sphere  plots  as  an  arc  of  a 
circle  on  the  equatorial  plane  of  the  sphere. 

Another  type  of  projection,  the  equal  area  projection,  is  used  for 
compiling  statistical  information  on  the  frequency  and  orientation  of 
lines  or  planes.  It  therefore  should  be  used  to  plot  and  evaluate  the 
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FIG.  4. 


Line  Parallel  To  Plane  A 
Pole  Of  Plane  A 

Equatorial  Plane 


Projection  Of  Plane  A 
On  Lower  Hemisphere 


I  PROJECTION  OF  PLANE  AND  LINES  ON 
A  SPHERE 
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raw  data  from  field  and  borehole  mapping  of  joints  and  other  geologic 
discontinuities.  The  equal  area  projection  of  great  circles  from  the 
sphere  to  the  equatorial  plane  results  in  a  distortion  from  the  circular 
arc>  and  therefore  is  not  quite  as  simple  to  use  for  stability  analyses 
as  the  equal  angle  projection. 

4.1.2  Equal  Angle  Projections 

Figures  4.2  and  4.3  show  the  lower  hemisphere,  egdal  angle  method 
for  projecting  a  point  from  the  surface  of  the  sphere  to  equatorial 
plane. 

A  line  is  drawn  from  point  P  on  the  sphere  to  the  upper  pole,  U, 
of  the  equatorial  plane  (dashed  line  in  Fig.  4.2  and  4.3).  The  inter¬ 
section  of  this  line  with  the  equatorial  plane  (P1 )  is  the  desired  pro¬ 
jection  of  point  P.  In  Figure  4.2,  the  projection  of  plane  A  and  point 
P  from  the  lower  hemisphere  to  the  equatorial  plane  is  shown;  the  projec¬ 
tion  of  plane  A  plots  as  an  arc  of  a  circle  (or  line  of  meridian)  on  the 
equatorial  plane. 

The  projection  of  a  vertical  plane  will  project  as  a  straight  line 
through  the  origin  of  the  equatorial  plane.  A  horizontal  plane  will  project 
as  a  line  of  meridian  having  a  radius  equal  to  the  radius  of  the  sphere, 
with  the  same  origin.  All  points  projected  from  the  lower  hemisphere 
will  plot  within  this  circle  on  the  equatorial  plane.  Points  from  the 
upper  hemisphere  projected  on  the  equatorial  plane  will  plot  outside  the 
radius  of  the  sphere,  as  can  be  seen  for  the  projection,  Q.*  ,  of  point  Q 
in  Figure  4.3. 

A  diagram  of  the  stereonet  obtained  from  an  equal  angle,  lower 
hemisphere  projection  is  shown  in  Fig.  4.4.  The  lines  of  meridian 
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O'  it  a  Projection  of  Point  0  ( 0  is  Above  *he  Equatorial 
Plane  )  P'  it  o  Projection  of  Point  P  (  P  is  Below  the 
Equotoriot  Plane  ) 


Upper  Pole 


FIG. 4.3  PROFILE  OF  SPHERE  SHOWING  METHOD  OF 
EQUAL  ANGLE,  LOWER  HEMISPHERE 
PROJECTION 
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FIG.  4.4  STEREONET  (  WULFF  NET)  (EQUAL  ANGLE  PROJECTION) 
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through  the  N  and  S  poles  of  this  diagram  represent  great  circles  resulting 
from  the  Intersection  of  N-S-  striking  planes  with  the  reference  spnere. 

Dip  angles  for  these  planes  are  shown  on  the  E-W  axis  of  the  stereonet. 

The  meridians  for  steeply  dipping  planes  will  approach  straight  lines 
on  this  plot,  while  the  meridians  for  flat-lying  planes  will  plot  as 
arcs  of  circles  having  radii  approaching  the  radius  of  the  reference 
sphere.  Each  of  the  meridians  is  divided  into  180  degrees  by  E-W  lines 
of  latitude,  which  plot  as  arcs  of  circles  on  the  equal  angle  stereonet. 

To  represent  a  plane  which  strikes  other  than  N-S.  the  stereonet  of 
Fig.  4.4  must  be  rotated  so  that  Its  N-S  axis  is  aligned  In  the  direction 
of  the  strike  of  the  given  plane.  The  meridian  can  then  be  traced  from 
the  stereonet  so  that  it  is  oriented  in  its  proper  strike  direction. 

Note  that  the  true  dip  of  a  plane  or  line  should  be  determined  by  orient¬ 
ing  the  E-W  axis  of  the  stereonet  so  that  it  is  in  the  direction  of  the 
dip  of  the  line  or  plane. 

Stereonets  similar  to  that  shown  in  Fig.  4.4  are  available  from 
graphic  aid  suppliers.  It  is  suggested  that  such  a  stereonet  be  used  for 
the  example  problems  of  this  report  by  overlaying  clear  vellum  on  the 
stereonet  and  rotating  the  stereonet  about  its  center,  beneath  the  sheet 
of  vellum,  to  plot  planes  and  lines  of  various  strikes  and  dips. 

In  Figure  4.4  the  great  circle  projection  of  Plane  A  (dipping  40° 
west  and  striking  N-S)  plots  as  a  line  of  meridian.  The  pole  (or  normal) 
of  plane  A  is  located  90°  from  the  plane.  Lines  parallel  to  plane  A 
plot  as  points  on  this  line  of  meridian.  The  angle  between  two  such 
lines,  OP  and  OR  is  100°  and  is  found  by  counting  the  lines  of  latitude 
along  the  meridian,  between  points  R'  and  P' . 
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Figure  4.5  shows  the  projection  of  two  planes  on  the  stereonet,  one 
striking  N-S,  the  other  N  42°  E.  The  orientation  of  the  line  of  inter¬ 
section  of  the  two  planes  is  determined  from  the  point  of  intersection 
of  the  two  meridians.  In  this  case  the  line  of  intersection  dips  at  an 
angle  24°  in  a  direction  of  S  32°  W.  All  of  this  information  can  be 
determined  by  using  the  stereonet,  rotating  it  as  required  to  plot  lines 
of  meridian  and  read  angles.  The  dip  angle  is  read  by  rotating  the 
stereonet  until  either  the  NS  or  EW  axis  coincides  with  the  direction  of 
dip.  The  dip  angle  is  then  read  in  degrees  from  the  outer  edge  of  the 
stereonet. 

4.2  Use  of  Stereonet  to  Evaluate  Driving  and  Resisting  Forces  on  a 

Potential  Sliding  Wedge  of  Rock 

The  use  of  the  stereonet  in  stability  analyses  has  been  described 
by  John  (1968),  Goodman  (1964).  The  stereonet  can  be  used  to  evaluate 
the  stability  of  a  three-dimensional  wedge  of  rock  resting  on  planes 
having  frictional  resistances.  The  method  is  very  similar  to  the  two- 
dimensional  graphical  force  polygon  used  to  sum  forces.  However,  only 
the  orientation  (and  not  the  magnitude)  of  forces  is  determined  directly 
from  the  stereonet.  If  the  resultant  driving  force  acts  at  an  angle 
further  from  the  normal  to  the  potential  failure  planes  than  the  angle  of 
the  maximum  resisting  reaction  on  the  planes,  then  sliding  will  occur. 
Note  that  the  location  of  the  forces  and  reactions  is  not  known,  and 
a  summation  of  moments  is  not  carried  out. 

The  stability  analysis  is  divided  into  two  distinct  parts.  In  the 
first  part  the  orientation  of  the  maximum  resisting  reaction  on  the 
potential  failure  planes  is  plotted  on  the  stereonet.  (For  sliding  on  a 
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single  plane,  the  maximum  reaction  would  be  oriented  at  0  degrees  to 
the  normal  of  the  plane.)  Zones  of  stability  and  instability  can  thus 
be  outlined  on  the  stereonet,  strictly  by  considering  the  orientation 
of  the  reactions  on  the  potential  sliding  planes. 

The  second  part  involves  determination  of  the  orientation  of  the 
resultant  driving  force  acting  on  the  wedge.  This  force  may  include 
the  weight  of  the  wedge  as  well  as  acceleration  forces,  uplift  water 
pressures  on  the  planes  of  failure,  and  driving  forces  on  the  wedge 
from  structures  such  as  dam  abutments.  Graphical  addition  of  vectors 
is  used  in  conjunction  with  the  stereonet  to  determine  the  orientation 
of  the  resultant  vector  force.  If  the  orientation  of  the  resultant 
driving  force  falls  within  the  zone  of  stability  on  the  stereogram, 
then  the  wedge  is  stable;  if  the  orientation  of  the  resultant  driving 
force  lies  outside  the  stable  zone,  then  the  wedge  is  unstable. 

Not  only  is  the  stereonet  method  of  evaluating  the  stability  of 
a  wedge  simple  and  rapid,  it  also  possesses  the  advantage  that  a  variety 
of  forces  required  to  cause  failure  or,  conversely,  to  ensure  stability 
can  be  clearly  visualized,  without  resorting  to  extensive  computations. 
4,3  Sliding  on  a  Single  Fricitional  Plane 

The  simple  case  of  sliding  on  a  single  plane  is  described,  to  illus¬ 
trate  the  use  of  the  stereonet  in  stability  analysis.  Of  course,  a  true 
two-dimensional  problem  (where  the  resultant  driving  vector  force,  R, 
acts  in  the  direction  of  the  dip)  is  more  simply  solved  using  a  conven¬ 
tional  two-dimensional  force  polygon.  However,  for  cases  where  the 
driving  vector  is  not  in  the  direction  of  the  dip  (such  as  might  occur 
when  an  abutment  load  acts  on  a  wedge),  the  stereographic  method  can 
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be  used  to  solve  problems  which  cannot  be  readily  solved  using  a  two- 
dimensional  force  polygon. 

4.3.1  Orientation  of  reaction  force  on  the  plane  of  failure 

The  reaction  force  at  failure,  R^  (summation  of  the  normal  force, 
N,  and  maximum  shear  force,  S)  is  oriented  at  the  angle  of  friction  0, 
from  the  normal  to  the  plane.  Should  the  tendency  for  sliding  be  down- 
dip,  then  S  acts  upslope  and  R^  is  as  shown  in  Fig.  4.6a.  A  friction 
cone  can  be  drawn  to  show  the  possible  orientations  of  R^  for  sliding 
in  other  directions.  The  sides  of  the  cone  are  oriented  at  0  degrees 
to  the  normal,  as  shown  in  Fig.  4.6a  and  b.  As  long  as  the  resultant 

driving  vector,  R,  acts  at  an  angle  less  than  0  degrees  to  the  normal, 

then  sliding  will  not  occur  in  any  direction.  When  R  =  R^,  sliding  is 
initiated. 

A  friction  cone  will  plot  as  a  circle  on  an  equal  angle  stereonet, 
as  shown  in  Fig.  4.6c.  The  position  of  the  normal  force  is  first 
located  on  the  stereonet.  (The  position  of  the  normal  force  is  located 

at  the  pole  of  the  plane.)  The  friction  circle  can  then  be  drawn  by 

marking  off  40  degree  angles  from  N,  on  great  circles  passing  through 
N.  (Note  that  N  is  not  in  the  center  of  the  circle  formed  by  the 
friction  cone.) 

4.3.2  Stability  of  wedge  of  weight  W  with  uplift  force,  U,  acting 
on  the  failure  plane 

It  is  immediately  apparent  that  a  wedge  of  weight  W  will  not  s'ide 
on  the  plane  of  failure  because  0  i s  40  degrees  and  exceeds  the  slope 
angle  of  30  degrees.  This  is  also  apparent  from  Fig.  4.6c,  where  the 
weight  vector,  W,  falls  within  the  friction  cone. 
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If  a  porewater  pressure  were  acting  on  the  plane  of  failure,  the 
stability  of  the  wedge  would  be  reduced.  The  porewater  vector  force,  U, 
acts  normal  to  the  plane  of  failure,  as  shown  in  Fig.  4.6a.  The  resultant 
driving  vector,  R  =  W  +  0,  can  be  determined  by  drawing  the  two  vectors  to 
scale  (see  Fig.  4.6a)  and  determining  the  angle  of  the  resultant.  In 
this  case,  the  magnitude  of  U  is  given  as  0.44  W  and  therefore  the  angle 
of  R  from  the  vertical  is  found  to  be  20  degrees.  R  is  thus  located 
10°  outside  the  friction  circle,  in  the  unstable  zone. 

The  factor  of  safety  for  the  two  cases,  with  and  without  the  uplift 
force  acting,  is  shown  in  Fig.  4.6d.  The  tangent  of  the  angle  between 
the  normal  and  the  resultant  driving  force  determines  the  denominator 
in  each  case. 

4.3.3  Graphical  procedure  for  determining  the  direction  of  resultant 
vector  force 

The  summation  of  a  series  of  vectors  cannot  be  performed  using  the 
stereographic  projection  alone,  because  there  is  no  method  for  showing 
magnitudes  of  forces  on  the  stereographic  projection.  However,  the 
orientation  of  the  resultant  vector  can  be  determined  using  the  stereo¬ 
graphic  projection  in  combination  with  the  graphical  addition  of  vectors, 
two  at  a  time.  Three  vectors,  W  and  U  of  the  preceding  example  and  an 
additional  vector  A, are  illustrated  in  Fig.  4.7.  The  graphical  addition 
of  these  vectors  is  performed  as  shown  in  Fig.  4.8.  As  described  in 
the  preceding  example,  vectors  W  and  U  are  added  graphically  thus  deter¬ 
mining  the  orientation  of  W  +  U,  which  is  found  to  be  20  degrees  from 
the  vertical  (Fig.  4.8a).  Vectors  W  and  A  are  then  added,  determining 
the  orientation  of  W  +  A,  30  degrees  from  the  vertical  (Fig.  4.8b). 


FIG.  4.7  GRAPHICAL  SUMMATION  OF  3  VECTORS 


GIVEN:  I.  Weight  of  Wedge  *_  w 

2.  Porewater  Foret,  U  ,  Equal  to  0.44  W,  on  30°  Plane 

3.  Force  A  ,  Equal  to  0.6  W,  Acting  S  45°  W.  Dip  10° 


(a)  Vortical  Plant  Oritnttd  N-S 
Containing  Vectors  W  and  U 


(b)  Vertical  Plant  Oritnttd 
S  45°  W  Containing  Vectors 
W  ond  A 


FIG.  4.8  GRAPHICAL  DETERM  I  NATION_OF_  ORIENTATION  OF 
RESULTANT  VECTOR  ,  W  -t-  A  +  U 
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The  orientations  of  vectors  W,  U,  A,  W  +  0,  and  W  +  A  are  then  plotted 
on  the  stereogram  (Fig.  4.8c,  solid  arcs). 

Once  these  vectors  have  been  plotted,  the  orientation  of  the 
resultant  vector,  R  =  W  +  A  +  U,  can  be  found  using  only  the  stereogram. 
This  is  accomplished  by  finding  the  line  of  intersection  of  two  planes. 
One  plane  contains  W  +  A  and  0,  the  other  contains  W  +  0  and  S  (Fig.  4.7) 
The  intersection  of  the  two  planes  is  the  resultant  vector,  R  «  W  +  K  +  0 
On  the  stereogram  in  Figure  4.8b,  a  great  circle  is  drawn  through  W  +  A 
and  U,  another  great  circle  is  drawn  through  W  +  0  and  A.  The  two 
great  circles  intersect  at  R  =  W  +  0  +  A,  which  is  thus  determined  as 
dipping  42  degrees  from  the  horizontal  in  a  direction  of  S  32°  W. 

4.3.4  Determination  of  direction  of  movement  and  factor  of  safety  for 
case  of  resultant  driving  vector,  W  +  U  +  A,  acting  on  the  wedge 
In  Fig.  4.9  the  resultant  driving  vector,  R  =  Q  +  0  +A  has  been 
combined  with  the  friction  cone  diagram.  R  =  W  +  A  +  U  lies  outside  of 
the  friction  cone,  therefore  sliding  of  the  wedge  will  occur.  The 
direction  of  sliding  on  the  plane  will  be  in  the  direction  of  the  shear 
force,  S.  Sliding  is  along  a  line  plunging  25°  in  a  S  27°  W  direction 
(down  an  apparent  dip  slope).  Note  that  this  direction  is  not  the  same 
as  the  S  32°  W  direction  of  the  resultant  driving  vector,  R. 

The  factor  of  safety  is  determined  from  the  angular  distances  along 
this  great  circle.  From  N  to  R^,  the  angle  is  40  degrees,  while  from 
N  to  R  the  angle  is  75  degrees.  The  factor  of  safety  is  therefore: 

tan  40°  =  q. 22 
tan  75° 
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4.3.5  Minimum  force  NW  required  to  cause  failure 


The  orientation  of  the  minimum  force,  NW,  required  to  cause  failure 
on  an  otherwise  stable  slope  can  be  rapidly  determined  from  the  stereonet. 
To  determine  the  magnitude  of  the  minimum  force,  one  auxiliary  graphical 
construction  is  required  (Fig.  4.10). 

For  the  single  30°  plane  shown  in  Fig.  4.6,  the  wedge  is  stable 
under  its  own  weight,  W.  To  reduce  the  factor  of  safety  to  unity,  the 
angle  between  the  weight,  W,  and  the  limiting  reaction,  R^  must  be 
closed.  The  minimum  angle  is  10  degrees  and  will  be  obtained  when 
the  driving  forces  cause  sliding  directly  down-dip  (to  the  south). 

Any  other  direction  of  sliding  will  result  in  a  larger  angle  between 
R^  and  W  and  therefore  a  larger  value  for  NW. 

The  minimum  force,  NW  will  therefore  be  directed  to  the  south  and 
will  be  directed  upward  10  degrees  so  that  it  is  normal  to  R^  (Fig.  4.10). 
The  minimum  force  will  be  almost  horizontal  for  the  case  of  frictional 
sliding  on  a  wedge  loaded  only  by  its  own  weight,  where  the  factor  of 
safety  is  near  unity. 

4.4  Sliding  on  Two  Frictional  Planes 
4.4.1  General 

The  possible  modes  of  failure  of  a  wedge  on  two  planes  can  be  rapidly 
determined  from  the  stereonet.  The  orientation  of  the  driving  forces 
determines  whether  sliding  along  the  line  of  intersection  of  the  planes 
or  sliding  on  either  one  of  the  planes  will  occur.  An  example  problem 
for  sliding  on  two  planes  has  been  used  to  clarify  the  following  dis¬ 
cussion.  The  problem  is  illustrated  in  Figs.  4.11  through  4.15. 


TOO 
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4.4.2  Orientation  of  line  of  intersection  of  the  two  planes 

The  orientation  of  the  line  of  intersection  of  two  potential 
failure  planes  is  determined  using  the  stereonet  as  illustrated  in 
Fig.  4.11.  The  great  circles  for  the  two  planes  are  drawn  on  the  stereo- 
net  and  their  intersection  is  determined  as  described  in  section  4.1. 

For  the  example  problem  illustrated  in  Fig.  4.11,  the  line  of  intersection 
is  oriented  S  27°  W  and  plunges  40  degrees  from  the  horizontal. 

Figure  4.12  is  a  block  diagram  of  the  two  planes,  showing  their 
line  of  intersection  and  the  friction  cones  acting  on  each  plane. 

For  convenience,  the  friction  cones  are  shown  above  the  sliding  plane. 

4.4.3  Reaction  forces  on  the  failure  planes 

Stable  and  unstable  zones  are  separated  on  the  stereonet  (Fig.  4.13) 
by  the  limiting  reaction  forces,  R.  and  R,  .  The  unstable  zones  include 

4  4 

zones  for  sliding  down  the  intersection,  sliding  up  the  intersection, 
sliding  on  single  planes,  and  lifting  of  the  wedge  off  the  planes.  For 
the  case  of  sliding  on  plane  1  alone,  the  orientation  of  R.  ,  as  defined 

4 

by  the  friction  cone  on  plane  1,  separates  the  stable  and  unstable  zones. 
For  sliding  along  the  intersection  of  planes  1  and  2,  the  orientation 
of  R.  +  R.  separates  the  stable  and  unstable  zones.  The  boundary  be- 

4  4 

tween  sliding  on  the  intersection  and  sliding  on  plane  1  is  the  great 
circle  which  passes  through  and  Sj,  the  normal  and  shear  forces, 
respectively,  on  plane  1.  This  great  circle  represents  a  plane  normal 
to  plane  1  and  parallel  to  the  line  of  intersection. 

4.4.4  Method  of  locating  boundary  between  stable  and  unstable  zones 
for  the  case  of  sliding  along  the  line  of  intersection 

The  location  of  the  resultant,  R,  +  R.  ,  must  be  determined  in 

4  4 


Plane  I 


Given  : 

Plane  I  :  Strike  N54°E,  Dip  62°SEI**20‘> 
Plane  2  :  Strike  N04°W,  Dip  59° SW,  *‘40° 


FIG.  4.11  SLIDING  ON  TWO  PLANES:  ORIENTATION  OF 
LINE  OF  INTERSECTION 
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order  to  outline  the  stable  and  unstable  zones  for  the  case  of  sliding 
along  the  line  of  intersection.  The  shear  forces  on  planes  1  and  2 
will  act  in  the  direction  of  sliding,  which  in  this  case  is  parallel  to 
the  line  of  intersection.  Therefore,  the  shear  forces,  Sj  and  S^.  will 
both  plot  on  the  stereogram  at  the  same  point  as  does  the  line  of  inter¬ 
section  (point  Sj,  §2  in  Fig.  4.13). 

The  direction  of  the  reaction  force  on  each  plane  is  known,  since  the 
direction  of  its  components,  the  normal  and  shear  forces  on  that  plane, 

are  fixed.  Tht  reaction  force,  R.  ,  must  act  within  the  plane  in  which 

1 

N.  and  S.  act.  Therefore,  the  direction  of  the  reaction  force,  R.  ,  can 

.  _  _  1 

be  located  by  drawing  a  great  circle  through  N^  and  Sp  R^  is  located 

where  this  circle  intersects  the  friction  cone  of  plane  1.  Similarly, 

Rl  is  located  where  the  great  circle  drawn  through  ii2  and  §2  intersects 

the  friction  cone  of  plane  2.  No  matter  what  driving  forces  act  on  the 

wedge,  as  long  as  the  limiting  case  of  sliding  along  the  intersection  is 

considered,  then  the  orientation  of  both  R.  and  R  are  fixed. 

Li  _  2 

If  R.  and  R.  are  summed,  their  resultant,  R.  +  R  ,  must  act  in 
1  U2  _  L1  L2 
a  plane  parallel  to  R.  and  R.  .  This  plane  can  be  located  on  the 

1  2 

stereonet  of  Fig.  4.13  by  drawing  a  great  circle  through  R.  and  R  . 

1  z 

For  sliding  along  the  intersection  of  plane  1  and  2,  the  reaction 

R.  +  R.  will  be  located  somewhere  on  this  great  circle,  but  its  position 

L1  L2 

along  the  great  circle  will  depend  on  the  orientation  of  the  driving 

forces  acting  on  the  wedge,  since  the  orientation  of  the  driving  vector 

affects  the  relative  magnitudes  of  R.  and  R.  .  Should  the  orientation 

1  2 

of  the  resultant  vector  forces  lie  outside  the  stable  zone  outlined  by 

R.  +  R  ,  then  sliding  will  occur  on  the  intersection  of  the  two  planes. 
L1  L2 
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In  Fig.  4.13,  the  weight  vector,  W,  is  located  just  within  the  stable 
zone.  Only  a  very  small  force  directed  toward  the  south  would  be  required 
to  move  the  driving  vector  out  of  the  stable  zone  and  cause  sliding  down 
the  line  of  intersection  of  planes  1  and  2. 

4.4.S  Minimum  force  (Ny)  required  to  cause  sliding  of  the  wedge 

In  order  to  cause  sliding  of  the  wedge,  the  resultant  driving  vector 
must  lie  outside  the  stable  zone.  The  minimum  force,  NW  required  to 
cause  sliding  can  be  determined  by  means  which  is  directly  analogous  to 
the  method  for  determining  the  minimum  force  for  sliding  on  a  single 
plane  (Refer  to  Section  4.3.4).  To  close  the  force  polygon  (and  obtain 
a  factor  of  safety  of  one)  a  force  must  be  added  which  connects  the  tip  of 
the  existing  vector  (Weight,  W,  in  this  case)  to  the  plane  of  the  re¬ 
action,  R.  +  R.  .  The  minimum  force  will  be  the  one  acting  normal  to 
1  2 

the  plane  of  R.  +  R.  as  shown  in  Fig.  4,14. 

4  L2 

The  or ien  tat  ion  of  the  minimum  force  can  be  determined  from  the 

stereonet.  Its  magn i tude  can  be  determined  by  graphical  construction 

of  the  force  polygon  (such  as  Fig.  4.14).  In  Fig.  4.14,  the  minimum 

angle  between  W  and  R.  +  R.  ,  which  must  be  closed  for  a  factor  of 

1  2 

safety  of  one,  is  4  degrees.  The  minimum  force  is  also  directed  upward 
(in  this  case  at  an  angle  of  4°)  in  order  to  intersect  the  plane  of 
R.  +  R.  at  right  angles.  Note  that  the  strike  of  the  minimum  force 

L1  4 

(in  this  case  S40W)  is  not  the  same  as  the  strike  of  the  line  of  inter¬ 
section  (in  this  case  S27W) .  In  general,  the  minimum  force,  NW,  will 
not  be  oriented  directly  along  the  strike  of  the  line  of  intersection 
unless  the  wedge  is  acted  on  only  by  its  own  weight  and  the  friction 
angle  or.  the  two  planes  are  the  same. 
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of  Planes  I  and  2  I  NW  *  07W,  is  Directed 

S  40*  W  and  Acts  Upward  4° 
from  the  Horizontal,  so  thot_it 
is  Normol  to  the  Plane  of  RL|+  R 


FIG. 4.14  SLIDING  ON  TWO  PLANES:  MINIMUM  FORCE 
REQUIRED  TO  CAUSE  SLIDING 


It  should  be  noted  that  the  difference  in  magnitude  between  the 
minimum  force,  NW  ,  and  the  force  NW  directed  horizontally  and  parallel 
to  the  strike  of  the  line  of  intersection  will  be  very  small  in  most 
cases,  particularly  for  cases  when  the  primary  force  acting  on  the 
Wedge  is  its  own  weight,  and  the  factor  of  safety  is  only  slightly 
greater  than  one.  In  these  cases,  a  reasonable  (but  slightly  unconserv¬ 
ative)  approximation  is  that  NW  acts  horizontally,  parallel  to  the 
strike  of  the  line  of  intersection. 

4.4.6  Factor  of  safety  and  minimum  forces  required  to  stabilize  the  wedge 
Two  separate  conditions  exist  for  determination  of  the  factor  of 
safety  of  the  wedge  and  the  forces  required  to  stabilize  the  wedge.  Con¬ 
sider  the  two  conditions  illustrated  in  Fig.  4.15:  case  1,  where  the 
wedge  is  acted  on  by  a  driving  force,  D,  causing  sliding  on  a  single 
plane  (plane  I),  and  case  2,  where  the  wedge  is  acted  on  by  a  driving 
force,  B,  causing  sliding  along  the  line  of  intersection  of  planes 
1  and  2. 

Case  1 :  The  wedge  is  acted  upon  by  the  driving  force,  D,  and  will 

slide  on  plane  1  alone.  There  will  be  no  normal  force  on  plane  2.  In 

this  case,  the  orientations  of  both  the  normal  force,  Nj,  on  plane  1, 

and  the  driving  force,  0,  are  known,  while  the  orientation  of  the  shear 

force,  S.,  and  the  reaction,  R.  ,  on  plane  1  remain  to  be  determined. 

_  L1 

S.  and  R.  are  known  to  act  within  the  plane  of  N  and  5.  Thus,  their 

1  Lj  i 

position  is  obtained  by  drawing  a  great  circle  (solid  line  in  Fig.  4.15a) 
through  Nj  and  D,  then  locating  at  90  degrees  from  Nj  and  locating 

R.  at  0  degrees  from  N.,  along  the  great  circle.  In  the  example 

1  _  .  Q 

shown,  the  angle  between  N.  and  R.  is  0  =  20  ,  and  the  angle  between 


"'•MM*.-'. 


R,  and  D  is  37°. 
L1 


Therefore,  the  factor  of  safety  is: 


Maximum  shear  force  available  tan  20° _ 

Actual  shear  force  mobilized  tan  (20°  +  37°) 


0.24 


(Refer  to  Fig.  4.15a.)  Sliding  of  the  wedge  will  be  in  the  direction  of 
the  shear  force,  S,  in  this  case  plunging  19°  in  a  direction  S  40°  W, 
on  plane  1 . 

The  magnitude  of  the  minimum  force,  P,  required  to  close  the  37° 

angle  between  R.  and  5  (and  thereby  increase  the  factor  of  safety  to 

1 

one)  can  be  determined  from  the  graphical  construction  in  Fig.  4.15a. 

If  the  magnitude  of  D  is  known,  then  the  minimum  force  is: 

P  .  =  D  sin  37° 

m  i  n 

Case  2:  The  wedge  is  acted  upon  by  the  driving  vector,  B,  and 

will  slide  on  the  intersection  of  planes  1  and  2.  The  direction  of  the 

shear  forces,  Sj  and  are  fixed  parallel  to  the  line  of  intersection 

of  planes  1  and  2,  while  the  positions  of  N,  +  N.  and  R.  +  R.  remain  to 

1  2 

be  determined.  They  cai,  be  found  by  drawing  a  great  circle  (solid  line 

in  Fig.  4.15b)  through  S, ,  S?  and  B.  R.  +  R.  is  located  at  the  inter- 

1  2 

section  of  this  great  circle  and  the  great  circle  through  R.  and  R  . 

L1  1-2 

Nj  +  Nj  is  located  at  the  intersection  of  the  great  circle  through 

Sj  +  and  B,  and  the  great  circle  through  Nj  and  f^. 

The  factor  of  safety  in  this  case  is  determined  by  the  51°  angle 

between  Nj  +  Nj  and  B  and  the  33°  angle  between  and  R^  +  R^  . 


The  factor  of  safety  is 
F.S. 


therefore : 

=  tan  (33°) 

tan  (33°  +  18°) 


tan  33° 
tan  51° 


.  53 


The  direction  of  sliding  is  along  the  line  of  intersection, 


S  27 


W, 


in 


downdi pat  40°. 

The  concept  of  a  factor  of  safety  is  somewhat  misleading  in  this 

because  the  force  required  to  stabilize  the  wedge  does  not  have 

to  close  the  18°  angle  between  R.  +  R.  and  B.  Instead,  the  minimum 

1  L 2 

force,  P  .  is  B  sin  16  .  Note  that  the  new  resultant  driving  vector, 
min 

P  .  +  B,  acts  in  a  plane  which  is  different  from  the  plane  in  which  B 

min 

originally  acted. 

4.5  Sliding  of  a  Wedge  Bounded  by  Three  Planes 

The  case  for  sliding  of  a  wedge  bounded  by  three  or  more  planes  is 
only  slightly  more  complicated  than  the  case  for  a  wedge  bounded  by  two 
planes.  With  three  planes  another  friction  circle  is  added  to  the 
stereonet.  Depending  on  the  orientation  of  the  driving  forces,  sliding 
will  occur  on  any  one  of  the  three  planes,  on  any  one  of  the  three  lines 
of  intersection,  or  the  wedge  will  lift  off  the  three  planes.  Methods 
for  determining  the  minimum  forces  to  cause  sliding,  for  determining 
factors  of  safety,  etc.,  are  identical  to  those  described  for  the  two 
plane  case  (section  4.4).  Prior  to  performing  the  stability  analysis, 
a  basic  decision  must  be  made  as  to  which  planes  are  potential  sliding 
planes  and  which  wedges  are  critical. 

Figures  4.16  and  4.17  illustrate  the  three  plane  case.  The  orienta¬ 
tions  and  friction  angles  for  the  three  planes  are  given  in  the  block 
diagram  of  Fig.  4.16.  The  corresponding  stereonet  is  illustrated  in 
Fig.  4.17.  For  this  case  it  is  readily  apparent  that,  regardless  of 
the  presence  of  the  third  plane,  the  wedge  is  still  closest  to  a  failure 
by  sliding  (under  its  own  weight)  along  the  intersection  of  planes 
1  and  2.  In  order  for  failure  to  occur  by  sliding  up  plane  3,  or  by 
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sliding  along  the  intersection  of  plane  3  and  1,  or  by  sliding  along 
the  intersection  of  plane  3  and  2,  an  appreciable  driving  force  acting 
upward  (toward  the  North)  would  be  required. 

4.6  Wedge  Bounded  by  Three  Planes  but  Dayliqhted  by  Cut  Face 

Although  the  wedge  bounded  by  planes  1,  2,  and  3  is  stable  for  the 
condition  illustrated  in  Fig.  4.17,  it  will  not  be  stable  if  plane  3 
is  daylighted  at  the  base  and  sides  of  the  cut,  as  shown  in  Fig.  4.18. 
Cases  similar  to  this  will  commonly  occur  in  rock  masses  where  joint 
sets  form  multiple  wedges,  rather  than  a  single  wedge. 

If  planes  1  and  2  are  present,  as  shown  in  Fig.  4.18a,  then  wedge 
A  would  still  be  stable  under  its  own  weight,  as  was  previously  shown  in 
Fig.  4.17.  However,  if  plane  3  is  daylighted  at  the  edge  of  the  cut, 
then  plane  2  no  longer  restrains  wedges  A  and  B,  thus  sliding  will  occur 
along  the  intersection  of  planes  1  and  3,  as  indicated  in  the  stereonet 
of  Fig.  4. 18b. 

Another  possible  mode  of  failure  would  be  for  wedges  A  and  B 
(acting  as  a  single  wedge)  to  rotate  away  f rom  the  line  of  intersection 
of  planes  1  and  3  and  slide  on  plane  3  alone.  This  is  likely  to  occur 
if  the  mass  of  the  wedges  is  concentrated  over  plane  3,  away  from  the 
line  of  intersection  of  planes  1  and  3.  (Rotational  wedge  failures 
are  analyzed  in  section  4.7.) 

Another  possibility  is  that  wedges  A  and  B  would  break  up  and  slide 
as  individual  blocks,  wedge  A  possibly  sliding  on  the  intersection  and 
wedge  B  sliding  on  plane  3  alone. 

When  the  planes  represent  joint  sets  rather  than  single  discontin¬ 
uities,  then  modes  of  failure  similar  to  those  described  above  become 
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I  Wedge  B  - 
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v  Sliding  up  ^ 
\  Intersection  of 
\  Plones  I  and  3 


Plane  3  is  Daylighted  at 
Edge  of  Block 


Sliding  Down 
Intersection  of 
Planes  I  and  3 


STABLES 


V\^N| 


Sliding  on  Plane  3 


Sliding  on  Plane  I 


FIG  4.18  WEDGE  BOUNDED  BY  THREE  PLANES!  PLANE  3 
DAYLIGHTED  AT  EDGE  OF  BLOCK 


quite  possible.  Some  of  the  joints  within  a  set  may  daylight  on  the 
sides  or  edges  of  slopes  and  allow  failure  of  wedges  on  single  planes 
rather  than  multiple  planes.  Joints  of  other  sets  may  also  free  the 
edges  of  wedges,  facilitating  sliding  on  a  single  plane.  In  general, 
a  lower  factor  of  safety  will  be  obtained  if  the  mass  has  a  tendency 
to  move  as  several  rigid  bodies  broken  up  by  the  joint  sets,  rather 
than  moving  as  a  single  rigid  body.  It  is  very  important  that  these 
possibilities  be  considered  during  the  exploratory  phase  of  the  slope 
stability  study,  prior  to  selecting  the  critical  wedges  on  which  the 
stability  analyses  will  be  performed. 

4.7  Rotation  of  Wedge  on  Plane  3 

The  factor  of  safety  against  rotation  of  a  wedge  can  be  estimated 
using  the  stereonet  if  the  point  of  application  of  the  driving  forces 
on  the  failure  plane  is  known.  Consider  the  example  of  Fig.  4.19. 

The  wedge  illustrated  here  is  the  same  as  the  wedge  of  Fig.  4.18  except 
for  an  increase  in  the  angle  of  friction  on  plane  1  from  34°  to  40°, 
so  that  the  wedge  is  now  stable  against  sliding  under  its  own  weight 
down  the  intersection  of  planes  1  and  3.  It  is  assumed  for  a  rotational 
failure  that  the  wedge  rotates  about  point  C,  at  the  base  of  the  slope. 
It  is  also  assumed  that  the  limiting  shear  and  normal  stresses  on  plane 
3  can  be  summed  as  forces  N3  and  S3 ,  acting  at  point  D  (at  the  point  of 
intersection  of  the  resultant  driving  vector,  W,  with  plane  3).  This 
implies  that  there  is  no  variation  in  the  angle  of  friction  on  plane  3. 
Based  on  these  assumptions  the  direction  of  the  shear  force,  S3>  will 
be  perpendicular  to  line  CO,  as  shown  in  Fig.  4.19. 
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Plane  1  to  CD  (  Containing  Ns  and  Sj ) 


(a) 

N 


FIG.  4.19  ROTATION  OF  WEDGE  ON  PLANE  3 


The  given  conditions  for  the  wedge  of  Fig.  4.19a  are: 

Plane  1:  Strike  N32°E,  dip  40°SE,  4  =  40° 

Plane  3:  Strike  E-W,  dip  50°S,  4  =  20° 

In  addition,  it  is  assumed  that  the  center  of  gravity  of  the  wedge  is 
known,  so  that  the  orientation  of  line  CD  is  kncwn.  Line  CD  is  the 
moment  arm  from  the  center  of  rotation,  C,  to  the  point  of  intersection 
of  the  weight  vector,  W,  with  plane  3  (point  D) .  Its  orientation  is 
assumed  to  be  S45°W  on  plane  3. 

It  is  immediately  apparent  from  the  given  conditions  on  planes  1 
and  3  that  the  wedge,  under  its  own  weight,  W,  is  stable  against  sliding 
on  plane  1  alone,  since  the  dip  of  plane  1  equals  the  angle  of  friction 
on  plane  1  (40°).  The  wedge  is  therefore  also  stable  against  sliding 
on  the  intersections  of  planes  1  and  3. 

The  wedge  is  not  stable  against  sliding  on  plane  3  alone,  since  the 
dip  of  that  plane  exceeds  its  angle  of  friction.  If  the  wedge  were  to 
extend  an  infinite  distance  to  the  right  of  the  diagram,  then  the  wedge 
would  behave  as  if  it  were  sliding  on  plane  3  alone,  since  all  of  the 
weight  of  the  wedge  would  be  over  plane  3  and  an  infinitesimal  portion 
of  the  weight  would  act  on  plane  1.  In  this  case  the  moment  arm  for 
rotation  of  the  wedge  would  extend  parallel  to  the  strike  of  plane  3 
(line  CD').  Thus  the  shear  force,  Sg,  would  act  directly  up  the  dip- 
slope  of  plane  3,  the  identical  condition  for  sliding  on  a  single  plane. 
The  factor  of  safety  against  rotational  sliding  would  be: 

r  c  RESISTING  MOMENT  CD1  (W  tan  20°)  tan  20° 

DRIVING  MOMENT  CD'  (W  tan  50  )  tan  50 

the  identical  factor  of  safety  for  sliding  directly  down  plane  3. 
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Another  extreme  would  be  the  case  where  the  weight  of  the  wedge 
was  concentrated  near  the  intersection  of  planes  1  and  3  so  that  the 
moment  arm  extended  directly  up-slope  on  plane  3  (line  CD* ' ) .  The 
limiting  shear  force,  S^,  in  this  instance  would  be  directed  along  the 
strike  of  plane  3.  The  driving  force,  W,  has  no  component  in  this 
direction,  so  that  rotation  would  not  occur,  reaardless  of  the  frictional 
resistance  on  plane  3.  The  factor  of  safet>  against  rotation  is  therefore: 

CD"  ( tan  20°) 

F.S.  =  - 1 - — ■*-  =  » 

CE' ' (tan  0°) 

An  intermediate  case,  where  line  CD  is  directed  S45°W  on  plane  3, 
is  shown  in  the  stereonet  of  fig.  4.19.  The  rotational  stability  of 
the  wedge  is  determined  as  follows:  The  great  circle  representing  plane 
3  is  drawn  on  the  stereonet,  as  well  as  the  normal  force  and  friction 
circle  for  plane  3.  Line  CD  is  located  in  its  given  direction  on  the 
stereonet  at  the  intersection  of  plane  3  and  a  vertical  plane  oriented 
S45°W.  S  is  then  located  on  plane  3,  at  a  90°  angle  from  line  CD. 

J 

A  great  circle  is  then  drawn  through  S3  and  Nj,  representing  the  plane 
containing  and  N^.  This  great  circle  will  also  be  located  90°  from 
CD,  along  a  great  circle  (in  this  case  a  straight  line)  passing  through 
CD  and  W.  The  component  of  W  which  is  normal  to  CD  is  located  at  this 
point.  The  factor  of  safety  against  rotation  is: 

F.S.  =  tan  d/tan  (angle  between  and  the  component  of  W  1  to  CD) 

_ tan  20°  =  tan  20°  _  ,-g 

tan  (20°  +  12°)  tan  32° 

Although  the  wedge  is  stable  against  sliding  along  the  intersection 
of  planes  1  and  3,  it  is  unstable  for  the  case  of  rotation  on  plane  3. 
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However,  the  factor  of  safety  against  rotation  (F.S.  =  0.58)  is  still 
appreciably  higher  than  the  factor  of  safety  against  sliding  on  plane 
3  alone  (F.S.  =  0.30). 

In  order  for  the  factor  of  safety  against  rotation  on  plane  3  to 
be  equal  to  one,  the  center  of  gravity  must  be  located  so  that  line 
CD' 1 1  is  oriented  S27°W  on  plane  3,  as  indicated  by  the  dashed  line 
tangent  to  the  friction  circle  in  the  stereonet  of  Fig.  4.19. 
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CHAPTER  FIVE 


DYNAMIC  STABILITY  OF  ROCK  SLOPES 


5.1  Introduction 

The  stability  of  rock  slopes  subjected  to  dynamic  loads  has  usually 
been  treated  as  a  pseudo-static  problem  by  engineers  assessing  the 
"dynamic  factor  of  safety."  In  this  approach,  the  dynamic  forces  on  a 
potential  sliding  mass  of  weight  W  are  assumed  to  be  equivalent  to  a 
horizontal  force  of  KW  acting  through  the  center  of  gravity  toward  the 
free  surface  of  the  slope.  The  constant  K  is  called  the  seismic 
coefficient,  the  value  of  which  is  commonly  taken  between  0.05  and  0.20 
for  earthquake  design  in  seismically  active  areas.  If  this  method  is 
used,  the  factor  of  safety  is  computed  by  the  methods  outlined  in  Chapters 
3  and  4  for  static  problems.  The  only  adjustment  which  has  to  be  made 
is  that  the  resultant  force  R  acting  on  a  three  dimensional  tetrahedron 
includes  not  only  the  weight  and  pore  pressures,  but  also  the  force  KW 
directed  in  a  horizontal  direction.  Wittke  (1965)  has  presented  this 
type  of  analysis  to  assess  the  dynamic  stability  of  slopes  and  has 
suggested  that  the  horizontal  force  KW  be  directed  along  the  line  of 
intersection  of  the  two  planes  on  which  sliding  takes  place  for  the  most 
critical  effect.  The  meaning  of  the  factor  of  safety  calculated  for 
dynamic  loading  is  somewhat  nebulous  however  because  the  dynamic  forces 
are  not  constant  static  forces  acting  in  one  direction.  The  actual 
factor  of  safety  of  a  slope  subjected  to  dynamic  loading  varies  with 
time,  and  movements  of  the  slope  only  occur  when  the  factor  of  safety 
is  momentarily  below  1.0.  Thus  the  average  factor  of  safety  computed 
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from  a  pseudo-static  analysis  using  the  seismic  coefficient  K  is  not 
meaningful  because  the  analysis  does  not  indicate  the  magnitude  of  the 
strains  or  displacements  which  may  develop  in  the  slope.  An  estimate 
of  the  displacement  or  strain  in  the  slope  after  shaking  enables  the 
engineer  to  make  a  judgment  on  whether  the  displacements  are  harmful  in 
terms  of  either  structural  damage  or  are  large  enough  to  cause  a  consid¬ 
erable  decrease  in  shear  strength  of  the  slope  materials. 

In  this  chapter  the  basic  analysis  given  by  Newmark  (1965)  for  the 
dynamic  analysis  of  earth  slopes  is  modified  to  assess  the  displacement 
rock  slopes  might  experience  under  dynamic  loadings.  A  criterion  is 
also  given  for  determining  if  the  computed  displacement  may  be  harmful 
to  the  stability  of  the  slope. 

5.2  Dynamic  Analysis  of  Rock  Slopes 

The  first  step  in  the  dynamic  analysis  of  a  rock  slope  is  to  evaluate 
the  dynamic  resistance  of  the  slope.  The  dynamic  resistance  is  defined 
as  the  minimum  force  applied  through  the  center  of  gravity  of  the 
potential  sliding  mass  which  will  just  begin  to  move  the  mass  above 
the  assumed  failure  surface.  The  dynamic  resistance  is  usually  denoted 
by  NW  (Fig.  3.4a)  where  N  is  a  coefficient  and  W  is  the  weight  of  the 
potential  sliding  mass.  Physically  the  dynamic  resistance  is  the  minimum 
shearing  resistance  which  can  be  mobilized,  in  addition  to  that  required 
for  static  stability,  to  resist  the  effects  of  a  dynamic  load. 

For  the  three  dimensional  analysis  of  rock  slopes,  as  shown  in  Figs. 


3.4  and  3.5,  the  problem  is  to  find  the  direction  which  will  minimize 
N  and  solve  for  the  magnitude  of  NW.  Graphical  and  analytical  procedures 
have  been  developed  for  finding  the  minimum  dynamic  resistance  in 


sections  3.2.2  and  3.4.2  respectively 


The  quantity  Ng  is  the  steady  acceleration  in  the  direction  of  NW 
which  will  just  overcome  the  resistance  of  the  sliding  mass.  If  the 
maximum  acceleration  Ag  in  the  area  of  the  slope  is  less  than  Ng  then 
the  slope  is  definitely  safe.  However,  if  Ag  exceeds  Ng,  the  slope 
does  not  necessarily  fail  because  the  ground  acceleration  may  only  exceed 
Ng  for  a  very  short  period  of  time.  During  this  time,  a  relative  dis¬ 
placement  occurs  between  the  portions  above  and  below  the  failure  surface. 
Figure  5.1  shows  a  plot  of  the  ground  motions  observed  during  the  1940 
El  Centro  earthquake.  Note  that  the  maximum  acceleration,  Ag  (  =  V  ) , 
of  0.32  g  occured  only  for  a  short  period  of  time  and  also  note  that  if  a 
slope  in  the  area  of  this  record  had  a  dynamic  resistance,  NW,  of  0.2  W 
that  the  dynamic  resistance  would  only  be  exceeded  for  very  short  periods 
of  time  during  6  pulses.  Although  some  relative  displacement  would 
occur  during  these  short  times  the  slope  would  not  necessarily  fail. 
Newmark' s  method  of  analysis  provides  a  means  of  calculating  the  displace¬ 
ments  which  occur  for  the  case  when  A>  N. 

The  following  example  illustrates  Newmark' s  analysis  for  the  case 
of  a  single  acceleration  pulse  acting  on  the  base  beneath  a  sliding 
mass  with  a  rigid  plastic  resistance  between  the  base  and  mass  with  a 
resistance  of  NW. 

Consider  the  rigid  body  shown  in  Fig.  5.2  having  a  weight  W,  mass 
M,  and  having  a  motion  x.  The  motion  of  the  ground  on  which  the  mass 
rests  is  designated  by  y(t),  where  y  is  a  function  of  time  t.  The 
relative  motion  of  the  mass,  compared  with  that  of  the  ground,  is  given 
by  u,  where 

u  ■  x  -  y  (5.1 ) 
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FIG  5. 1  EL  CENTRO,  CALIFORNIA,  EARTHQUAKE  OF  MAY  18,  1940, 
N-S  COMPONENT 


The  resistance  to  motion  is  accounted  for  by  a  shearing  resistance, 
which  can  be  expressed  as  being  proportional  to  the  weight  W,  and  having 
a  magnitude  of  NW.  This  shearing  resistance  corresponds  to  an  acceler¬ 
ation  of  the  ground  of  magnitude  Ng  that  would  cause  the  mass  to  move 
relative  to  the  ground. 

The  accelerating  forces  acting  on  the  mass  M  are  shown  in  Fig.  5.3. 
The  acceleration  considered  Is  a  single  pulse  of  magnitude  Ng,  lasting 
for  a  time  interval  t  .  The  resisting  acceleration,  Ag,  is  shown  by 
the  dashed  line  in  Fig.  5.3.  The  accelerating  force  lasts  only  for  the 
short  time  interval  indicated,  but  the  decelerating  force  lasts  until 
the  direction  of  motion  changes. 

In  Fig.  5.4  the  velocities  are  shown  as  a  function  of  time  for 
both  the  accelerating  force  and  the  resisting  force.  The  maximum 
velocity  for  the  accelerating  force  has  a  magnitude  V  given  by  the 
express  ion 


V  =  Agto  (5.2) 

After  the  time  t  is  reached,  the  velocity  due  to  the  accelerating  force 

remains  constant.  The  velocitiy  due  to  the  resisting  acceleration  has 

the  magnitude  Ngt.  At  a  time  t  ,  the  two  velocities  are  equal  and  the 

m 

relative  velocity  becomes  zero,  or  the  body  comes  to  rest  relative  to 
the  motion  of  the  ground.  The  formulation  for  t  is  obtained  by  equating 
the  velocity  V  to  the  quantity  Ngt,  which  results  in  the  expression 


Sn  ~  Ng 

The  maximum  displacement  of  the  mass  relative  to  the  ground, 
obtained  by  computing  the  shaded  triangular  area  in  Fig.  5.4. 


(5.3) 

,*  is 
The 
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FIG.  5.2  RIGID  BLOCK  ON  A  MOVING  SUPPORT. 
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rIG.5.3  RECTANGULAR  BLOCK  ACCELERATION  PULSE 


FIG. 5.4  VELOCITY  RESPONSE  TO  RECTANGULAR 
BLOCK  ACCELERATION 
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calculation  is  made  as  follows: 


'  ?  V  '»  -  T  V  'o 

_  i  jd  i  vi 

or  u  “  n  . *  v 

m  2  Ng  2  Ag 

hence  ^  (1  .  N}  (5.4) 

The  acceleration  pulse  shown  in  Fig.  5.3  corresponds  to  an  infinite 
ground  displacement.  During  a  real  earthquake  a  mass  would  undergo  a 
number  of  pulses  occurring  in  random  order,  some  positive  and  some 
negative  as  shown  in  Fg,  5.1.  If  we  now  consider  a  second  pulse  of  a 
negative  magnitude  that  is  sufficient  to  bring  the  velocity  to  zero 
even  without  the  resisting  force,  then  it  can  be  shown  that  the  net 
displacement  occurring  with  the  resistance  generally  cannot  exceed  that 
which  would  occur  without  the  resistance. 

The  result  of  using  Eq.  (5.4)  Is  to  generally  overestimate  the 
relative  displacement  for  an  earthquake  because  the  equation  does  not 
take  into  account  the  fact  that  the  pulses  occur  in  opposite  directions. 
However,  Eq.  (5.4)  should  give  a  reasonable  order  of  magnitude  for  the 
relative  displacement  and  it  does  indicate  that  the  displacement  is 
proportional  to  the  square  of  the  maximum  ground  velocity. 

The  result  derived  above  is  also  applicable  to  the  case  for  a  group 
of  pulses  in  which  the  resistance  in  either  direction  of  possible  motion 
is  the  same.  For  a  situation  in  which  the  body  has  a  resistance  to 
motion  greater  in  one  direction  than  in  another,  one  must  take  into 
account  the  cumulative  effect  of  the  displacements.  A  simple  example 
where  this  effect  must  be  considered  is  found  by  rotating  Fig.  5.2 
clockwise  so  that  the  body  has  a  tendency  to  slide  downhill.  In  this 
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situation,  ground  motions  in  an  upslope  direction  leave  the  mass  without 
any  additional  relative  motion  except  where  the  magnitudes  of  the  motions 
are  extremely  large.  One  may  consider  this  model  applicable  to  a  slope. 

Similar  calculations  using  wave  forms  of  4  different  earthquake 
records  were  made  for  a  resistance  of  NW  for  downhill  movement  and  an 
infinite  resistance  for  uphill  movement  on  a  digital  computer  at  the 
University  of  Illinois.  These  results  are  shown  in  Fig.  5.5.  A  con¬ 
servative  upper  bound  to  the  displacement  is  given  by 
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(5.5) 


For  *.4  a  reasonable  upper  bound  is  given  by 

,  _  vL  /.  .  H\  A 

5  “  2gN  A'  N  (5.6) 

N 

And  for  ^  <  0.2  a  reasonable  upper  bound  for. the  displacement  is  given  by 


6 


(5.7) 


The  simplified  calculation  presented  for  one  acceleration  pulse 

and  the  calculation  of  displacements  for  actual  earthquake  records  show 

that  the  slope  movements  relative  to  the  base  are  proportional  to  the 

2 

square  of  the  maximum  particle  velocity,  V  ,  at  a  given  ratio  of  N/A. 
Since  the  calculations  of  Newmark,  presented  in  Fig.  5.5,  are  for  four 
different  earthquake  records,  it  is  conservative  to  use  the  relations 
given  in  Eqs.  5.5,  5.6,  and  5.7  for  ground  motions  from  nuclear 
explosions  because  the  duration  of  shaking  is  significantly  shorter. 

5.3  Permissible  Movement  of  Rock  Slopes 

Ultimately,  the  engineer  must  decide  if  the  dynamic  displacement 
calculated  from  Eqs.  5.3  -  5.5  is  acceptable.  Many  slopes  in  soil  and 
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STANDARDIZED  MAX.  DISPLACEMENT  ,  in. 


FIG.  5.5  STANDARDIZED  DISPLACEMENT  FOR 
NORMALIZED  EARTHQUAKES . 
(UNSYMMETRICAL  RESISTANCE) 


, ,  ,v 


soft  shales  have  undergone  considerable  movement  (as  much  as  6-10  ft) 
under  earthquake  loading  without  catastrophic  consequences.  On  the  other 
hand,  there  have  been  catastrophic  failures  of  some  rock  slopes  during 
dynamic  loading,  such  as  the  Madison  Valley,  Montana  slide  and  catas¬ 
trophic  slope  failures  have  also  been  observed  in  sensitive  marine  clays 
in  Anchorage,  Alaska  under  earthquake  loading.  Jointed  rock  slopes  and 
slopes  composed  of  sensitive  marine  clays  are  similar  in  that  they  are 
composed  of  materials  which  are  strain  softening  for  displacements  be¬ 
yond  those  required  to  develop  the  maximum  shearing  resistance.  A 
qualitative  diagram  of  shearing  strength  mobilized  versus  displacement 
parallel  to  the  discontinuity  is  shown  in  Fig.  5.6b  for  rough  rock 
surfaces.  The  peak  shearing  strength  given  by  point  C  on  this  diagram 
is  given  by 

t  =  a  tan  ($  +  i)  (5.8) 

n  r 

where  i  is  the  geometrical  component  of  resistance  given  by  the  roughness 
along  the  discontinuities  and  *  is  the  residual  angle  of  shearing 
res  i  s  tance. 

With  further  displacement  the  asperities  are  sheared  off  to  a 
certain  extent  along  the  discontinuity  and  eventually  at  larger  displace¬ 
ments  the  shear  strength  will  be  reduced  to  a  value  given  by 

i  =  o  tan  i>  (5.9) 

n  r 

where  f  is  the  residual  angle  of  shearing  resistance  along  the  ei • s- 
continuity.  Thus  in  any  given  case,  the  potential  fall  off  in  strength 
between  the  peak  strength  and  the  residual  strength  is  on  tan  i,  where 
i  is  the  angle  which  the  roughness  makes  with  the  average  direction  of 
lovemen  t  alone  the  discontinuity,  as  shown  in  Fi>i.  5.6c.  The  value  of 
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FIG,  5,6  RELATIONSHIP  BETWEEN  PEAK  SHEAR  STRENGTH  AND 

THE  COMPONENT  OF  STRENGTH  DUE  TO  SURFACE  ROUGHNESS 
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^  is  a  function  of  the  type  of  rock  of  which  the  slope  is  composed  and 
can  be  relatively  easily  determined  from  smooth  samples  in  the  laboratory. 
The  selection  of  the  i  value,  however,  is  somewhat  difficult  in  that  there 
are  usually  several  groups  of  undulations  on  a  di scont i nu i t i y  which 
have  different  i  values.  For  instance,  there  may  be  broad  undulations 
(first  order  irregularities)  with  wave  lengths  on  the  order  of  8-10  ft 
which  may  have  an  i  associated  with  them  which  may  only  be  on  the  order 
of  5-15°  (Fig.  5.7).  Whereas  there  are  shorter  undulations  (second  order 
irregularities),  which  may  have  higher  i  values  (10-46°)  „s  shown  in 
Fig.  5.7.  If  the  dynamic  resistance  NW  is  calculated  on  the  basis  of 
a  peak  shearing  strength  such  as  point  C  on  diagram  5.6b,  then  it  is  of 
utmost  importance  to  know  approximately  the  wave  length  of  the  asperity 
associated  with  the  value  of  the  angle  i  chosen  in  the  analysis.  For 
instance,  if  the  value  of  i  is  associated  with  a  quarter  wave  length 
denoted  by  H  in  Fig.  5.6c,  then  it  is  obvious  that  the  dynamic  displace¬ 
ment  as  computed  by  Newmark' s  method  must  be  less  than  H  or  the  shear 
strength  value  upon  which  the  calculation  of  NW  was  computed  is  no  longer 
valid.  The  displacement  in  this  case  would  have  been  enough  to  roll  up 
and  over  the  asperity  shown  in  Fig.  5.6c  and  the  shear  strength  would 
have  been  reduced  to  some  value  lower  than  the  peak  shear  strength. 

The  full  value  of  i  would  not  be  effective  because  of  the  lower  slope 
of  the  surface  roughness  near  the  top  and  possioly  because  part  of  the 
roughness  could  have  been  sheared  off  by  the  dynamic  movement.  On  the 
other  hand,  a  va'ue  of  i  used  in  estimating  the  peak  shear  strength 
used  in  the  calculation  of  NW  could  be  a  rela>  low  value  ot  about 

r-10°;  and  this  value  of  i  could  be  <  <  ms  i  s  t  en 1  wi'1'  a  1‘ii-itli  H  (l"iq.  5.6c) 
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FIG.  5.7  AN  EXAMPLE  OF  A  DISCONTINUITY  ILLUSTRATING  FIRST 
AND  SECOND-ORDER  IRREGULARITIES 


on  the  order  of  5  or  6  ft.  In  this  case  If  the  calculated  dynamic 
displacement  caused  by  a  given  earthquake  turns  out  to  be  something 
like  6  inches,  then  we  would  definitely  say  that  the  slope  Is  probably 
stable  or  no  problems  would  result  from  dynamic  loading.  This  Is  true 
because  the  6  inch  displacement  will  not  significantly  reduce  the  strength 
assumed  in  calculations  of  NW.  The  relative  displacement  would  have  to 
be  on  the  order  of  3  or  4  ft  to  s I gn i f i can tlv  lower  the  peak  shearing 
resistance  in  this  case.  Thus,  in  general,  the  criterion  which  should 
be  used  to  decide  if  a  certain  displacement  is  detrimental  or  not  is 
the  wave  length  of  the  asperities  giving  the  geometrical  component  of 
resistance  because  this  i s  the  resistance  which  can  be  destroyed  by  the 
dynamic  displacement.  If  preliminary  calculations  Indicate  that  the 
dynamic  displacements  will  not  be  large  in  comparison  with  the  displace¬ 
ments  necessary  to  significantly  lower  the  sheer  strength,  then  the 
dynamic  displacement  calculated  will  probably  he  acceptable  and  the 
slope  can  be  judged  safe.  However,  if  the  dynamic  displacement  computed 
is  on  the  same  order  of  magnitude  as  the  wave  length  of  the  discontin¬ 
uities  or  t he  order  of  displacement  necessary  to  significantly  reduce 
the  '..hear  strength  along  a  discontinuity,  then  the  slope  may  not  be  safe. 


135 


CHAPTER  SIX 


SUMMARY  AND  CONCLUSIONS 


6 . 1  Static  Stability  of  Rock  Slopes 

In  this  report  the  methods  of  analysis  for  assessing  the  static 
stability  of  rock  slopes  in  three  dimensions  have  been  described  and  a 
method  for  calculating  the  dynamic  resistance  and  displacement  of  rock 
slopes  subjected  to  earthquake  loading  has  been  given. 

The  general  procedures  for  determining  the  static  factor  of  safety 
of  rock  slopes  analytically  by  vector  analysis  and  graphically  with 
slereonets  are  given  in  Chapter  3  and  Chapter  4  respectively.  Although 
the  mechanics  of  the  calculation  are  different  in  each  of  these  approaches, 
the  same  basic  steps  are  followed  in  each  method  for  determining  the 
static  factor  of  safety.  The  steps  in  the  analysis  are  as  follows: 

(1)  the  intersection  of  the  various  joint  sets  with  each  other 
and  with  the  slope  face  must  be  inspected  to  determine  the 
tetraherdra  which  may  be  potential  failure  wedges.  These  wedges 
must  then  be  analyzed  in  detail. 

(2)  the  forces  tendinn  to  disturb  the  equilibrium  of  the  wedge 
should  be  added  vectorially  to  give  a  resultant  driving  force. 
These  disturbing  forces  are  the  weight  of  the  wedge,  W,  the 
external  load  applied  to  Lhe  wedge  by  a  structure,  Q.,  and 

the  porewater  forces  acting  on  various  faces  of  the  tetrahedron 
given  as  Uj,  U2,  and  U^ .  This  step  is  illustrated  by  equations 
3.10  and  3.68  for  the  vector  analysis  and  is  shown  in  Figs. 

4.7  and  4.8  for  the  graphical  method  using  stereonets. 
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(3)  the  mode  of  failure  must  then  be  determined.  For  example 
a  wedge  supported  on  two  base  planes  can  either  slide  along 
the  line  of  intersection  of  the  two  planes,  slide  on  either 
plane  or  rotate  on  either  plane.  The  kinematics  of  failure 
will  depend  upon  the  orientation  of  the  disturbing  force  in 
relation  to  the  orientation  of  the  supporting  planes.  These 
kinematic  tests  for  sliding  are  illustrated  by  equations  3.16, 
3.17,  3.20,  3.21,  3.22  and  3.23  for  a  tetrahedron  supported 
on  two  base  planes.  The  kinematic  tests  for  rotation  are  given 
In  Table  3.1.  The  orientation  of  the  resultant  disturbing 
force  which  will  cause  various  modes  of  failure  for  a  wedge 
supported  on  two  base  planes  as  shown  in  Fig.  4.13  by  the 
graphical  method  using  stereonets.  Equations  3.81-3.93  are 
kinematic  tests  to  determine  the  sliding  mode  for  a  tetrahedron 
bounded  by  3  base  planes  by  means  of  vector  analysis.  Fig.  4.17 
illustrates  the  method  of  determining  the  mode  of  sliding  of 
a  tetrahedron  bounded  by  3  base  planes  by  utilization  of 
s  te  reone  ts  . 

(4)  after  the  mode  of  failure  is  determined  the  maximum  shearing 
resistance  which  can  bo  mobilized  in  the  direction  of  movement 
is  compared  to  the  shearing  forces  necessary  for  equilibrium 
to  obtain  a  factor  of  safely.  This  step  has  been  illustrated 
by  the  many  example  problems  worked  in  Chapters  3  and  4. 

The  detailed  analyses  given  in  Chapters  3  and  4  can  be  used  to 
solve  most  of  the  problems  wh'ch  arise  in  the  calculation  of  rock  slope 
stability.  A  majority  of  the  real  problems  which  arise  however  are 
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wedges  acted  on  by  their  own  weight,  partially  submerged  beneath  a 
phreatic  surface,  and  resting  on  two  base  planes.  The  influence  of 
external  loads  is  small  relative  to  the  weight  of  the  wedge.  For  this 
case,  which  has  been  the  most  common  field  case  encountered  by  the 
authors,  there  are  several  approximate  generalizations  which  can  be  made 
such  that  ail  the  details  of  analyses  presented  in  Chapters  3  and  4 
are  not  necessary  to  obtain  a  fairly  accurate  answer  to  the  problem. 
First  of  all  for  a  wedge  resting  on  two  base  planes,  acted  on  by  only 
its  own  weight,  sliding  will  occur  along  the  line  of  intersection  of 
the  two  planes  if  a  line  drawn  down  the  dip  in  both  planes  tends  to 
intersect  the  line  of  intersection.  If  in  either  one  of  the  planes  a 
line  drawn  down  the  dip  is  directed  away  from  the  line  of  intersection 
then  sliding  will  occur  on  that  plane  only  and  the  wedge  will  move 
away  from  the  line  of  intersection.  If  a  wedge  is  acted  on  by  its  own 
weight,  it  will  slide  down  the  maximum  dip  if  sliding  occurs  on  one 
plane  and  the  factor  of  safety  can  be  easily  computed.  If  it  is  deter¬ 
mined  above  that  sliding  will  take  place  along  the  line  of  intersection, 
the  slope  of  the  line  of  intersection,  .t,  as  shown  in  Fig.  6.1,  should 
be  determined  immediately  by  means  of  graphics.  The  angle  of  friction  if 
required  for  stability  will  always  be  equal  to  or  less  than  a  it  there 
are  no  pore  pressures  on  the  joint  surfaces.  fhe  next  step  is  to  de¬ 
termine  the  angle,  3,  included  between  planes  1  and  2  in  a  plane  per¬ 
pendicular  to  the  line  of  intersection  OA  as  shown  in  Fig.  6.1.  The 
smaller  anqle  6 ,  the  lower  the  value  of  <f  required  for  stability. 

As  6  approaches  zero,  the  v  lue  of  required  for  stability  approaches 
zero;  and,  as  3  approaches  180°  the  value  of  ♦  required  for  stability 


Line  of  Intersection 
Dip  of  OA 


FIG.  6.1.  FRICTION  ANGLE  REQUIRED  FOR  STABILITY  OF  WEDGE  OF  WEIGHT,  W,  FOR  C 


approaches  a.  The  next  step  is  to  determine  the  skewness  Y  of  the  wedge 

as  shown  in  Fig.  6.2.  For  the  non-skewed  or  symmetrical  case  (Cj  =  C2 

and  y  *>  0,  as  shown  in  Fig.  6.1)  the  4>  required  for  a  factor  of  1.0 

for  various  values  of  0  and  a  are  shown  in  Fig.  6.1.  If  the  angle  of 

shearing  resistance  $  is  the  same  on  both  planes  1  and  2,  the  value  of  <f> 

required  for  stability  is  less  for  the  symmetrical  case  (Cj  =  C^>  Fig.  2) 

than  for  the  skewed  case  (Cj  11  C2)  for  the  same  values  of  a  and  0. 

Figure  6.2  illustrates  the  effect  of  skewing  the  planes  on  the  value 

of  if  required  for  a  factor  of  safety  of  1.0.  From  Fig.  6.2  the  value 

of  <f  required  for  a  factor  of  safety  of  1.0  can  be  determined  from  the 

value  of  tan  <f/tan  a  for  various  values  of  3  and  y  where  Y  is  a  measure 

of  the  skewness  of  the  wedge  as  shown  in  Fig.  6.2.  The  curve  labeled 

C1  =  C2  =  °)  in  F'9-  summarizes  the  curves  presented  in  Fig.  6.1 

for  the  symmetrical  case.  The  additional  curves  presented  in  Fig.  6.2 

illustrate  the  effect  of  skewing.  These  curves  illustrate  the  sensitivity 

of  the  value  of  required  to  skewing.  For  values  of  Y  less  than  20 

the  tan  <!>  .  ./tan  a  values  are  increased  by  only  6  percent  above  the 

requ 1  red 

values  for  Y  =  0.  However  if  t'ie  wedge  is  skewed  60°  (l  =  60  )  then 

values  of  tan  <}>  .  ,/tan  a  are  approximately  twice  the  values  for 

required  rr 

o 

Y  =  0.  For  y  =  40  tan  4  .  ,/tan  o  values  are  increased  by  approx- 

requ 1  red 

imately  30  percent  from  the  case  where  y  =  0. 

The  five  example  problems  shown  below  in  Table  6.1  illustrate  all 
the  conditions  which  are  considered  in  Fig.  6.2, 


FOP  STABILITY  OF  A  WEDGE  OF  WEIGHT,  W  FOR  VARIOUS  VALUES  OF 


TABLE  6.1 

Case  1  6-180°  Cj-^-O 

Case  2  C^O 

Case  3  6-0 

Case  4  Ci“C2*  Y“° 

Case  5 


Analysis  of  each  of  these  five  cases  by  the  use  of  stereonets  is 
illustrated  in  Figs.  6.3,  6.4,  6.5,  6.6,  and  6.7. 

6.2  Dynamic  Stability 

The  analysis  of  rock  slopes  in  a  static  fashion  by  considering  the 
inertia  force  as  a  static  load  (Wittke,  19o5a)  is  not  considered  adequate 
for  assessment  of  the  dynamic  behavior  of  rock  slopes.  In  this  report 
a  method  was  given  in  sections  3.2.2  and  3.4.2  for  calculating  the  dynamic 
resistance  NW  for  a  rock  slope  by  means  of  vector  analysis.  The  calcu¬ 
lation  of  the  dynamic  resistance  by  means  of  stereonets  is  given  in 
section  4.3.4,  The  dynamic  resistance  Nrf  should  then  be  used  in  the 
dynamic  analysis  proposed  by  Newmark  (1965),  which  is  explained  and 
illustrated  in  Chapter  5.  From  this  analysis  the  dynamic  displacement 
is  computed.  This  displacement  should  then  be  compared  with  the  wave¬ 
length  of  the  asperities  on  the  failure  planes  as  shown  in  Fig.  5.6  to 
determine  if  the  dynamic  displacement  would  be  detrimental  to  the  static 
stability  of  the  slope.  The  calculation  of  the  dynamic  factor  of  safety 
using  a  pseudo  static  analysis  has  little  meaning. 
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STABLE 


Line  of  interaction 

A 


I,  =  l80°,  Cx  *  o..  «  0 
Kor  factor  of  safety  *  1 
^required  a 


MG.  6.3.  WEDGE  ACTING  UNDKK  OWN  WEIGHT 
CASE  It  SINGLE  PLANE,  6  -  180° 
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APPENDIX  A 


COMPUTER  PROGRAMS 
for  the 

STABILITY  OF  ROCK  SLOPES 


A-l 


MR  XSJMWWiWIIMlWSaM*  wksbss 


PROGRAM  #1 


This  program  is  for  the  analysis  of  a  rock  tetrahedron  supported  on 
two  planes  and  acted  on  by  its  own  weight.  The  program  is  coded  in  FORTRAN 
IV  language.  The  format  and  definition  of  parameters  are  given  below. 


INPUT  PARAMETERS 


CARD  1 

Format  (10F7.2) 

Parameter 

Card  Column 

BETA  1 

1  -  7 

BETA  2 

8-14 

GAMMA  1 

15—  21 

GAMMA  2 

22  -  28 

PHI  1 

29  -  35 

PHI  2 

36  42 

ALPHA 

43-49 

DELTA 

50  -  56 

H 

57-63 

Q 

64  -  70 

OUTPUT  PARAMETERS 


Parameter 

UX(1) 

UY(1) 

UZ(1) 

VX(1) 

VY  (1) 
VZ(1) 

WX(1) 

WY(1) 

WZ(1) 

UX(2) 

UY  (2) 
UZ(2) 

VX(2) 

VY(2) 


Definition  of  Parameter 

Strike  of  plane  1 
Strike  of  plane  2 
D i p  of  plane  1 
Dip  of  plane  2 

Angle  of  Shearing  Resistance  on  plane  1 
Angle  of  Shearing  Resistance  on  plane  2 
a.  Slope  angle  as  shown  in  Fig.  3.9 
6 t  Slope  angle  as  shown  in  Fig.  3.9 
Height  of  slope  as  shown  in  Fig.  3.9 
Desired  factor  of  safety  for  design 


Definition  of  Parameter 

x,  y,  and  z  components  of  a  unit  vector  U  in 
the  direction  of  the  strike  of  plane  1. 


x,  y,  and  z  components  of  a  unit  vector  V  in 
the  direction  of  the  dip  of  plane  1. 


x,  y,  and  z  components  of  a  unit  vector  W  in 
a  direction  normal  to  plane  1. 


x,  y,  and  z  components  of  a  unit  vector  U  in 
the  direction  of  the  strike  of  plane  2. 


x,  y,  and  z  components  of  a  unit  vector 
the  direction  of  the  dip  of  plane  2. 


<• 


WX(2) 

x,  y,  and  z  components  of  a  unit  vector  W  in 

WY  (2) 

a  direction  normal  to  plane  2. 

WZ(2) 

XX 

x,  y,  and  z  components  of  a  vector  X  directed 

XY 

along  the  line  of  intersection  of  planes  1 

XZ 

and  2 . 

X 

The  magnitude  of  the  vector  X. 

FS 

Factor  of  safety. 

Definition  of  Parameters  Calculated  in  Program  in  Intermediate  Steps 

ODX  x,  y,  and  z  components  of  line  OD,  Fig.  3.9 

ODY 

ODZ 


OCX  x,  y,  and  z  components  of  line  OC ,  Fig.  3.9 

OCY 

OCZ 

OBX  x,  y,  and  z  components  of  line  OB,  Fig.  3.9 

OBY 

OBZ 

OSX  x,  y,  and  z  components  of  line  OS,  Fig.  3.9 

OSY 

OSZ 

RX  x,  y,  and  z  components  of  resultant  force  action  on  tetrahedron 

RY  (Note:  since  this  program  is  for  a  slope  acted  on  by  its  own 

RZ  weight  RX  and  Ri  are  always  zero  and  RZ  is  taken  as  -1  unit 

since  the  magnitude  of  the  weight  does  not  affect  the  factor  of 
safety.) 

EMX  Moment  about  the  line  of  intersection  of  planes  1  and  2  of 

the  weight  of  the  tetrahedron  applied  through  the  center  of 
gravity  of  the  tetrahedron. 

EN(j)  Magnitude  of  the  normal  component  of  the  vector  R  on  plane  J 

TX(j)  x,  y,  and  z  components  of  the  tangential  component  of  the  vector 

TY(j)  R  on  plane  J 

TZ(J) 

T(j)  Magnitude  of  the  tangential  component  of  the  vector  R  on  plane  J 

SXl  x,  y,  and  z  components  of  the  vector  S . 2  shown  in  Fig.  3.5(b) 

SY1 

SZ1 

SX2  x,  y,  and  z  components  of  the  vector  S , 2  shown  in  Fig.  3.5(b) 

SY2  Z 


COMPUTER  PROGRAM  #1 


11 

12 


STABILITY  ANALYSIS  OF  ROCK  SLOPES 

DIMENSION  B(2),C(2),PHI(2),EN(2),UX(2),UY(2),UZ(2),VX(2),VY(2), 

1 VZ(2) ,WX(2) ,WY(2) ,WZ(2) ,TX(2) ,TY(2) ,TZ(2),T(2) 

READ  (5,11)  BETA1 ,BETA2,GAMMA1 ,GAMMA2,PHI1 ,PHI2, ALPHA, DELTA, H,0 
FORMAT  (10F7.2) 

WRITE  (6,12)  BETA!, BETA2.GAMMA1 ,GAMMA2,PHI1 ,PHI2, ALPHA, DELTA, H,Q 
FORMAT  (1H1,9X,8HBETA1  =  ,F7. 1 ,8H  DEGREES, 10X,8HBETA2  =,F7.1, 

18H  DEGREES/ 1  OX, 8HGAMMA1  =,F7.1,8H  DEGREES, 1  OX, 8HGAMMA2  =,F7.1, 
28H  DEGREES/ 1  OX, 8HPH II  =  ,F7.1 ,8H  DEGREES, 1  OX, 8HPHI2  =,F7.1, 
38H  DEGREES/ 1  OX, 8HALPHA  =,F7.1 ,8H  DEGREES,! OX, 8HDELTA  =,F7.1, 
48n  DEGREES/ 1 OX, 8HHEIGHT  = ,F7. 1 ,4HFEET,14x,8HQ  =,F7.2) 

8(1)  =  BETA1  /(180. 0/3. 141593) 

B(2  =  BETA2  /  180.0/3.141593) 

C(1J  =  GAMMA1/(180. 0/3. 141593) 

C  ( 2 )  =  GAMMA2/O80.0/3.141593) 

PHI(l)  =  PH11  / ( 1 80 . 0/3. 141593) 

PHI (2)  =  PHI2  /( 180.0/3. 141593) 

A  =  ALPHA  /(180. 0/3. 141593) 

0  =  DELTA  / ( 1 80 . 0/3. 141593) 

DO  101  K=1 ,2 
UX(K)  =  COS(B(K)) 

UY(K)  =  SIN(B(K) ) 

UZ(K)  =  0.0 
VX(  K) 

VY  ( K) 


=  COS(C(K))*SIN(B(K) 
=-COS(C(K) )*COS(B(K) 


101 

13 


14 


ffi*VZ(K))-^UZ(K)*VY(K)) 
*VX(K))-(UX(K)*VZ(K)) 


VZ(K)  =-SIN(C(K)) 

WX(K)  . 

WY(K)  . . . .  . . 

WZ(K)  =  (UX(K)*VY(K))-(UY(K)*VX(K)) 

WRITE  (6,13)  UX(1 ) ,UY(1 ) ,UZ(1 ) ,VX(1 ) ,VY(1 ) ,VZ(1 ) ,WX( 1 ) ,WY( 1 ) ,WZ(1 ) 
FORMAT (lH0,9X,4HUl  =,F6.3,1H,,2X,F6.3,1H,,2X,F6.3/ 

1  10X.4HV1  =,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3/ 

2  10X.4HW1  =,F6.3,1H,,2X,F6.3.1H,,2X,F6.3) 

WRITE  (6,14)  UX(2) ,UY(2) ,UZ(2) ,VX(2) ,VY(2) ,VZ(2) ,WX(2) ,WY(2) ,WZ(2) 
FORMAT (1 HO, 9X.4HU2  =,F6.3,1H,,2X,F6.3,1H, ,2X,F6.3/ 


1  10X.4HV2  =,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3/ 

2  10X.4HW2  =,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3) 

XX  =  (WY(2)*WZ(1)-WZ  2)*WY(1)) 

XY  =  (WZ(2)*WX(1  -WX(2)*WZ(1)) 

XZ  =  (WX(2)*WY(1)-WY(2)*WX(1)) 

X  =  S0RT(XX**2  +XY**2  +  XZ**2) 

WRITE  (6,15)  XX.XY.XZ.X 

15  FORMAT (l HO, 9X,4HX  =,F6. 3,1H, ,2X ,F6.3,1H, ,2X,F6.3,5X,8HABS(X)  =, 
1F6.3) 


ODX 

1 

ODY 

ODZ 

OCX 

1 

OCY 

OCZ 


1 


1siN°A)*SIN(B(l  )^j  )(1  ^*SIN(B(1  ^ )WH*COS(A)*COS(B(  1 ))/ 


H*COS(A)/SIN(A) 

H 


fSTS?IS£|f^^i))S  > (2)  )*SIN<B<  2)  ”  5+<H*C0S(''’’ 

H*COS(A)/SIN(A) 

H 


COS ( B ( 2 ) ) / 


AAA  =(SIN(A)/C0S(A)KXZ/XY) 

866  =(XZ/XY)-(SIN(D)/COS(D) ) 

FFF=(SIN(D)*COS(A) )/(COS(D)*SIN(A) ) 

HH  =  H*AAA*FFF/DBb 

CCC  =(H+HH)/XZ 

OBX  =  CCC*XX 

OBY  =  CCC*XY 

nR7  =  rrr*x7 

OSX  = (OBX+OC  X+ODX ) / 4 , 0 

OSY  = (OBY+OCY  i-ODY )  /4 . 0 

OSZ  =(0BZ+0C7+ODZ)/4.0 

RX  =  0.0 

RY  =  0.0 

RZ  =-1.0 

EMX  =  XX*(0 YY  R/’-US/1,i!Y)+XY*(OSZ*RX--OSX*R?)iXZ*(OSX*RY-0$Y*RX) 
IF  (EMX)  102,ll)3,M4 

104  WRITE(6,105) 

105  FORMAT ( 1 1 10 , 9X , 29HRESULTANT  INTERSECTS  PLANE  El) 

J  =  1 

110  EN(J)=  RX*WX(.1)+RY*WY(J)  :RZ*WZ(J) 

TX(J)=  RX-EN(J)*HX(J) 

TY(J)=  RY-EN(J)AWY(.l) 

TZ(J)=  RZ-EN(J)*WZ(.l) 

T(J)  =  SQRT(TX(J)**2  +TY(J)**2  +TZ(J)**2  ) 

FS  =  EN(J)*SIN(PHI(J))/(T(.i)*COS(PHI(J))) 

IF  (FS.LT.Q)  GO  TO  106 

108  WRITE(6,107)  FS,J 

107  FORMAT ( 1 HO, 9X , 4HFS  =,F5.2,10X,25HSLIDING  OCCURS  ON  PLANE  E,I1) 

GO  TO  200 

106  IF(J.EQ.l)  GO  TO  111 
GO  TO  112 

111  SX1  =  (XY*WZ(1 )-XZ*WY(l ) )/X 
SY1  =  (XZ*WX(1)-XX*WZ(1))/X 
SZ1  =  (XX*WY(1 )-XY*WX( 1 ) )/X 

DDD1=  TX ( 1 )*SX1+TY(1 )*SY1+TZ(1 )*SZ1 
IF(DDDl.LE.O.O)  GO  TO  108 
GO  TO  103 

102  WRITE  (6,109) 

109  FORMAT ( 1  HO , 9X , 291 IRESULTANT  INTERSECTS  PLANE  E2) 

J  =  2 

GO  TO  no 

112  SX2  =-(XY*WZ(2)-XZ*WY(2) )/X 
SY2  =-(XZ*WX(2)-XX*WZ(2))/X 
SZ2  =-(XX*WY(2)-XY*WX(2) )/X 

DDD2=  TX(2)*SX2+TY(2)*SY2+TZ(2)*SZ2 
IF  (DDD2.LE.0.0)  GO  TO  108 

103  T12  =  (RX*XX+RY*XYhRZ*XZ)/X 
EN12X=  RX-T12*XX/X 

EN12Y=  RY-T12*XY/X 
EN12Z=  RZ-T1 2*XZ/X 
EEE  =-WX(l )*WY(2)+WY( 1 )-WX(2) 

EEE1  =-ENl 2X*WY(2)+EN12Y*HX(2) 

EEE2  =+WX( 1 )*EN1 2Y-WY( 1 )*EN12X 
ENN1  +  EEE 1 /TEE 
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.Amir 


ENN2  =  EEE2/EEE 

FS-CCENN1 *SIN(PHI(1 } )/COS(PHI (1 ) ) )+(ENN2*SIN(PHI (2) )/C0S(PHI(2) ) ) ) 
1/T12 

WRITE  (6,113)  FS 

113  F0RMAT(1H0,9X,4HFS  =,F5.2,10X,45HSLIDING  OCCURS  ALONG  THE  LINE  OF 
2INTERSECTI0N) 

200  STOP 


PROGRAM  ft 2 


This  program  is  for  the  analysis  of  a  rock  tetrahedron  supported  by 
3  planes  and  acted  on  by  its  own  weight,  an  external  force,  and  porewater 
forces  on  each  of  the  three  support  planes.  The  program  is  coded  in 
FORTRAN  IV  language.  The  format  and  definitions  of  parameters  are  given 
below. 


INPUT  PARAMETERS 


CARD  1 

FORMAT  (9F5.1) 

Parameter 

Card  Colur 

BETA  1 

1  -  5 

BETA  2 

6-10 

BETA  3 

11-15 

GAMMA  1 

16  -  20 

GAMMA  2 

21  -  25 

GAMMA  3 

26  -  30 

PHI  1 

31  -  35 

PHI  2 

36  -  40 

PHI  3 

41  -  45 

CARD  2 

FORMAT  (7F7.2) 

WT 

1  -  7 

OX 

8-14 

QY 

15  -  21 

« 

22  -  28 

UP  1 

29  -  35 

UP  2 

36  -  42 

UP  3 

43  -  49 

Definition  of  Parameter 

Strike  of  plane  1 
Strike  of  plane  2 
Strike  of  plane  3 
Dip  of  plane  1 
Dip  of  plane  2 
Dip  of  plane  3 

Angle  of  Shearing  Resistance  on  plane  1 
Angle  of  Shearing  Resistance  on  plane  2 
Angle  of  Shearing  Resistance  on  plane  3 


Weight  of  the  rock  tetrahedron  being  analyzed 
z,  y  and  z  components  of  a  force  Q  applied 
to  the  tetrahedron 

Magnitude  of  the  porewater  force  on  plane  1 

Magnitude  of  the  porewater  force  on  plane  2 

Magnitude  of  the  porewater  force  on  plane  3 


OUTPUT  PARAMETERS 

Pa rameter  Definition  of  Parameter 

UX(1)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction  of 

UY(1)  the  strike  of  plane  1 

UZ  (1) 


VX(1)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction  of 

VY(1)  the  dip  of  plane  1 

VZ(1) 


x,  y,and  z  components  of  a  unit  vector  in  the  direction  normal 
to  plane  1 
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WX(1) 

WY(1) 

WZ(1) 


UX  (2) 
UY(2) 
UZ(2) 


x,  y,  and  z  components  of  a  unit  vector  in  the  direction  of 
the  strike  of  plane  2 


VX(2)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction  of 

VY(2)  the  dip  of  plane  2 

VZ  (2) 

WX(2)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction 

WY(2)  normal  to  plane  2 

WZ(2) 

UX(3)  x,  y  and  z  components  of  a  unit  vector  in  the  direction  of 

UY(3)  the  strike  of  plane  3 

UZ(3  ) 


VX(3)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction 

VY(3)  of  the  dip  of  plane  3 

VZ(3) 

WX(3)  x,  y,  and  z  components  of  a  unit  vector  in  the  direction 

WY(3)  normal  to  plane  3 

WZ(3) 


X12X 

X12Y 

X12Z 

Xl2 

X23X 

X23Y 

X23Z 

X23 


x,  y,  and  z  components  of  a  vector,  X^>  in  the  direction 
of  the  line  of  intersection  of  planes1  1  and  2  as  shown 
in  Fig.  3.14 

the  magnitude  of  the  vector  *12 


x,  y,  and  z  components  of  a  vector,  X23,  in  the  direction  of 
the  line  of  intersection  of  planes  2  and  3  as  shown  in  Fig.  3.14 


the  magnitude  of  the  vector  X23 


X31X 

X31Y 

X31Z 

X31 

S121X 
S 1 2 1 Y 
S121Z 

S122X 

S122Y 

S122Z 

S232X 

S232Y 

S232Z 


x,  y,  and  z  components  of  a  vector,  X31  in  the  direction  of  the 
line  of  intersection  of  planes  1  and  3  as  shown  in  Fig.  3.14 

the  magnitude  of  the  vector  *31 

x,  y,  and  z  components  of  vector  jSjg  shown  in  Fig.  3.14 
x,  y,  and  z  components  of  vector  2^12  stlown  'n  Fig.  3.14 


x,  y,  and  z  components  of  vector  shown 


in  Fig.  3.14 
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S233X 

S233Y 

S233Z 

S313X 
S3 13Y 
S313Z 

S3 1 IX 
S3 11Y 
S3 11Z 


x,  y,  and  z  components  of  vector  -S„-  shown  in  Fig.  3.14 


x,  y,  and  z  components  of  vector  shown  in  Fig.  3.14 


x,  y,  and  z  components  of  vector  shown  in  Fig.  3.14 


(Sliding 


(Sliding 


(SI iding 


(SI iding 


(Sliding 


(SI iding 


tends  to  occur  along  X12  on  planes  1  and  2) 

Factor  of  Safety 
Normal  force  on  plane  1 

Normal  force  on  plane  2 

Component  of  disturbing  forces  parallel  to  the  line  of  inter 
section  of  planes  1  and  2 

tends  to  occur  along  X23  on  planes  2  and  3) 

Factor  of  Safety 
Normal  force  on  plane  2 

Normal  force  on  plane  3 

Component  of  disturbing  forces  parallel  to  the  line  of  inter 
section  of  planes  2  and  3 

tends  to  occur  along  X31  on  planes  1  and  3) 

Factor  of  Safety 
Normal  force  on  plane  3 

Normal  force  on  plane  1 

Component  of  disturbing  forces  parallel  to  the  line  of  inter 
section  of  planes  1  and  3 

tends  to  occur  on  plane  1) 

Factor  of  Safety 
Normal  force  on  plane  1 
Tangential  force  on  plane  1 

tends  to  occur  on  plane  2) 

Factor  of  Safety 
Normal  force  on  plane  2 
Tangential  force  on  plane  2 

tends  to  occur  on  plane  3) 

Factor  of  safety 
Normal  Force  on  plane  3 
Tangential  force  on  plane  3 


Note:  Tests  1  through  12  in  program  ft 2  are  used  to  determine  the  mode  of 
sliding  as  can  be  determined  from  equs.  3.69,  3.70,  3.71,  3.81, 
3.82,  3.83,  3.84,  3.85,  3.86,  3.87,  3.88,  3.89,  3.90,  3.91,  3.92, 
3.93,  3.94,  3.95,  3.96,  3.97,  and  3.98. 


COMPUTER  PROGRAM  #  2 


C  STABILITY  ANALYSIS  OF  ROCK  SLOPES-THREE  INTERSECTING  JOINT  SETS 

DIMENSION  B(3),C(3),GAMMA(3),VX(3) ,VY(3) ,VZ(3),UX(3),UY(3) , 
1UZ(3),WX(3),WY(3),WZ(3) 

DO  200  IJK=1,2 

READ  (5,11)  BETA1 ,BETA2,BETA3,GAMMA(1 ) ,GAMMA(2) ,GAMMA(3) ,PHI1 , 
1PHI2,PHI3,WT,QX,QY,QZ,UP1 ,UP2,UP3 

11  FORMAT  (9F5.1/7F7.2) 

WRITE  (6,12)  BETA1 ,BETA2,BETA3,GAMMA(1 ) ,GAMMA(2) ,GAMMA(3),PHI1 , 
1PHI2,PHI3,WT,QX,QY,QZ,UP1,UP2,UP3 

12  FORMAT ( 1  HI ,8X ,8HBETA1  =,F5.1,5H  DEG. ,14X,8HBETA2  =,F5.1,5H  DEG., 

1  14X.8HBETA3  =,F5.1 ,5H  DEG. .//9X.8HGAMMA1  =,F5. 1 ,5H  DEG. ,14X, 
28HGAMMA2  =,F5.1,5H  DEG. ,14X,8HGAMMA3  =,F5.1,5H  DEG.//9X,8HPHI1  = 
3.F5.1.5H  DEG..14X.8HPHI2  =,F5. 1 ,5H  DEG. .14X.8HPHI3  =,F5.1,5HD 

4EG.//9X.4HWT  =,F7.2,5X,4HQX  =,F7.2,5X ,4HQY  =,F7.2,5X,4HQZ  =,F7.2, 
55X.5HUP1  = , F7 . 2 , 5X , 5MUP2  =,F7.2,5X,5HUP3  =,F7.2) 

B(1 )  =  BETA1/(180. 0/3. 141593) 

B( 2)  =  BETA2/(180. 0/3. 141593) 

B( 3)  =  BETA3/{ 180.0/3. 141593) 

C(l)  =  GAMMA(1)/(180. 0/3. 141593) 

C(2)  =  GAMMA(2)/ (180.0/3. 141 593) 

C ( 3)  =  GAMMA(3)/(180.0/3. 141593) 

PHI1  =  PHI1/(180. 0/3. 141593) 

PHI2  =  PHI2/(180. 0/3. 141593) 

PHI 3  =  PHI3/(180. 0/3. 141593) 

DO  201  K=1 ,3 
UX(K)  =  C0S(B(K) ) 

UY(K)  =  SIN(B(K) ) 

UZ(K)  =  0.0 

VX(K)  =  COS(C(K) )*SIN(B(K)) 

VY(K)  *-COS(C(K))*COS(B(K)) 

VZ(K)  =-SIN(C(K)) 

WX(K)  =  (UY(K  *VZ(K))-(UZ(K)*VY(K)) 

WY(K)  =  (UZ(K)*VX(K))-(UX(K)*VZ(K)) 

WZ(K)  =  (UX(K)*VY(K) )-(UY(K)*VX(K)) 

IF<GAMMA(K) .GT.90.0)  GO  TO  201 
WX(K)  =  -WX(K) 

WY(K)  =  -WY(K) 

WZ(K)  =  -WZ(K) 

201  CONTINUE 

WRITE  (6,13)  UX( 1 ) ,UY( 1 ) ,UZ( 1 ) , VX ( 1 ) , VY( 1 ) , VZ( 1 ) ,WX( 1 ) ,WY ( 1 ) ,WZ( 1 ) 

13  FORMAT (1 HO, 9X,4HU1  =,F6.3,1H,,2X,F6.3,1H,,2X,F6.3/ 

1  10X.4HV1  =,F6.3,1H,,2X,F6.3,1H,,2X,F6.3/ 

2  10X.4HW1  =,F6.3,1H, ,2X,F6.3,1H, ,2X ,F6.3) 

WRITE  (6,14)  UX(2) ,UY (2) ,UZ(2) , VX( 2) , VY (2) ,VZ(2) ,WX(2),WY(2) ,WZ(2) 

14  FORMAT ( 1  HO , 9X , 4HU2  ■,F6.3,1H,,2X,F6.3,1H,,2X,F6.3/ 

1  10X.4HV2  *,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3/ 

2  10X.4HW2  =,F6.3,1H,,2X,F6.3,1H, ,2X,F6.3) 

WRITE( 6 ,115)  UX(3) ,UY (3) ,UZ( 3 ) ,VX(3) ,VY(3) ,VZ(3) ,WX(3) ,WY(3),WZ(3) 
115  F0RMAT(1H0,9X,4HU3  =,F6.3,1H, ,2X,F6.3,1H, .2X.F6.3/ 
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1  10X.4HV3  =,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3/ 

2  10X.4HW3  =,F6.3,1H, ,2X,F6.3,1H, ,2X,F6.3) 

W1X  =  WX(1) 

W1Y  =  WYM ) 

W1Z  »  WZ(1) 

W2X  =  MX  2 
W2Y  =  MY  2 
W2Z  =  WZ(2) 

W3X  =  MX  3 
W3Y  =  WY (3) 

W3Z  =  WZ(3) 

XI 2X=W2Y*W1 Z-W2Z*W1 Y 
X 1 2Y=W2Z*W1 X-W2X*W1 Z 
X 1 27=U2X*W1 Y-W2Y*W1 X 
XI 2=SQRT {XT  2X**2+X1 2Y**2+X1 2Z**2 ) 

X23X=W3Y*W2Z-M3Z*W2Y 
X23Y=W3Z*W2X-W3X*W2Z 
X23Z=W3X*W2Y-W3Y*W2X 
X23  =SQRT(X23X**2+X23Y**2+X23Z**2) 

X31 X=W1 Y*W3Z-W1 Z*W3Y 
X31Y=W1Z*W3X-W1X*W3Z 
X31 Z=W1 X*W3Y-W1 Y*W3X 
X31  =SQRT ( X31 X**2+X31 Y**2+X31 Z**2 ) 

WRITE(6,15)  X12X,X12Y,X12Z,X12,X23X,X23Y,X23Z,X23,X31X,X31Y,X31Z, 
1X31 

15  FORMAT ( 1  HO, 9X.5HX  12  =,2(F6.3,lH,,2X),F6.3,5X,10HABStX12)  =,F6.3// 

1  10X.5HX23  =,2(F6.3»1H, ,2X) ,F6.3,5X,10HABS(X23)  =,F6.3// 

2  10X.5HX31  =,2(F6.3,1H, ,2X),F6.3,5X,10HABS(X31 )  =,F6.3) 

RX=  QX  +  UP1*W1X  +  UP2*W2X  +  UP3*W3X 

RY=  QY  +  UP1*W1Y  +  UP2*W2Y  +  UP3*W3Y 

RZ=  QZ  +  UP1*W1Z  +  UP2*W2Z  +  UP3*W3Z  +  WT 

SI  21 X=  W1Y*X12Z-W1Z*X12Y 

SI  21 Y=  W1Z*X12X-W1X*X12Z 

SI  21 Z=  W1X*X12Y-W1Y*X12X 

SI 22X=  XI 2Y*W2Z-X1 2Z*W2Y 

S122Y=  XI 2Z*W2X-X1 2X*W2Z 

S122Z=  XI 2X*W2Y-X1 2Y*W2X 

S232X=  W2Y*X23Z-W2Z*X23Y 

S232Y=  W2Y*X23X-W2X*X23Z 

S232Z=  W2X*X23Y-W2Y*X23X 

S233X=  X23Y*W3Z-X23Z*W3Y 

S233Y=  X23Z*W3X-X23X*W3Z 

S233Z=  X23X*W3Y-X23Y*W3X 

S313X=  W3Y*X31Z-W3Z*X31Y 

S313Y=  W3Z*X31 X-W3X*X31 Z 

S313Z*  W3X*X31Y-W3Y*X31X 

S311X=  X31Y*W1Z-X31Z*W1Y 

S31 1 Y=  X31 Z*W1 X-X31 X*W1 Z 

S311Z=  X31X*W1Y-X31Y*W1X 

WRITE( 6 ,121 )  S121X,S121Y,S121Z,S122X,S122Y,S122Z,S232X,S232Y, 
1S232Z,S233X,S233Y,S233Z,S313X,S313Y,S31 3Z,S311X,S311Y,S311Z 
121  F0RMAT(1H0,9X,5HS121=,2(F6.3,1H,,2X},F6.3,5X,5HS122=,2CF6.3,1H,, 

12X) ,F6.3//10X,5HS232=,2(PG.3,1H, ,2X) ,F6.3,5X,5HS233=,2(F6.3,1H, , 
22X),F6.3//10X,5HS:i3-,2(F6.3,1H,,2X),F6.3,5X,5HS3n=,2(F613,1H,, 
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32X),F6.3) 

TEST1  =RX*X12X+RY*X12Y+RZ*X12Z 
TEST2  =RX*X23X+RY*X23Y+RZ*X23Z 
TEST3  =RX*X31X+RY*X31Y+RZ*X31Z 
TEST4  =RX*W1X+RY*W1Y+RZ*W1Z 
TEST5  =RX*W2X+RY*M2Y+RZ*W2Z 
TEST6  =RX*W3X+RY*W3Y+RZ*W3Z 
TEST7  =RX*S121X+RY*S121Y+RZ*S1 21 Z 
TEST8  =RX*S1 22X+RY*S1 22Y+RZ*S1 22Z 
TEST9  =RX*S232X+RY*S232Y+RZ*S232Z 
TEST10=RX*S233X+RY*S233Y+RZ*S233Z 
TEST1 1 =RX*S31 3X+RY*S31 3Y+RZ*S31 3Z 
TEST1 2=RX*S31 1 X+RY*S31 1 Y+RZ*S31 1 Z 

IF(TEST1. GE. 0.0. AND. TEST7.LE. 0.0. AND. TESTS. LE. 0.0)  GO  TO  101 
I F(TEST2. GE.O.O. AND. TEST9  .LE.O.O.AXD.TESTIO.LE.O.O)  GO  TO  102 
I F(TEST3. GE.O.O. AND. TEST1 1 .LE.O. 0.AND.TEST12.LE.0.0)  GO  TO  103 
IF(TEST4.LE.0.0.AND.TEST7  .GE. 0.0. AND. TEST12. GE.O.O)  GO  TO  104 
I F ( TEST5 .LE.O. 0 .AND . TEST8  .GE.0.0.AND.TEST9  .GE.O.O)  GO  TO  105 
I F ( TEST 6 . LE . 0 . 0 . AND . TEST! 0 . GE . 0 . 0 . AND . TEST1 1 . GE . 0 . 0 )  GO  TO  106 
IF(TEST4. GE. 0.0. AND. TEST5. GE.O.O. AND. TEST6. GE.O.O)  GO  TO  119 

101  WRITE ( 6,1 07 ) 

107  FORMAT ( 1  HO , 9X , 50HSL IDI NG  TENDS  TO  OCCUR  ALONG  X12  ON  PLANES  1  AND 
12) 

T12  =  TEST1/X12 
EN1 2X=RX-T1 2*X1 2X/X1 2 
EN1 2Y=RY-T1 2*X1 2Y/X1 2 
EN1 2Z=RZ-T1 2*X1 2Z/X1 2 
EEE  =W1X*W2Y-W2X*W1Y 
IF  (EEE. EQ. 0.0)  GO  TO  122 
EEE1 =W2X*EN1 2Y-W2Y*EN1 2X 
EEE2=W1 Y*EN1 2X-W1 X*EN1 2Y 
EN1  =EEE1/EEE 
EN2  =EEE2/EEE 
GO  TO  123 

122  EEE=W1Y*W2Z-W1Z*W2Y 
EEE1=W2Y*EN12Z-W2Z*EN12Y 
EEE2=W1 Z*EN1 2Y -W1 Y*EN1 2Z 
EN1=EEE1/EEE 
EN2=EEE2/EEE 

123  FS  =((EN1*SIN(PHI1 ) /COS ( PHI  1 ))+(EN2*SlN(PHI2)/C0S(PHI2)))/T12 
WRITE(6,108)  FS,EN1 ,EN2,T12 

108  F0RMAT(1H0,9X,4HFS  =,F5.2,10X,5HEN1  =,F7.2,5X,5HEN2  =,F7.2,5X,5FT1 
12  =,F7. 2) 

GO  TO  200 

102  WR I TE( 6 , 1 09 ) 

109  F0RMAT(1H0,9X,50HSLIDING  TENDS  TO  OCCUR  ALONG  X23  ON  PLANES  2  AND 
13) 

T23  =  TEST2/X23 
EN23X=RX-T23*X23X/X23 
EN23Y=RY-T23*X23Y/X23 
EN23Z=RZ-T23*X23Z/X23 
EEE  =W2X*W3Y-W3X*W2Y 
IF(EEE.EQ.O.O)  GO  TO  124 
EEE2  =W3X*EN23Y-W3Y*EN23X 
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EEE3  =W2Y*EN23X-W2X*EN23Y 
EN2  =EEE2/EEE 
EN3  =EEE3/EEE 
GO  TO  125 

124  EEE=W2Y*W3Z-W2Z*W3Y 
EEE2=W3Y*EN23Z-W3Z*EN23Y 
EEE3=W2Z*EN23Y-W2Y*£N23Z 
EN2=EEE2/EEE 
EN3=EEE3/EEE 

1 25  FS  =( (EN2*SIN(PHI2)/COS(PHI2) )+(EN3*SIN(PHI3)/C0S(PHI3) ) )/T23 
WRITE(6, 110)  FS,EN2,EN3,T23 

110  F0RMAT(1H0,9X,4HFS  =,F5.2,10X,5HEN2  =,F7.2,5X,5HEN3  =,F7.2,5X,5HT2 
13  =,F7.2) 

GO  TO  200 

103  WRITE{ 6 ,111) 

111  F0RMAT(1H0,9X,50HSLIDING  TENDS  TO  OCCUR  ALONG  X31  ON  PLANES  3  AND 
11) 

T31  =  TEST3/X31 
EN31 X=RX-T31 *X31 X/X31 
EN31Y=RY-T31*X31Y/X31 
EN31 Z=RZ-T31 *X31 Z/X31 
EEE  =W3X*W1Y-W1X*W3Y 
IF(EEE.EO.O.O)  GO  TO  126 
EEE3  =W1X*EN31Y-W1Y*EN31X 
EEE1  =W3Y*EN31 X-W3X*EN31 Y 
EN3  =EEE3/EEE 
EN1 =EEE1 /EEE 
GO  TO  127 

126  EEE=W3Y*W1 Z-W3Z*W1 Y 
EEE3=W1 Y*EN31 Z-Wl Z*EN31 Y 
EEE1 =W3Z*EN31Y-W3Y*EN31 Z 
EN3=EEE3/EEE 
EN1=EEE1/EEE 

127  FS  =((EN3*SIN(PHI3)/C0S(PHI3))+(EN1*SIN(PHI1 )/COS(PHIl )))/T31 
WRITE(6,112)  FS,EN3,EN1,T31 

112  FORMAT ( 1  HO , 9X , 4HFS  =,F5.2,10X,5HEN3  =,F7.2,5X,5HEN1  =,F7.2,5X,5HT3 
11  =,F7.2) 

GO  TO  200 

104  WRITE(6,113) 

113  F0RMAT(1H0,9X,33HSLIDING  TENDS  TO  OCCUR  ON  PLANE  1) 

EN1  =  -TEST4 

T1X  =  RX  +  EN1*W1X 
T1Y  =  RY  +  EN1*W1Y 
T1Z  =  RZ  +  EN1*W1Z 
T1  =  SQRT(T1 X**2+T1 Y**2+T1 Z**2 ) 

FS  =  EN1*SIN(PHI1 )/(COS(PHIl )*T1 ) 

WRITE( 6 ,114)  FS,EN1,T1 

114  FORMAT (1 HO, 9X.4HFS  *,F5.2,10X,5HEN1  =,F7.2,5X,5HT1  =,F7.2) 

GO  TO  300 

105  WRITE(6,215) 

215  FORMAT { 1  H<7, 9X , 33HSL I DI NG  TENDS  TO  OCCUR  ON  PLANE  2) 

EN2  =  -TESTS 
T2X  =  RX  +  EN2*W2X 
T2Y  =  RY  +  EN2*W2Y 
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T2Z  =  RZ  +  EN2*W2Z 

T2  =  SQRT(T2X**2+T2Y**2+T2Z**2) 

FS  *  EN2*SIN(PHI2)/(C0S(PH12)*T2) 

WRITE(6,116)  FS.EN2.T2 

116  F0RMAT(1H0,9X,4HFS  =,F5.2,10X,5HEN2  =,F7.2,5X,5HT2  =,P7.2) 

GO  TO  300 

106  WRITE(6,117) 

117  FORMAT (1H0,9X,33HSL ID ING  TENDS  TO  OCCUR  ON  PLANE  3) 

EN3  =  -TEST6 

T3X  =  RX  +  EN3*W3X 
T3Y  =  RY  +  EN3*W3Y 
T3Z  =  RZ  +  EN3*W3Z 
T3  =  SQRT(T3X**2+T3Y**2+T3Z**2) 

FS  =  EN3*SIN(PHI3)/(C0S(PHI3)*T3) 

WRITE(6,118)  FS.EN3.T3 

118  FORMAT ( 1 HO , 9X , 4HFS  =,F5.2,10X,5HEN3  =,F7.2,5X,5HT3  =,F7.2) 

GO  TO  300 

119  WRITE(6,120) 

120  FORMAT ( 1  HO , 9X , 40HR0C  K  WEDGE  IS  LIFTED  OFF  THE  BASE  PLANES) 

300  IF(TEST1 .LT.O.O.AND.TEST?  .LE.O.O.AND.TEST8  .LE.O.O)  GO  TO  101 
IF(TEST2.LT.0.0.AN0.TEST9  . LE.O.O. AND. TEST10. LE.O.O)  GO  TO  102 
IF(TEST3. LT. 0.0. AND. TEST!!. LE.O.O. AND. TEST12.LF. 0.0)  GO  TO  103 

200  CONTINUE 
STOP 
END 


*  DATA 
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In  accordance  with  letter  from  DAEX-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Hendron,  Alfred  Joseph 

Analytical  and  graphical  methods  for  the  analysis  of 
slopes  in  rock  masses  /  by  A.  J.  Hendron,  Jr.,  E.  J. 

Cording,  and  A.  K.  Aiver,  Department  of  Civil  Engineering, 
University  of  Illinois,  Urbana,  Ill.  Vicksburg,  Miss.  : 

U.  S.  Waterways  Experiment  Station  ;  Springfield,  Va.  : 
available  from  National  Technical  Information  Service,  1980. 

xi,  148,  lL  p.  :  ill.  ;  27  cm.  (Technical  report  -  U.  S. 

Army  Engineer  Waterways  Experiment  Station  ;  GL-80-2) 

Prepared  for  U.  S.  Army  Engineer  Nuclear  Cratering  Group, 
Livermore,  Calif.,  under  Contract  No.  DACW  39-67 -C-0097. 

Reprint  of  NCG  Technical  Report  No.  36. 

References:  p.  llt8. 

1.  Dynamic  slope  stability.  2.  Graphical  methods.  3.  Rock 
masses,  h.  Rock  mechanics.  5-  Slope  stability.  6.  Stereonet. 
7.  Vector  analysis.  I.  Aiyer,  Arunachalam  Kulathu,  joint 
author.  II.  Cording,  Edward  J.,  joint  author.  III.  Illinois. 

(Continued  on  next  card) 
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University.  Dept,  of  Civil  Engineering.  IV.  United  States. 
Army  Engineer  Nuclear  Cratering  Group.  V.  Series:  United 
States.  Waterways  Experiment  Station,  Vicksburg,  Miss. 
Technical  report  ;  GL-80-2. 
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